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ABSTRACT 

We present results from an adaptive optics survey for substellar and stellar companions to Sun-like 
stars. The survey targeted 266 F5-K5 stars in the 3 Myr to 3 Gyr age range with distances of 10- 
190 pc. Results from the survey include the discovery of two brown dwarf companions (HD 49197B 
and HD 203030B), 24 new stellar binaries, and a triple system. We infer that the frequency of 0.012- 
0.072 M© brown dwarfs in 28-1590 AU orbits around young solar analogs is 3. 2~^\% (2a limits). 
The result demonstrates that the deficiency of substellar companions at wide orbital separations from 
Sun- like stars is less pronounced than in the radial velocity "brown dwarf desert." We infer that the 
mass distribution of companions in 28-1590 AU orbits around solar-mass stars follows a continuous 
dN/dM 2 oc M^ QA relation over the 0.01-1.0 M@ secondary mass range. While this functional form 
is similar to the that for <0.1 Mq isolated objects, over the entire 0.01-1.0 Mq range the mass 
functions of companions and of isolated objects differ significantly. Based on this conclusion and 
on similar results from other direct imaging and radial velocity companion surveys in the literature, 
we argue that the companion mass function follows the same universal form over the entire range 
between 0-1590 AU in orbital semi-major axis and ~0. 01-20 Mq in companion mass. In this context, 
the relative dearth of substellar versus stellar secondaries at all orbital separations arises naturally 
from the inferred form of the companion mass function. 

Subject headings: stars: binaries: visual — stars: imaging — stars: low-mass, brown dwarfs — stars: mass 
function 



1. INTRODUCTION 

The properties of brown dwarf companions to stars are 
important for understanding the substellar companion 
mass function (CMF), the formation of brown dwarfs, 
and the formation and evolution of low-mass ratio binary 
systems. Widely-separated brown dwarf companions, in 
particular, are an important benchmark for studying the 
properties of substellar objects because of their accessi- 
bility to direct spectroscopic characterization and their 
relative ease of age-dating — from assumed co-evality with 
their host stars. 

However, brown dwarf companions have been an elu- 
sive target for direct imaging. The main challenge has 
been the need to attain sufficient imaging contrast to 
detect secondaries that are >10 3 fainter than their host 
stars at angular separations spanning solar system-like 
scales (<40 AU = 074 at 100 pc). 

The problem is alleviated at young ages when brown 
dwarfs are brighter. In addition, nearby stars offer an 
additional advantage because the relevant angular scales 
are correspondingly wider and more accessible to direct 
imaging. Young nearby stars are thus the preferred tar- 
gets for substellar companion searches through direct 
imaging. 

Nevertheless, early surveys for substellar compan- 
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ions, performed with seeing-limited or first-generation 
high- c ontrast imaging technology ( Oppcnheimcr et all 
[20011: iHinz et all 12002 : IMcCarthv fc Zuckermanl l2004h 
had very low detection rates, < 0.5%. This low brown 
dwarf companion detection rate was similar to that in- 
ferred from pr ecision radial velocity surveys (< 0.5% 
over 0-3 AU: | Marcv fc Butlerl I2000T ) . and prompted 
IMcCarthv fc Zuckermanl (|2004l ) to conclude that the so- 
called "brown dwarf desert" extends far beyond the or- 
bital separations probed by radial velocity surveys, out 
to at least «1200 AU. 

Over the past few years, advances in adaptive op- 
tics (AO) technology and high-contrast imaging meth- 
ods have improved the chances for the direct imaging 
of substellar secondaries. Modern AO systems, with 
>200 corrective elements spread across the beam of a 5- 
10 m telescope, are able to deliver high order rectification 
(<250 nm r.m.s. residual error) of wavefronts perturbed 
by Earth's turbulent atmosphere at up to kHz rates. In 
addition, our empirical appreciation of the local young 
stellar population has improved over the past decade, as 
demonstrated by the recent discoveries of a large number 
of young (< 500 M yr) stellar associations within 200 pc 
from the Sun fe.g.. iKastner et al.lll997l: iMamaiek et alj 
19991: IZuckerman fc Webb! 120001: IZuckerman et al.ll2001l : 
Montes et al7 l2001: Zu ckerman & Songll2004l . and refer- 
ences therein). These have allowed us to select more 
suitable targets for direct imaging searches for substellar 
companions. 

Several recent direct imaging surveys of nearby 
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young stellar associ ations conducted on high-orde r AO- 
cquippe d telescopes (iNeuhauser fc Guentherll2004l . 25 A- 
M stars; IChauvin et al.ll2005allbl , 50 A-M stars) or with 
the Hubble Space Telescope (HST: [Lowrance et al.ll2005l 
45 A-M stars) have enjoyed higher detection rates (2- 
4 %) than the first generation of surveys. In addition, 
at very wide (>1000 AU) separations, where the detec- 
tion of brown dwarf companion s to solar-neighbo rhood 
stars is not hindered by contrast. iGizis et' al. (2001) have 
found that the frequency of substellar companions to 
F-K dwarfs is fully consistent with that of stell ar com- 
panions to G dwarfs (jDuauennov fc Mavorlll99iT ). Thus, 
while the radial velocity "brown dwarf desert" remains 
nearly void within 3 AU even after the discovery of nu- 
merous extra-solar planets over the past decade, brown 
dwarf secondaries at > 100-1000 AU separations seem to 
not be as rare. 

The precise frequency of substellar companions in 
direct imaging surveys remains controversial. Sev- 
eral highly sensitive s urveys performed with the HST 
(ISchroeder et all I2000L 23 A-M stars; i Brandner et "all 
120001 28 G-M stars: iLuhman et al.ll2005l 150 B- M stars) 
and with high-order AO (|Masciadri et al.ll2005l 28 G-M 
stars; IBiller et al.ll2007l 54 A-M stars) have reported null 
detections of substellar companions. Given the low (few 
per cent) detection rate of substellar companions in di- 
rect imaging surveys, it is now clear that, with < 50 
targets per sample, some of these surveys were too small 
to expect to detect even a single brown dwarf compan- 
ion. However, the non- detection of substellar secondaries 
in tw o largest surveys (|Luhman et alj l2005t IBiller et al.l 
l2007f ) is potentially significant. 

Given current understanding of the importance of stel- 
lar mass for (1) stellar multiplicity rates (see review in 
ISterzik fc D urisenl |2004f) and (2) binary mass ratio dis- 
tributions (see review in lBurgasser e t al. 2007), it is im- 
perative that any study of the substellar companion fre- 
quency is considered in the context of the mass distribu- 
tion of primary stars in the sample. Indeed, a large sur- 
vey sample comprising primaries with identical masses is 
ideal. 

The problem of the brown dwarf companion frequency 
is perhaps most comprehensively dealt with in the con- 
text of solar mass primaries. For these a uniquely large 
body of stellar and substellar multiplicity data exist on 
all orbital scales. On one hand, the exhaustive spec- 
troscopic and imaging stud y of G dwarf multiples by 
iDuquennov fc Mavorl (fl991) provides an important an- 
chor to the properties of 0.1-1.0 M© stellar companions 
to Sun- like stars. On the other hand, the results from 
more than a decade of precision radial velocity surveys 
for planets around G and K stars allow a comparison 
with the planetary-mass end of the substellar companion 
mass range. 

A large uniform sample of young Sun-like stars has 
been compiled by the Formation and Evolution of Plan- 
etary Systems (FEPS) Spitzer Legacy team. The pur- 
pose of the FEPS Legacy campaign with Spitzer was to 
study circumstellar disk evolution in the mid-IR. How- 
ever, the sample is also well-suited for a high-contrast 
imaging survey for substellar companions. Seventy per- 
cent of the FEPS stars are younger than ~500 Myr, and 
all are within 200 pc. 

As an auxiliary component to the FEPS program, we 



imaged most of the northern FEPS sample with the high- 
order AO systems on the Palomar 5 m and the Keck 
10 m telescopes. We further expanded our AO survey 
by observing several dozen additional nearby and mostly 
young solar analo gs. Preliminary results from the project 
were published in IMetchev fc Hillen brand! (|2004l ) and in 
IMetchev fc Hillenbrand! ( 2006f ). including the discover- 
ies of two brown dwarf companions: HD 49197B and 
HD 203030B. The survey has now been completed, and 
no further brown dwarf com panions have b een found. 
The results were analyzed in IMetchev! ([2006). Here we 
present the AO survey in its entirety and focus on the 
statistical interpretation of the data. 

The paper is organized as follows. A full description of 
the survey sample is given in § [2j The Palomar and Keck 
AO observing campaigns and the data reduction and cal- 
ibration techniques are described in § [3] The candidate 
companion detection approach and the survey detection 
limits are discussed in § |H The various methods used 
for bona fide companion confirmation are presented in 
§ Section [5] summarizes the results from our survey, 
including all of the newly-discovered and confirmed sub- 
stellar and stellar secondaries. Section [7] contains a brief 
discussion of the various sources of incompleteness and 
a full discussion of the biases in the survey. (A full- 
fledged incompleteness analysis is presented in the Ap- 
pendix.) In § [8] we estimate the frequency of wide sub- 
stellar companions to young solar analogs, and present 
evidence for trends in the companion mass and compan- 
ion frequency with semi-major axis and primary mass. 
In §[5] we consider the results of the current investigation 
in the broader context of stellar multiplicity, and suggest 
the existence of a universal CMF. Section ITOl summarizes 
the findings from our study. 

2. TARGET SAMPLE 

The main criteria used for selecting stars for the survey 
were Sun-like mass, youth, proximity, and visibility from 
the Northern hemisphere. In this Section we describe 
how they were applied to generate our AO survey sample. 

2.1. Selection 

The selection of the AO survey sample was largely 
based on the target selection criteria empl oyed in the 
const ruction of the FEPS program sample (jMever et al.l 
2006). The FEPS selection criteria will not be repro- 
duced here. The final FEPS target list comprises 328 
F5-K5 stars within 200 pc distributed uniformly in loga- 
rithmic age intervals between 3 Myr and 3 Gyr. Approx- 
imately a third of these are members of open clusters and 
stellar associations, and the remainder are field stars. We 
observed 228 of the 240 FEPS stars north of S = -30° 
with AO at Palomar or Keck. 

A further 38 solar analogs were added to the AO sur- 
vey toward the end of the first epoch of observations to 
bolster the sample size, mirroring FEPS target-selection 
policy. The additional stars were selected from three 
sources: (1) the broader compilation of FEPS candi- 
date targets, including stars that had been eliminated 
from the final FEPS sample based on infrared back- 
ground or age redundancy conside rations; (2) th e com - 
pila tions of nearby young st ars bv iMontes et al.l (|2001l ) 
and IWichmann et alj (|2003l ); a nd (3) our own P alomar 
echelle survey of nearby stars (|White et al.l [2007). The 
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final set of 266 targets in our AO sample has spectral 
type and age distributions similar to those of the FEPS 
sample. 

2.1.1. Deep and Shallow Samples 

To optimize sensitivity to substellar companions, we 
chose to observe a portion of the youngest and nearest 
AO sample stars with deep coronagraphic exposures. We 
applied the following additional guidelines to select stars 
for the deep coronagraphic sub-sample: 

1. stellar activity and lithium levels indicating ages 
less than 500 Myr; 

2. no AK S < 4 objects between 0'.'8 and 13'.'0, as 
determined from real-time flux ratio measurements 
during survey observations; 

3. nearby stars were given priority over more distant 
stars; 

4. isolated stars, not belonging to one of the young 
open clusters or stellar assocciations, were given 
priority for deep observations. 

The first criterion was motivated by the fact that sub- 
stellar companions should be intrinsically brightest at 
the youngest ages. The second constraint was aimed at 
avoiding the loss of sensitivity to faint objects over a 
large portion of the detector field of view (FOV) because 
of the presence of a bright neighboring star. 1 Binaries 
with separations < 0.8" had both their components suf- 
ficiently well-covered by the 1" coronagraphs in the Palo- 
mar and Keck AO cameras that they were allowed in the 
deep sample. The motivation for the third constraint 
was to optimize sensitivity to substellar companions at 
the smallest physical separations. The last criterion was 
applied to avoid duplication with previous sensitive high- 
angula r resolution studies of open clusters: iBouvier et al.1 
(119971 the Pleiades, A O). iKohler et all (Up per Scor- 
pius: 12000 . speckle), and lPatience et all (|2002t a Persei, 
speckle). 

Based on the additional criteria outlined above, 84 of 
the 228 stars selected from the final FEPS sample and 
16 of the 38 additional targets were included in our deep 
sample. The deep sub-sample thus consists of 100 F5-K5 
stars with ages less than 500 Myr. 

All remaining stars were observed primarily in short se- 
quences of non-coronagraphic images to establish stellar 
multiplicity. These will be referred to as the "shallow" 
sample. The shallow sample includes 11 stars older than 
500 Myr that were also observed with long coronagraphic 
exposures: 2 Hyades (~600 Myr) members and 9 other 
stars whose subsequent age-dating showed that they were 
older than originally estimated. Although these 11 stars 
were observed coronagraphically, for the purpose of lim- 
iting our deep sample only to the observations with the 
highest sensitivity to substellar mass, they are not con- 
sidered as part of the deep sample. 

1 Following more accurate post-reduction photometry, a 3'.'1 
companion to one of the stars in our deep sample, HD 31950, was 
found to be only AKg = 3.70 mag fainter (Table|8]l. Although this 
companion violates criterion [2] we have chosen to keep HD 31950 
as a member of the deep sample. 



The deep and shallow sample stars and their character- 
istics are listed in Tables [T] and [2l Median age, distance, 
and spectral-type statistics for the deep, shallow, and 
complete ( deep + shallow) samples are given in Table |3l 

2.2. Stellar Properties 

Our sample stars are near-solar (G2 V) analogs, rang- 
ing in spectral type between F5 and K5 (6300 K > 
T ff > 4400 K) and, depending on stellar age, be- 
tween IV and V in luminosity class (3.4 < log<? < 
4.5 in cgs units). The corresponding mass range, 
based on dynamical mass estimates in bi nary systems 
and on stellar thermodynamic models (jBaraffe et al.l 
ll998HD'Antona fc Mazziteliilll994f l is approximately 0.7- 
1.3 Mo, fol l owing the design of the FEPS sample 
(|Meyer et al.l 12006). For greater detail in the assign- 
ment of spectral types, effective temperatures, and sur- 
face gravities to FEPS sample stars we refer the reader to 
ICarpenter et alj ((2008). Histograms of the distribution 
of stellar effective temperatures and masses of all stars 
in our AO survey sample are shown in Figure [T] 

Seventy-nine of our sample stars are members of 
known young stellar associations: Upper Scorpius, 
a Persei, the Pleiades, and the Hyades. For these we 
have adopted ensemble ages from the lit erature: 5 Myr 
for th e Upper Scorpius OB assoc iation (iPreibisch et ail 
120021 1 . 80 Myr for a Per sei (IVentura et al.l 1 199811 . 
120 Myr for the Pleia des (IVentura et all 1 19981) . and 
600 Myr for the Hyades (jPerryman et al.lll998l ). Sample 
stars that do not belong to any known associations 
were age-dated following one of two approaches: (1) 
based on the strength of the chromosphcric Ca II H 
and K (3968A and 3933A) line emission for > 30 Myr- 
old stars, and using the r e cent activity-age relation 
of iMamaiek fc Hillenbrand] (|2008H . or (2) isochrone 
fitting for < 30 Myr-old pr e -main sequence stars using 
the tracks of [Baraffe et al.l (|1998| l . Where data from 
high-dispersion optical spectra were previou sly available 
(jStrassmeier et al.l 120001: IWhite et all 120071) . these were 
also reviewed for the strength of the lithium 6708 A 
absorption line to put additional constraints on the 
ages. All chromospheric and isochronal ages were 
assumed uncertain to within a factor of two. Ages 
for a remaining set of 20 stars not present in the 
exten ded FEPS sample were taken f rom the litera- 
ture dBarrado v Navascues et ail Il997t iGutierrez et al.1 
1999t iMontes et al l 120011 : IWichmann et all 120031 : 
Nordstrom et al.l 12004( 1. Histograms of the age distri- 
bution of the complete survey sample and of the deep 
sub-sample are shown in Figure [2] 

Distances to 166 sample stars with individual Hippar- 
cos parallaxes were ta ken from the Hipparcos catalog 
(jPerryman et all 11997( 1 . For a further 55 known mem- 
bers of young open clusters and OB associations, we 
adopted the corresponding mean cluster distance, calcu- 
lated from a combination of trigonometric, orbital, sec- 
ular, and cluster parallaxes in the literature, a s inferred 
from Hipparcos and Tycho-2 (|Hog et al.l 12000) astrome- 
try, long-baseline intcrfcromctry, or high-resolution spec- 
troscopy. The adopted distances for open cluster mem- 
bers were: 133 ±6 pc for stars in the Pleiades (a weighted 
mean of the distances to seven members presented in 
Pan et al.ll2004 iMunari et al.ll2004 TZwahlen et al]|2004 
Southworth et al.ll2005l andassuming ~ 1° cluster angu- 
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lar extent from I Adams et ai1l200l!). and 190 ± 11 pc for 
stars in a Persei ( van Leeuwenll 19991 assuming 1° cluster 
radius). For stars belonging to t he Upper Scorpius as- 
sociation we adopted 145 ± 40 pc (|de Zeeuw et alJH999t 
iMamajek et al.l r2002T ) . All of these distances agree with 
estimates from main-sequence fitting for the correspond- 
ing clusters. For 18 more stars we ado pted secular par- 
allax es from lMamaiek et all ((20021) and lMamaiekl (|2004 
12007ft - Finally, for 27 remaining >30 Myr-old stars we 
obtained approximate distances based on a combination 
of moving group association, secular parallax, and spec- 
troscopic parallax, with care to avoid redundancy in dis- 
tance and age derivation. The errors on the distances in 
these cases were conservatively assumed to be 25%-50%. 
More refined distance and age estimates for these stars 
will be included in a future publication from the FEPS 
program. 

Accurate proper motions for the sample stars are essen- 
tial in identifying bona fide companions through multi- 
epoch astrometry. Proper motions for the 166 stars with 
individual Hipparcos parallaxes were taken from the Hip- 
parcos database. For the remaining 100 stars proper mo- 
tions were adopted from The Second U.S. Naval Observa- 
tory CCD Astrograph Catalog (TJCAC2; Zacharia s et al.l 
120041) and from the Tycho-2 Catalog (|Hpg et alll2000f ). 
The three catalogs provided similar astrometric accu- 
racy (±1.0 mas yr _1 ) for the sample stars, although the 
UCAC2 and Tycho-2 catalogs went deeper. 

Figure [3] presents histograms of the heliocen- 
tric distances (panel a) and total proper motions 
(y/ (fi a cosS) 2 + fij; panel b) of the stars in the complete 
sample and in the deep sub-sample. The bi-modal distri- 
bution of the distances is a combined effect of the large 
heliocentric distances (130-190 pc) of the youngest (3- 
120 Myr) stars in the sample, and of the preference given 
to closer systems at older ages. 

3. OBSERVATIONS, DATA REDUCTION, AND 
CALIBRATION 

3.1. Observing Strategy 

A complete description of the observing strategy of our 
AO su rvey was given in iMetchev fc Hillenbrand! ((2004) 
and in lMetchevI J2006). Here we briefly review the ap- 
proach and summarize the survey observations. 

The full sample of 266 stars was observed in the near- 
IR with AO at the Palomar and Keck II telescopes on 
47 clear nights over the course of 3 years: between 2002 
January 31 and 2005 January 24. Additional astrometric 
follow-up was obtained during 2006 and 2007 in a few 
individual cases. 

We opted to conduct the entire survey in the K$ 
band to take advantage of the much better AO perfor- 
mance at >2/i m. Although cool T-type brown dwarfs 
(T e s < 1400 K; lGolimowski et al.l[200l iVrba et al.ll200"l 
are faintest at K band in the near-IR, warmer (poten- 
tially younger) L-type brown dwarfs are brightest at K. 
Thus, given superior imaging contrast and the relative 
youth of our deep sample, the 2/^m region was seen as the 
best choice for optimizing sensitivity to close-in young 
substellar companions. 

The majority of the science targets were observed 
first at Palomar. Only seven of the targets (all be- 
longing to the deep sample) were observed initially and 
only at Keck. The Palomar campaign was conducted 



with the PALAO system dTrov et alj |2000D a nd the 
PHARO near-IR camera (jHavward et al.l |200H ) in its 
25 mas pix -1 mode, providing a 25.6" x 25.6" FOV. At 
Keck , we used the facility AO system (|Wizinowich et ah! 
2000) on Keck II and the NIRC2 near-IR camera in its 
40 mas pix" 1 mode, offering an FOV of 40.6" x 40.6". To 
improve overall sensitivity and contrast the 100 targets 
in the deep sample were observed coronagraphically with 
the opaque 0'.'97-diameter occulting spot in PHARO and 
the partially transmissive l'.'O-diameter occulting spot in 
NIRC2. All of the sample stars were sufficiently bright 
to allow use of the AO systems in NGS mode, i.e., to 
have the wavefront sensing performed on the primaries 
themselves. 

3.1.1. First-Epoch Observations at Palomar 

We spent 24 min of net exposure time per target dur- 
ing first epoch at Palomar, attaining an imaging depth 
of K$ ~ 19.7 mag on average for stars in the deep sam- 
ple. The 24 min of exposure were divided in 4 sets of 
6 min taken at different orientations of the telescope 
Cassegrain ring (CR) rotator. The 6 min of net expo- 
sure at each CR rotator angle consisted of two sets of 
three one-minute on-target exposures, with three one- 
minute sky exposures in between. The purpose of the 
CR rotation approach was to improve the quality of 
point-spread function (PSF) subtractions for data taken 
with an equatorial- mount telescope (Palomar) , in a man- 
ner similar to that attained with angular differential 
imaging (API) on a ltitude-azimuth-mounted telescopes 
(jMarois et al.l 120061 ). Stacking images taken at differ- 
ent CR angles also averages out detector and sky noise, 
much like mosaicking dithered images. Unfortunately, 
in addition to being less efficient, the CR rotation ap- 
proach was later found to also produce notable smear- 
ing of the PSF in the co-added de-rotated images at 
> 5" from the star, leading to degradati on both in imag - 
ing depth and in astrometric precision (jMetchev|[2006). 
We have since demonstrated that judicious matching of 
nearby science targets to use as PSFs for one another — a 
suitable approach for surveys of target-rich young stel- 
lar associations — enhances the contrast attainable with 
PALAO by 0.5-1.0 mag over the one repor ted here with- 
out in curring the overhead of CR rotations ((Tanner et al.l 
120071 

We used two different undersized Lyot stops to block 
the secondary obscuration and the secondary mirror sup- 
port structure at Palomar: the "medium" and the "big" 
cross, obscuring 40 % and 76% of the tota l telescope aper- 
ture, respectively (|Havward et al.l (2001!) . The use of an 
appropriately sized Lyot stop was expected to notice- 
ably improve the dynamic range achievable in high-order 
AO coronagraphy by sup pressing light diffracted by th e 
edge of the coronagraph (jSivaramakrishnan et alj [2001). 
Early experiments with th e PALAO/PHARO system by 
lOppenheimer et ail (|2000f ) had suggested that the big 
cross provided the best contrast in single exposures of 
up to several seconds, outperforming the medium and 
"standard" (no undersizing) Lyot masks by up to 0.5 mag 
between 0.5"-2.0" from bright stars. However, our expe- 
rience from observing each star in multiple longer expo- 
sures was that the less oversized Lyot stops allowed bet- 
ter real-time monitoring of the star-coronagraph align- 
ment and more accurate post-processing image registra- 
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tion and astrometry. With the medium and the standard 
Lyot stops the position of the star behind the corona- 
graph could be monitored by the location of a Poisson- 
like spot within the dark area of the coronagraph: the re- 
sult of constructive interference of high spatial frequency 
light pushed by the coronagraph to the periphery of the 
Lyot plane. The big Lyot stop likely shutters incom- 
ing starlight too aggressively to allow the formation of 
a sufficiently bright Poisson spot. Because image regis- 
tration of multiple exposures was crucial for obtaining 
greater overall exposure depth, we stopped using the big 
cross after March 2002. Given the adopted technique of 
rotating the Cassegrain ring to four mutually orthogonal 
orientations during the imaging of each star, the medium 
cross provided the best compromise between registration 
ability for the final images and consistency with which 
it would obscure the telescope secondary mirror support 
at each CR orientation. At the end of our survey, only 
seven of the 100 stars in the deep sample had their deep- 
est observations obtained with the big Lyot cross setup. 

In addition to the long coronagraphic iCg-band expo- 
sures, we also observed each deep sample target in short 
(1.4-10 s) unocculted exposures. These were taken to 
check for stellar multiplicity within the 0'.'5 coronagraph 
radius and to allow relative photometric calibration of 
the deep occulted exposures. The short exposures were 
obtained at J, H, and Ks bands, using a 1% neutral 
density (ND) filter to prevent saturation whenever nec- 
essary. For these observations we used a five-point dither 
pattern as is standard for infrared imaging. 

The 166 targets in the shallow sample were observed 
only in short dithered JHKs exposures at Palomar to 
check for stellar multiplicity. The imaging depth of the 
shallow survey varied greatly from star to star, depend- 
ing on whether the ND filter was used or not, and was 
generally in the 12 < Ks < 17 mag range. 

3.1.2. Follow-up Observations 

After an examination of the initial Palomar images, 
target stars which contained other objects in the same 
image — candidate companions — were followed up with 
additional imaging at later epochs to test for common 
proper motion between the candidate companions and 
the host stars. Upon establishing common proper mo- 
tion, candidate companions were observed spectroscopi- 
cally to confirm their physical association with the pri- 
mary. 

The imaging and spectroscopic follow-up was per- 
formed at both Palomar and Keck. Imaging at Keck was 
done in 6 x 1 min coronagraphic integrations per target, 
with an additional 3x1 min spent on sky. We used the 
"inscribed circle" NIRC2 pupil mask (90.7% throughput) 
to occult the telescope mirror outer edge. (None of the 
available NIRC2 pupil masks occult the Keck segment 
edges and the secondary support structure.) In most 
cases the 6 min-long exposures at Keck were «0.5 mag 
deeper (Ks ~ 20.2 mag) than the 24 min Palomar expo- 
sures, and occasionally revealed new candidate compan- 
ions. Nearly half (48/100) of our deep sample stars were 
observed at Keck in addition to at Palomar, including 
the seven targets observed only at Keck. Because of the 
marginal difference in the depths of the Keck and Palo- 
mar components of the deep survey, and for the sake of 
preserving the integrity of our well-defined 100-star deep 



sample, we analyze the Palomar and Keck AO campaigns 
together as a single survey. 

We obtained near-IR spectra of several bona fide 
and candidate companions for the purposes of fur- 
ther confirmation of their physical association and 
characterization of their photospheres. The spectro- 
scopic observations and data reductio n were described in 
IMetchev fc Hillenbrand! (|2004l [2006). Here we present 
spectroscopy of only one additional companion candi- 
date, to ScoPMS 214. The observations and data re- 
duction for that are briefly described in § 16.3.11 

3.2. Imaging Data Reduction 

The imaging data reduction procedure for the survey, 
including flat-fielding, sky-subtraction, bad-pixel correc- 
ti on, image registration, and imag e stacking was detailed 
in IMetchev &: H illenbrand (2004j). We have since aug- 
mented the procedure to include a correction for the non- 
linear flux response of the PHARO and NIRC2 detectors. 
Near-infrared detectors often have non-linear response 
even at small flux levels, that is important to take into ac- 
count when seeking accurate photometry. We measured 
the non-linearity of the PHARO and NIRC2 detectors 
from series of variable-length exposures of the brightly 
illuminated telescope dome interiors, interspersed with 
multiple dark frames to mitigate charge persistence ef- 
fects. The response of the PHARO camera, which em- 
ploys an HgCdTe detector that does not support multiple 
non-destructive reads, was found to be >1% non-linear 
beyond 10,000 counts/pix and >5% non-linear beyond 
45,000 counts/pix. The InSb detector on NIRC2, which 
supports non-destrictive read-outs, was found to be >1% 
non- linear beyond 3000 counts/pix/read and >5% non- 
linear beyond 7000 counts/pix/read. We created custom 
IDL routines 2 to linearize the PHARO and NIRC2 flux 
response. The linearization was applied to all images 
before any other data reduction steps. 

To enhance our ability to detect faint candidate com- 
panions in the deep coronagraphic exposures we at- 
tempted various methods of PSF removal, including: (1) 
subtracting a median-combined PSF of the star formed 
from the individual images taken at all four CR angles 
at Palomar, (2) subtraction of a 180°-rotated version of 
the image centered on the star from itself, (3) high-pass 
filtering by subtracting a Gaussian-smoothed (Gaussian 
FWHM = 1-3 x PSF FWHM) version of the image from 
itself, and (4) simple subtraction of an azimuthally medi- 
aned radial profile. We found that (1-3) gave comparable 
results, while (4) did not perform as well as the rest be- 
cause of the four-cornered or six-spoked symmetry of the 
PALAO or Keck AO PSFs. Even though (3) is arguably 
the most widely used method for PSF subtraction when 
separate PSF observations are not available and when 
the observations were not taken using ADI, we found 
that bec ause of the central symmetry of the brig htest AO 
speckles (|Boccaletti et al.ll2002l : lBloemhoj[2003[ ) method 
(2) worked almost as well. Method (2) also did not al- 
ter the photon statistics of the PSF-subtracted image 
in the spatially correlated manner incurred by Gaussian 
smoothing. Therefore, for PSF subtraction we relied on 

2 The PHARO and NIRC2 detec- 
tor linearization routines are available at 

http : //www . astro . caltech . edu/palomar/200inch/palao/Pharo/pharo . html 
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method (2) the most. 



3.3. Precision Astrometry 

Multi-epoch astrometry is essential for demonstrating 
physical association of bound pairs. This is the principal 
method employed for candidate companion confirmation 
here. Below we discuss the calibration steps that we 
undertook to ensure self-consistent astrometric measure- 
ments throughout our campaign. 

We calibrated our astrometry by obtaining repeated 
measurements of the positions of well-known visual bi- 
naries at each observing epoch. We selected binary 
stars wit h well-known ephemeris from the Sixth Orbi t 
Catalog (jHartkopf et alJl200H : IHartkopf fc Masonll2003D . 
combining binaries with grade 1 (accurately determined, 
short-period) and grade 4 (less accurately known, longer- 
period) orbital s olutions, as recomm e nded for astromet- 
ric calibration by IHartkopf fc Ma son ( 2003). Despite the 
lower quality of the orbital solutions for the grade 4 bina- 
ries, their periods are generally much longer, so that their 
motions are predicted with sufficient accuracy for many 
years into the future. The selected calibration binaries 
and their orbital parameters are given in Tabled] 

The above astrometric calibration was adequate for de- 
tecting astrometric signals Sp/p > 1% with PHARO. 
Such accuracy allowed the confirmation of the first 
brown dwarf companion in our survey HD 49197B 
(jMetchev fc Hillenbrandl I2004D . However, that initial 
calibration assumed that the pixel scale and field ori- 
entation over the entire PHARO detector were well- 
determined from measurements taken near the center 
of the array, ignoring possible image distortion in the 
focal plane. In reality, the PHAR O beam is known to 
be distorted (|Havward et al.|[2001[ ). Accurate character- 
ization of this distortion was necessary before consider- 
ing the results from our complete survey, which focused 
on stars with small proper motions (10-100 mas yr _1 ; 
Fig-EJs) and detected candidate companions over the en- 
tire 25.6" x 25.6" PHARO FOV. 

We arrived at an improved astromet ric calibra t ion of 
the PHARO 25 mas pix -1 camera in iMetchevI (|2006l 
§4), where we determined the full extent of the focal 
plane distortion over the entire array and solved for its 
dependence on telescope hour angle, declination, and 
orientation of the CR rotator. For that calibration we 
used a custom-made astrometric mask with pinholes 
distributed on a rectangular grid that we inserted in 
the telescope beam path at the Cassegrain focus. From 
exposures taken with the mask in place we measured 
the variations in the spacing among the pinhole images 
with changes in the instrument gravity vector. We 
found that the PHARO pixel scale varied by up to 
Sp/p = 0.9% from the center to the corner of the 
array in the 25 mas pix -1 camera. After fitting two- 
dimensional polynomials to the distortion, we calibrated 
the variation to within 0.15% over the entire chip. 
The polynomial fits to the focal plane distortion on 
the PHARO 25 mas pix -1 camera and its dependence 
on telescope pointing are given by Equations 4.1-4. 4 
and 4.7-4.11, and Tables 4.4-4.5 in IMetchevI (|2006f ). 
An IDL program that corrects for the distortion at 
an arbitrary coordinate on the PHARO 25 mas pix -1 
camera is available at the PHARO instrument web page 



http : //www . astro . caltech.edu/palomar/200inch/palao/Pharo 

A similar astrometric calibration has already been per- 
formed for all three NIRC2 cameras during the pr e-ship 
testing of the instrument ([Thompson et all 12001). Be- 
cause NIRC2 sits on the Keck II Nasmyth platform and 
thus has a constant gravity vector, the distortion of the 
camera pixel scales does not change with telescope point- 
ing. We implemented the existing astrometric calibration 
of the NIRC2 cameras in the analysis of our Keck AO 
imaging data. 4 

3.4. Photometry 

We used l-2xPSF FWHM-diameter apertures for ob- 
ject photometry, with the smaller apertures used on 
fainter sources for higher signal-to-noise measurements. 
The diffraction-limited FWHM of the X s -band PSF of 
PAL AO was consistently «0"1 while for the Keck AO 
system it was «0"05. The local background was mea- 
sured around each object in an annulus with a wide 
enough inner radius so that the halo of the point source 
did not affect the background measurement. The inner 
radius was as small as 1.5 x PSF FWHM for faint sources 
embedded in the halos of bright stars, or as large as 25- 
30 x PSF FWHM for the target primaries. The variations 
in the sizes of the apertures and of the background an- 
nuli resulted in photometric uncertainties on the order 
of 0.10-0.30 mag. Uncertainties of >0.5 mag were found 
in a few isolated cases involving very faint point sources 
and/or point sources near the edges of the FOV, where 
the PSF was noticeably distorted by anisoplanatism and 
circular apertures did not produce accurate photometry. 

PSF-fitting, rather than aperture photometry was used 
to measure the fluxes of closely-separated point sources. 
The photometric uncertainties in such cases were gener- 
ally < 0.20 mag. 

For absolute calibration we relied on the 2MASS fluxes 
of the primaries. Photometric measurements were always 
obtained relative to the fluxes of the target primaries, as 
measured from the unocculted, short exposures, often 
taken with the PHARO ND 1% filter in place. We cali- 
brated the near-IR extinction of the ND filter from pho- 
tometric measurements of three program stars on images 
taken with and without the filter in place. Images with 
Keck/NIRC2 were obtained only in coronagraphic mode, 
using predominantly the l"-diameter spot, although the 
2"-diameter spot was used during 16-18 May 2003. Un- 
like the PHARO coronagraphic spots, the NIRC2 spots 
are transmissive, offering the possibility to obtain rela- 
tive photometry with respect to the primary. A mea- 
su rement of the throughput of the 2" spot was given 
in iMetchev fc Hillenbrandl (|2004h . Subsequent observa- 
tions showed that such measurements were dependent on 
the quality of the AO correction, possibly because of the 
amount of additive background caused by light from the 
stellar halo diffracted within the area of the coronagraph. 
Thus, approximate relative photometry with the NIRC2 
coronagraph is likely feasible only with good AO cor- 

3 The PALAO/PHARO astrometric calibration was performed 
in March 2005. The optics on the PALAO system have since been 
realigned to accommodate recent and future science instrument 
upgrades. The astrometric calibration presented here is not appli- 
cable to PALAO data taken since 2007. 

4 A more precise astrom etric calibration of the NIRC2 cameras 
has since been obtained by Cameron et al. (2008). 
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rection (usually at H or K bands), when the amount of 
scattered ( "spill-over" ) light within the area of the coron- 
agraph is minimized. Table [5] lists the measured near-IR 
extinction in magnitudes for the PHARO ND 1% filter 
and for the 1" and 2" NIRC2 coronagraphs. The large 
apparent difference in the J-band transmissivity of the 
two NIRC2 coronagraphic spots is a probable effect of 
spill-over (more significant for the smaller spot), aggra- 
vated by poorer AO performance at J. 

4. OBJECT DETECTION AND DETECTION LIMITS 

4.1. Object Detection 

Object detection is a straightforward matter to auto- 
mate in point-source-rich images where the PSF is ra- 
dially symmetric, approximately constant in time, and 
has a well-characterized dependence on image location. 
Unfortunately, none of these qualifications describe the 
sparsely populated high-contrast images in our deep sur- 
vey, in which the main (and frequently only) point source 
is occulted by the coronagraph. In addition, auto- 
mated source finding in AO images of bright stars is 
hindered by large numbers of speckles. Speckles are 
individual images of the star that form from uncor- 
rected and/or induced (by the telescope optics) aber- 
rations in the wavefront, and appear indistinguishable 
from point sources to automated detection routines. 
As a result, even though certain source d etection algo- 
rithm s have been developed (StarFinder; Diolai ti et al. 
20001 ), or adapted (DAOPHO T II, IDAC; |Stetson 
1992t Uefferies fc Christoul [l993). for diffraction -limited 
image restoration, they did not produce satisfac- 
tory results on our images. Our experiments wit h 
DAOPHOT, WAVDETECT (jFreeman et all [2002] ). 
and StarFinder produced large numbers of spurious 
detections, the vast majority of which could be identified 
with speckles around the coronagraph. If the signal-to- 
noise threshold in the source-finding algorithms was ad- 
justed to a correspondingly higher level, the algorithms 
would miss bona fide point sources far from the central 
star. The performance of the automated algorithms did 
not change whether we used various methods of PSF sub- 
traction f£j I3.2[) or not. Similar ex periences and conclu - 
sions were drawn independently bv lCarson et~afl (f2005h . 
who also used the PALAO /PHARO system for their sub- 
stellar companion search. Therefore, after some experi- 
mentation, and despite an understanding that automated 
source detection has the potential to offer greater re- 
peatability and conceptual clarity, we abandoned the ap- 
proach. 

Instead, we opted for visual point source identification, 
which, barring subjective factors, produces superior re- 
sults compared to automated detection. We carefully 
inspected all of the final coronagraphic images for can- 
didate companions. The visual inspection was repeated 
multiple times during the steps of image reduction, pho- 
tometry, and astrometry to reduce the effect of subjective 
factors to a minimum. 

The high-contrast imaging literature abounds 
with examples where the authors have resorted 

to by -eye identificati on of candidate companions 

(e.g., iTokovinin et al l 119991: iBrandner et all 120001: 
Schroeder et al.l 120001: iLuhman &: Javawardhanal 
McCarthy fc Zuckermanl 12004 iMasciadri et al.l 



2002 
200E 



Luhman et al.l 120051 : lLafreniere et al.ll2007f) . A notable 



exception is the study of lLowrance et aTl ((2005), who 
apply a rigorous custom-made automated detection 
scheme to thei r .ffST/NICMOS data. However, the 
lLowrance et al.l survey benefits from the well-behaved 
PSF of space-borne HST im a ging. In a separate in- 
stance, IShatskv fc TokovinTnl |2002) use DAOPHOT ii 
for their non-coronagraphic AO data. Still, they do 
not discuss an application of the approach to their 
set of coronagraphic data, which are likely to be 
speckle-dominated. 

4.2. Determination of Detection Limits 

We quantified our ability to visually detect faint 
objects by introducing artificial point sources in the 
Palomar and Keck Kg images of one of our targets, 
HD 172649, for which data were taken under good ob- 
serving conditions with Stre hl ratios of ss 5 0%. The 
method was first descri bed in IMetchev et aTl (|2003f ) and 
developed more fully in IMetchev! (|2006f) . We summarize 
it here briefly. 

We introduced 1000-5000 artificial point sources of 
constant brightness at random locations over the entire 
25'.'6 x 25'.'6 area of the image and counted the fraction 
of them that were retrievable by eye in 0.25"-1.0"-wide 
concentric annuli centered on the star. We recorded both 
the minimum point source Kg magnitude, at which 100% 
of the artificial point sources were visible at the given 
angular separation, and the maximum Kg magnitude, 
at which only a few artificial point sources were visi- 
ble. We took the mean of the K$ magnitude range as 
the representative limiting magnitude at the given sep- 
aration. We repeated the experiment for a range of ar- 
tificial star magnitudes, at steps of 0.5 mag, on both 
the coronagraphic and the non-coronagraphic images of 
HD 172649. The PSF for artificial stars in the coron- 
agraphic image was obtained from a fit to the brightest 
field object (AK$ = 6.4 mag), whereas in the unocculted 
image the PSF was obtained from a fit to HD 172649 it- 
self. 

The inferred detection limits based on the artificial 
point source experiments are shown in Figure |4] We 
see that the 6 min long Keck AO coronagraphic images 
offered 0.5-1.5 mag higher contrast and up to 0.5 mag 
greater depth than the 24 min PALAO images. The 
greatest difference in contrast is in the l'.'0-l'/5 angu- 
lar separation range, where the presence of waffle-mode 
distortion in the PALAO PSF limits the attainable con- 
trast. 

For the purpose of estimating the completeness of our 
survey, it was important to determine sensitivity limits 
on a per-star basis. However, the above approach was 
too tedious to apply to all observations. Instead, we 
employed a simpler strategy based on the r.m.s. scatter 
of the pixel counts in the radial profile of each sample 
star. To match the approximate spatial correlation scale 
in the PALAO and Keck AO images, we normalized the 
r.m.s. scatter to an aperture with radius equal to the 
O'.'IO FWHM of the PALAO PSF. That is, we multi- 
plied the r.m.s. profile by the square root of the number 
N of pixels in the photometry aperture; N = 50.3 pix 
for PALAO/PHARO with the 25 mas pix _1 camera and 
N = 20.0 pix for Keck AO/NIRC2 with the wide 
(40 mas pix^ 1 ) camera. This procedure imposed a more 
stringent requirement on the significance of the detec- 
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tion of a candidate point source by raising the effective 
multiple of the pixel-to-pixel r.m.s. scatter used as a 
threshold by an additional factor of 4.5-7.1. We show 
the thus-obtained 4cr aperture-normalized r.m.s. noise 
profile of the halo for our PALAO coronagraphic image 
of HD 172649 in Figure H We found that the 4er line 
best approximated the visually determined PALAO de- 
tection limits. The strongest systematic deviation of the 
4(7 r.m.s. profile from the visually-determined contrast 
limits is at angular separations >7". This is to be ex- 
pected, to some extent, because in this region we have 
adjusted the visual detection limits to account for CR 
angle image mis- registration (§ I3.1.1|) . 

The agreement between the detection limits from vi- 
sual inspection and from r.m.s. statistics is dependent on 
a number of factors, such as the radius of the normaliza- 
tion aperture, the PSF pixel sampling, the treatment of 
point source photon statistics (ignored in our r.m.s. anal- 
ysis), and the appropriate functional treatment of non- 
Gaussian sources of error (speckles, shape of the PSF 
core and halo; also ignored here). As a check on whether 
the adopted r.m.s. detection limit approach was a valid 
approximation across the range of PSF and image char- 
acteristics encountered in our survey, we repeated the 
artificial point source experiment on six additional im- 
ages of targets observed both at Palomar and at Keck. 
These images were taken under a range of seeing con- 
ditions, resulting in PSF Strehl ratios between 10% and 
50%. We found that, on average, the by-eye detection 
limits varied between three and five times the level of the 
aperture-normalized r.m.s. noise profile over the entire 
range (0'.'5-12'.'5) of probed angular separations. Thus, 
the additional experiments confirmed that our choice of 
the 4cr level was an adequate detection threshold. 

In closing, we note that because our image noise statis- 
tics in the contrast-limited regime are not Gaussian, the 
adoption of a 4tr threshold does not carry the statisti- 
cal significance of a confidence level at which 99.997% 
of random fluctuat ions are rejected. Only recently have 
iMarois et al.l (|2007h shown that quasi-Gaussian behavior 
of the PSF subtraction residuals can be obtained with the 
ADI technique, allowing such formal estimates on the de- 
tection limits. Because our data were not taken in ADI 
mode, the same formalism can not be applied here. 

4.3. Illustrative Detection Limits for the Deep Sample 

Table [6] lists the attained point-source magnitude sen- 
sitivity for each star in the deep sample at angular sep- 
arations of 1", 2", and 5". Beyond 5" the detection lim- 
its are constant to within 0.5 mag. In the cases where 
multiple images of the same star were taken at differ- 
ent epochs, we have listed the sensitivity only for the 
epoch with the deepest image. We thus formed a set 
of 58 Palomar and 42 Keck images that represented the 
deepest observations of the 100 stars in the deep sample. 

Figured depicts the range of attained iCs-band con- 
trast for the coronagraphic observations in the entire sur- 
vey (thick solid line), and from the Palomar (dotted line) 
and Keck (dashed line) portions of it. Figure [SJd uses 
the same notation to depict the imaging depth of the 
survey in terms of apparent Kg magnitude (i.e., with 
the magnitude of the primary added in each case). The 
median sensitivities of the combined survey range from 
AK S = 8.4 mag at 1" to AK S ~ 12.5 mag over 4"- 



12'.'5 in contrast and from K$ = 15.4 mag at 1" to 
Ks ~ 19.7 mag in depth. These detection limits will 
be used in the Appendix to estimate the completeness of 
the deep survey to substellar and stellar companions. 

We obtained the detection limits for the shallow sam- 
ple in a manner similar to that used for the deep sample: 
from the 4cr dispersion of the radial profile of each star, 
normalized to an aperture with radius equivalent to the 
FWHM of the PALAO PSF. The shallow sample detec- 
tion limits are given in Tabled where we have in addition 
listed the sensitivity at Of! 5. 

In some cases close binary companions elevate the dis- 
persion in the radial profile of the primary, resulting in 
unusually low sensitivities at certain angular separations: 
e.g., for HD 172649 at 5" in Tableland for HD 224873 at 
2" in Table El We have retained these lower sensitivities 
in Tables [S] and [7] as an indication that part of the images 
around the sample stars in question were compromised 
by a nearby bright companion. 

5. CONFIRMATION OF CANDIDATE COMPANIONS 

5.1. Detected Candidate Companions 

In the course of the three year survey we discovered 
287 candidate companions brighter than Ks = 20.6 mag 
within 12.5" of 130 from the 266 sample stars. Of these 
candidate companions 196 were around 61 of the 100 
stars in the deep sample. The remaining 91 were in the 
vicinity of 70 of the 166 shallow-sample targets. All can- 
didate companions around stars in the deep and shal- 
low samples are listed in Tables E] and M respectively. 
Figure [5] shows all detected candidate companions as a 
function of magnitude difference AKs and angular sep- 
aration p. Thirty-nine stars in the deep sample and 96 
in the shallow sample showed no projected companions 
within 12.5". 

5.2. Deciding Physical Association 

The physical association of each candidate companion 
was decided based on one of the following criteria: (1) 
common proper motion with the candidate primary, (2) 
a combination of the position on a J — Ks versus Mk s 
color — absolute magnitude diagram (CAMD; assuming 
equi-distance with the primary) and background star 
density arguments, or (3) extent of the radial profile of 
the candidate companion beyond that of a point-source 
PSF (which suggests an extragalactic object). Candidate 
substellar companions that satisfied the common proper 
motion test were also observed spectroscopically to con- 
firm that their spectral types were in agreement with 
their projected substellar masses. 

5.2.1. Proper Motion 

Proper motion is usually the criterion of choice in 
companion studies, as it provides nearly unambiguous 
evidence of association between two objects: whether 
as components of a gravitationally bound system or as 
members of a multi-star moving group sharing a common 
origin. We used the common proper motion criterion 
through the combined application of two requirements: 
(i) that the change in the position of the candidate com- 
panions relative to the primaries was within 3cr of zero 
in all of right ascension (a) , declination (<5) , angular sep- 
aration (p), and position angle (9) , and (ii) that the 
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expected change in relative positions of the candidate 
companions, had they been stationary background ob- 
jects, was more than 3cr discrepant in either a, S, p, or 
8 from the observed change. Often in cases of candidate 
close (<20 AU) binaries, criterion (i) was not satisfied 
because of appreciable orbital motion. In such situations 
we instead made sure that (iii) the observed change in 
relative position was much smaller (and less significant) 
than the expected change if the components of the can- 
didate binary were not gravitationally bound. A detailed 
example of the impleme ntation of the ab ove astrometric 
criteria is worked out in iMetchevI (|2006l §5.4.1). 

When a relatively bright field star (4 mag < AKs < 
8 mag) was present in the deep coronagraphic exposures 
at Palomar, its position in the shallow non-coronagraphic 
images was used as an additional astrometric reference. 
In cases where the subsequent astrometric measurements 
with respect to the primary and to fainter field objects 
showed such bright field stars to be approximately sta- 
tionary, they could be used to bootstrap the association 
of other candidate companions with the primary, and 
thus circumventing the somewhat higher positional un- 
certainty associated with locating the primary behind the 
opaque PHARO coronagraph. This technique was par- 
ticularly important in determining the association of sys- 
tems in the distant Upper Scorpius (145 pc) and a Persei 
(190 pc) regions, where the primaries have small proper 
motions (< 40 mas yr _1 ) and the images contain multi- 
ple background stars because of the low galactic latitude 
(5° < |b| < 25°). 

5.2.2. Near-IR CAMD and Background Object Density 

Systems with bright (AKs < 5 mag) close-in candi- 
date secondaries often lacked dual-epoch astrometry in 
our survey. Such systems were given lower priority in 
follow-up observations because the companions were con- 
sidered to be stellar and almost certainly bound. Multi- 
epoch astrometric analysis was inapplicable in these 
cases. However, the candidate stellar secondaries in these 
systems were bright enough to be seen in the shallow non- 
coronagraphic JHKs exposures of our targets. Hence, 
for the majority of the candidate stellar systems lack- 
ing astrometric confirmation, physical association could 
be estimated based on the near-IR colors and expected 
absolute magnitudes of the components. 

In evaluating the association of a candidate companion 
based on its near-IR photometry, we placed it on a J—K$ 
versus Mk s CAMD, and checked whether it laid on the 
same isochrone as the primary. In the substellar regime, 
especially near the L/T transition (12 < Mk < 14) 
where the isochrones are not well-constrained, we relied 
on the empirical main sequence as traced by nearby M-T 
dwarfs dLeggett et alj2002HReid et all2004D with k nown 
parallaxes ffrom lDahn et al. ll2002HVrba eHd1l200l . All 
candidate companions with available J-band photometry 
for which the astrometry was inconclusive had their phys- 
ical associations with their candidate primaries evaluated 
in this manner (Fig. [7|. To limit the probability of mis- 
classifying field stars as bona fide companions, positive 
associations were adopted only for candidate companions 
within a 5"-radius field of interest from the primary. 

This approach was successful mostly for stellar-mass 
companions bluer than J — K$ = 0.8 mag, i.e., ear- 
lier than spectral type MO. The main sequence for red- 



der M0-M6 dwarfs is nearly degenerate in J — Ks over 
nearly 4 mag in Mk s (see Fig. [7]) and does not allow 
reliable association estimates from the J — K$ color. At 
even later spectral types, potentially representative of 
young brown dwarf companions, the higher photomet- 
ric uncertainties and the larger empirical color scatter 
at substellar masses prevented the conclusive determina- 
tion of physical association in all but a handful of faint 
projected companions, ii-band photometry, where avail- 
able, did not improve the analysis because of the smaller 
wavelength range sampled by the H — K$ vs. J — Ks 
colors. Thus, no candidate substellar companions were 
confirmed through near-IR photometry. However, a few 
could be rejected. 

In addition to using near-IR colors, it was also pos- 
sible to obtain a probabilistic estimate of the physical 
association for a candidate companion to its correspond- 
ing primary by comparing the number of detected ob- 
jects within the 1275 survey radius to the surface den- 
sity of stars at the relevant galactic coordinates down to 
the limiting magnitude of the survey. Because of the 
lack of large-area deep (Ks < 20 mag) near-IR sur- 
vey data, we limited this type of analysis only to can- 
didate companions in the shallow survey. Although the 
depth of the shallow survey varied depending on the use 
of the ND filter at Palomar, it was roughly compara- 
ble to the 99% completeness limit of the 2MASS cata- 
log: Ks < 14.3 mag in unconfused regions of the sky. 
Therefore, for all candidate companions brighter than 
Ks = 14.3 mag, an empirical estimate of the association 
probability was possible based on 2MASS. Given that the 
faintest primaries in the sample have Ks magnitudes of 
9.6, such a probabilistic analysis could be performed on 
all candidate companions with AKs < 4.7 mag. 

To estimate the contamination from Ks < 14.3 mag 
field stars, we counted the number of 2MASS objects 
within a 5'-radius circular area offset by 12' from each 
sample star (to avoid bright artifacts) , and from that ob- 
tained the expected number of background objects in the 
5"-radius field of interest. We used this as an estimate 
of the purely geometrical chance alignment probability 
(CAP): 

7r(5") 2 

CAP = (number of 2MASS sources within 5')x \ ,,, . 

7r(5 Y 

(1) 

Table fTUl lists the separations, AKs and Ks magnitudes, 
and the CAPs for all sample stars with color compan- 
ions (i.e., the ones with "yes(c)" entries in Tables [8] and 
El). Most color companions have chance alignment prob- 
abilities <2%, with the exception of HD 155902B and 
HE 935B. However, both of these are very close (< 071) 
to their candidate primaries, and are thus almost cer- 
tainly physical companions. (These two systems are in 
fact below the resolution limits of the 95 mas PALAO 
-Ks-band PSF. Their binarity was only appreciated from 
PALAO J-band images, where the PSF is 50 mas wide.) 
The ensemble probability of at least one of the 17 color 
companions being a false positive is 33%, or 16% if 
HD 155902 and HE 935 are excluded. 

5.2.3. Source Extent 

Any bona fide companions to our sample stars were 
expected to be point sources. Apparent source exten- 
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sion could in principle be used to exclude background 
galaxies seen in projection. However, the determination 
of source extent is not a trivial task when the quality 
of the AO correction and hence, the size and shape of 
the PSF change throughout the course of a single night 
depending on guide star brightness and on atmospheric 
stability. In addition, anisoplanatism may radially elon- 
gate point sources PSFs far away from the central AO 
guide star. Therefore, departures from the nominal, 
diffraction-limited PSF size and from a centrally sym- 
metric PSF shape were regarded with caution. These 
were used to classify an object as an extended source 
only when they were in disagreement with the size and 
radial behavior of the profiles of other point sources in 
the same image, if such were present. 

5.2.4. Physical Association Summary 

Using the above criteria, we were able to determine 
the physical association for 198 of the 287 companion 
candidates. The proper motion criterion (§ I5.2.1[) was 
used to establish the majority of associations or non- 
associations: 166 out of 197 (84.3%). These included 55 
bona fide common proper motion companions and 111 
non-common proper motion background objects seen in 
projection. 

The CAMD and chance alignment criteria (§ I5.2.2j) 
are not conclusive in proving physical association. They 
were invoked only when astrometric follow-up was not 
obtained or the proper motion data were ambiguous, 
and when additional J-band images were taken (§ I5.2.2|) . 
These criteria were used in tandem to establish with high 
fidelity the physical association of candidate stellar com- 
panions in the 18 cases listed in Table [lOl The CMD 
criterion alone was used in seven cases to exclude back- 
ground interlopers. 

The source extent criterion (§ 15 . 2 .3[) was used to weed 
out faint galaxies, which may otherwise have red near- 
IR colors, partially due to line-of-sight extinction, and 
may thus pose as candidate substellar objects for the 
preceding criterion. This criterion was applied in four 
cases. 

None of the above criteria were applicable to 89 candi- 
date companions (31.0% of the total) that remained "un- 
decided." The vast majority of these were faint objects in 
the fields of distant (> 100 pc) stars with small apparent 
proper motions (< 50 mas yr _1 ), often at low galactic 
latitudes (b < 15°). These were often discovered only in 
follow-up deeper imaging with Keck and thus lack the 
full time-span of astrometric observations. Judging by 
the large number of such candidate companions per star, 
and based on expectations of the background star con- 
tamination rate at low galactic latitudes, probably none 
of these candidate companions are associated. Through- 
out the rest of the analysis, we shall assume that all 89 
of the undecided candidates are unassociated field stars. 

6. SURVEY RESULTS 

Prelimi nary results from the survey were alre ady pub- 
lished in iMetchev fe Hillenbrand! (|2004 120061 ). includ- 
ing the discovery of two brown dwarf companions, 
HD 49197B and HD 203030B. In this paper we re- 
port the results from the full survey. We found no 
more substellar companions in our sample. We summa- 
rize the findings on the two previously discovered sub- 



stellar companions in § 16.11 We also report 21 new 
stellar companions, in addition to three (HD 129333B, 
H E 373B, and RX J0329.1+0118 B) already announced 
in IMetchev k Hillenbrand! (|2004h . Four of the newly- 
discovered stellar companions, HD 9472B, HE 373B, 
HD 31950B, PZ99 J161329.3-231106B, in addition to 
RX J0329.1+0118B announced previously, have masses 
of only «0.1 Mq, and reside in very low mass ratio 
q = M2/M1 w 0.1 systems. The results on the stellar 
companions are detailed in § 16.21 

A proper motion companion to the star ScoPMS 214, 
considered to be a brown dwa rf based on its apparent K$ 
ma gnitude in IMetchev! ((2006), was found to most prob- 
ably be an unasociated foreground M star after spectro- 
scopic follow-up. This companion, ScoPMS 214"B", is 
discussed in § 16.31 as an example of a pathological case 
where spectroscopic analysis argues against the physical 
association in an apparent common proper motion sys- 
tem. 

Independently of the unbiased survey for substellar 
companions, we also observed and establi shed the phys- 
ical association of a previously known (jBouvier et al.l 
Il997h candidate companion to HII 1348. The estimated 
mass of HII 1348B is near the limit for sustained hy- 
drogen burning. Because of our deliberate inclusion of 
HII 1348 in our observing program based on known bina- 
rity, it is excluded from the present analysis. HII 1348B 
will be the subject of a future publication. 

6.1. Brown Dwarf Companions 

Both brown dwarf companions, HD 49197B and 
HD 203030B, were found in the 100-star deep sur- 
vey. The observed photometric and astrometric prop- 
erties of the two and their inferred masses are listed 
in Table QT] alongside those of the stellar secon- 
daries observed in our survey. The spectral types 
of HD 49197B and HD 203030B are L4 ± 1 and 
L7.5 ± 0.5, respectiv ely, and their masses are estim ated 
at 0.0601^20 M Q (IMetchev fc Hillenbrand! l200l and 
0.0231°;°?* M Q (IMetchev fc Hillenbrand! I200&T) . 

Because of their association with main sequence stars, 
the ages of HD 49197B and HD 203030B are known with 
relative certainty. Their moderate youth (250-500 Myr) 
makes them valuable as benchmarks for substellar prop- 
erties at log g « 5 surface gravities, ^0.5 dex lower than 
the gravities expected of ~3-5 Gyr-old brown dwarfs in 
the field. 

At the time of its discovery, HD 49197B was only the 
fifth known L dwarf younger than 1 Gyr. At a projected 
separation of only 43 AU from its host star, HD 49197B 
was also one of the closest-in resolved su bstellar compan - 
ions, second only to HR 7672B (14 AU: lLiu et aTll2002t ). 
Both HR 7672B and HD 49197B provided early indica- 
tion that the bro wn dwarf desert may not ext end much 
outside of 3 AU (jMetchev fc Hillenbrandll2004f K 

HD 203030B was the first young brown dwarf with a 
spectral type unambiguously as late as the L/T transi- 
tion. Its surprising underluminosity, by w0.5 dex com- 
pared to theoretical predictions for ^1400 K brown 
dwarfs at its age, indicated that its effective tempera- 
ture was w200 K cooler (i.e., «1200 K) than expected 
at the L/T transition. That is, the spectrophotomet- 
ric properties of HD 203030B indicated that either the 
effective temperature at the L/T transition had a hereto- 
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fore unappreciated dependence on surface gravity, or that 
the entire population of field substellar objects had had 
their effective temperatures and ages significantly over- 
estimated. 

In fact, undcrluminosity and <1400 K expected 
effective temperatures are observed in all known 
substellar companions near the L/T-transition 
(|Metchev fc Hillenbrand! 12006ft . including both 
the recently discove red T2.5 dwarf HN PegB 
(|Luhman et al.l 12007 a) and the presumed plan- 
etary mass 2MAS S J1207334-39 3 254B (L5-L9; 
IChauvin et "all l2005al iMohantv et~aT1 l2007f ). With 
the mean ages of their respective primaries ranging 
between 8 Myr and 2 Gyr, all six known L/T- 
trans ition compan ions (GJ 584C. iKirkpatrick et all 



2001 1 G J 337CD, iWilson et alj 120011: iB urgassc r et all 
200a 2MASS J1207334-393254B. IChauvin et alj|2005a 
HP 203030B iMetchev fc Hillenbrandl 120061 HN PegB, 
iLuhman et alj l2007a|) are likely younger than the 
2.9 Gyr model-dependent mean ag e of L/T-transitio n 
dwarfs in the solar neig hborhood (|Allen et all I2005T ). 
Therefore, the theory may indeed be overestimating 
the ages of field brown dwarfs, by a factor of at least 
1.5. This hypothesis has now been independently 
reinforced by the first measurem ent of the dyna mical 
mass of a binary field T dwarf. iLiu et all (|2008l ) find 
that the components of the T5.0+T5.5 dwarf binary 
2MASS J15344984-2952274AB are about 100 K cooler 
than derived for similar field objects: a fact that they 
interpret as evidence for a factor of « 6 ± 3 overestimate 
in the adopted ages of field brown dwarfs. Future 
high-contrast imaging and astrometric observations and 
discoveries of benchmark brown dwarfs with known ages 
and dynamical masses will shed important light on these 
surprising results. 

6.2. Stellar Secondaries 

The entire survey produced 24 new stellar com- 
panions, including the thre e already announced in 
iMetchev fc Hillenbrand! pOOl HD 129333B, HE 373B, 
and RX J0329.1+0118B). HD 129333 had previously 
been identified as a pro bable long-period sin gl e-line d 



1991), 



2005) 



spectroscopic binary by iDuauennov fc Mavorl 
and was independently resolved by iKonig et al.l 
Four other new binaries have since been indepen 
dently co nfirmed in analyses of Hipparcos "problem" 
stars by iMakarov fc Kaplan! (120051. H P 26990 and 
HD 1353631. iGoldin fc Makarovl (120061. HD 152555), 
and IGoldin fc Makarovl (|2007[ HD 155902). In addi- 
tion, PZ99 J161329.3-231106, resolved by us, has since 
been suggested clS cL possible spectroscopic binary by 
iGuenther etU] (|2007f ). 

In addition to the 24 new systems, the physical associ- 
ation of 51 known binary stars was confirmed astromet- 
rically. The star H P 91962, a previously known binary 
(|Mason et al.ll200H and references therein) , was resolved 
into a triple system. No higher-order multiples were re- 
solved. A higher fraction of multiple systems might have 
been expected, especially given the high rate of occur- 
rence (34-96%) of visu al companions to close (spectro- 
scopic) binary systems (jTokovinin et al.ll2006[ ). The vi- 
sual companions in such multiple systems must have ei- 
ther fallen outside of our 12'.'5 survey radius, or been 
removed from our AO sample by the design of the FEPS 



target list, which discriminates against visual compan- 
ions (see § 17. 2|) . 

The majority (57 out of 74) of the binaries plus the 
triple system are members of the shallow survey, as a 
result of the requirement that no AK$ < 4.0 mag candi- 
date companions were present at > 0'.'8 from deep sample 
stars (§ !2.ip . Hence, the binaries found in the deep survey 
either have very low mass ratios, such that the secondary 
is >4 mag fainter than the primary at Kg, or have high 
mass ratios but <0'.'8 angular separations, so that both 
components were fit under the 1" coronagraph. 

We derive ifs-band absolute magnitudes Mk s for the 
companions using the known distances to the primaries 
(Tables [1] and [5]). We estimate the stellar companion 
masses directly from Mk s and from the primary star 
age us ing stellar evolutionary models from iBaraffe et al.l 
(1998). The mass ratios of the resolved stellar binaries 
ranged between 0.1 and 1.0. Including the two substel- 
lar companions, the mass ratios covered the full 0.02- 
1.0 range. Table [IT] lists Mk s and the mass for each 
bona fide companion, along with projected separations 
and system mass ratios. 

6.3. The Apparent Proper Motion Companion to 
ScoPMS 214 

We detected seven projected companions within 12'.' 5 
of ScoPMS 214 (Fig.® Table®). Among these, can- 
didate companion 1 (CC1) is brightest and closest to 
the star, and shares the proper motion of ScoPMS 214 
to within 3tr limits over the course of 4.8 years (specif- 
ically, Aa/o{Aa) = 0.7, A5/a(A5) = 2A,Ap/a(Ap) = 
0.7, A9/a(A6) = 2.3). The apparent proper motion of 
CC1 is significantly different from the remaining three 
candidate companions (2, 3, and 4) to ScoPMS 214 
for which we have sufficiently precise astrometric so- 
lutions. The proper motion of ScoPMS 214 is pre- 
dominantly to the south (p a cosS — —5.6 mas yr" 1 , 
fig = —22.1 mas yr -1 ), and candidate companions 2- 
4 systematically lag behind in their declination motion 
AS, at a level of 4.5-5.3ct(A<5) over 4.8 years. These 
three candidates, along with candidate companion 6, for 
which the astrometry is insufficiently precise to decide its 
proper motion association status, are consistent with be- 
ing stationary distant objects seen in projection (Fig.®). 
Astrometric data for candidate companions 5 and 7 does 
not exist over the entire 4.8-year period, and hence they 
are not plotted on the proper motion diagram in Figure[9l 
However, they are also consistent with being background 
stars. 

In summary, CC1 satisfies all of the proper motion 
association criteria established in § 15.2. 1[ whereas none 
of the other candidate companions to ScoPMS 214 do. 
Therefore, CC1 has a high likelihood of being a bound 
companion to ScoPMS 214, although it could also be an 
unrelated member of Upper Scorpius — the parent asso- 
ciation of ScoPMS 214 — seen in projection. 

It is in principle possible to distinguish between the 
above two possibilities in a probabilistic manner, by fol- 
lowing a two-point cor relation function analysis, as done 
for Upper Scorpius by iKraus fc Hillenbrand! l|200l . We 
find that the probability to find at least one chance align- 
ment among our 23 deep sample targets that belong to 
young stellar associations with an unrelated >M4 dwarf 
(<0.1 Mq) within the same association is ~2.5%. That 
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is, if CC1 is a member of Upper Scoprius, then there 
is a 97.5% probability that it is physically bound to 
ScoPMS 214. Similar reasoning lead us to conclude in 
iMetchevI ff2006h that CC1 was a bona-fide companion to 
ScoPMS 214, which we named ScoP MS 2 14B. 

However, as we shall see in § 16.3.21 the spectro- 
scopic evidence argues against membership of CC1 
(ScoPMS 214 "B") in Upper Scorpius. 

6.3.1. Spectral Type and Effective Temperature of ScoPMS 

214"B" 

We obtained a R « 1200 if -band spectrum of 
ScoPMS 214 "B" (CC1) with Keck AO/NIRC2 on 14 July 
2005. We used the 80 mas- wide slit and the medium 
(20 mas pix -1 ) NIRC2 camera. We integrated for a total 
of 7.5 min on the companion, following an ABC pointing 
sequence with 2.5 min integrations per pointing. We ob- 
served a nearby AO star for telluric correction. Exposures 
of Ne and Ar lamps were obtained for wavelength cali- 
bration. The individual 2.5 min exposures were pair- wise 
subtracted and the spectrum of ScoPMS 214"B" was 
traced and extracted from each exposure in a 280 mas- 
wide (« 5.6 PSF FWHM) aperture. The three individ- 
ual spectra were median-combined and smoothed to the 
resolution set by the instrument configuration using a 
Savitsky-Golay smoothing algorithm. 

The resultant if -band spectrum of ScoPMS 214"B" is 
shown in Figure ITUI where it is compared to IRTF/SpeX 
if-band spectra of field M dwarf a nd M giant standards 
from the IRTF Sp ectral Library 5 (|Cushing et all 120051 : 
iRavner et al.l [20081. smoothed to the same resolution. 
The dominant atomic and molecular absorption features 
due to Na I, Ca I, and CO are identified. 

The overall if -band continuum slope of 
ScoPMS 214 "B" is much closer to the continuum 
slopes of the M dwarfs than to those of the M giants, 
although ScoPMS 214"B" is redder than both sets 
of standards. With an extinction of Ay ~ 2 mag 
towards Upper Scorpius the expected reddening of 
ScoPMS 214 "B" at if band is negligible. Instead, 
the discrepancy between the continuum slopes of 
ScoPMS 214 "B" and the M standards may be due to 
instrumental systematics between the Keck AO/NIRC2 
and IRTF/SpeX spectra. In particular, accurate con- 
tinuum slopes are difficult to extra c t from classical 
AO spectroscopy (|Goto et al.l 120031: iMcElwain et all 
l2007f ) because of the chromatic behavior of the AO 
PSF and because of the narrow slits (here 80 mas) 
used to match the width of the AO PSF. Nevertheless, 
an independent indication that ScoPMS 214"B" has a 
dwarf- like surface gravity (logg ~ 5) comes from the 
relatively shallow depth of the CO bandheads in the 
spectrum of ScoPMS 214"B": comparable in strength 
to the CO bandheads of the M dwarfs and weaker than 
the CO bandheads of the M giants. This is not unusual 
despite the possibility that ScoPMS 214"B" may still be 
contracting toward the main sequence. Even at 5 Myr 
ages M stars are expected to have surface gravities that 
are much more similar to those of dwarfs than to those 
of giants (logg ~ 1). 

From a visual examination of the spectra in Figure [10] 



we estimate that ScoPMS 214 "B" has an M3-M5 spec- 
tral type, based on the relative strengths of the Na I and 
Ca I absorption features compared to the other M dwarfs. 
Unfortunately, a more accurate classification based on 
the if -band spectrum alone is not possible. On one hand, 
the Na 12.21 (im doublet is known to be se nsitive to both 
effect ive temperature and surface gravity (jGorlova et al.l 
120031) . On the other hand, while the Ca I 2.26 /Ltm triplet 
is considere d to be a good temperature indicator for G 
and K stars (|Ali et al.lll995f ) , t he scatter at early- t o mid- 
M spectral types is significant (jGorlova et al.l l2003) . This 
is evidenced by the non-monotonic change in the depth 
of the Ca I triplet in the M3-M6 spectral type sequence 
in Figure [TUl Therefore, we adopt M3-M5 as our final 
estimate of the spectral type of ScoPMS 214 "B." 

The effective temperature corresponding to the M3- 
M5 range is 3250-2800 K (within errors of ±100 K), 
according to the field M dwarf temperature scale of 
iReid fe Hawlevl (|2005l see their Table 4.1). More re- 
cent work on M dwarf effective temperatures, sup- 
ported by highly ac curate photometric and i nterferomet- 
ric measurements (Casagrande et al.l [2008D . finds that 
the IReid fc Hawlevl scale systematically overestimates 
the temperatures of <3000 K fiel d M dwarfs by about 
100 K. However, ILuhmanl ([1999ft finds that young M 
dwarfs specifically are signifi cantly hotter than their older 
field counterparts. iLuhmanl 's conclusion is based on the 
requirement that all components of the GG Tau quadru- 
ple system lie on the same 1 Myr theoretical isoc hrone 
from the NextGen models of iBaraffe et al.l (jUsM), and 
is supported by population age analyses in ot her young 
associations, such as I C 348 (ILuhmanl 11999ft and the 
Orion Ne bular Cluster (jSlesnick et al.ll2004h . Although 
Luhman]'s conclusion relies on theoretical isochrones from 



Baraf feet al.l (|1998l ). the models in question have been 
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shown to most successfully and, on average, fairly accu- 
rately predict t he fundamental properties of pre-main- 
sequence stars (|Hillenbrand fe White! l200l . The effec- 
tiv e temperature range of 1 Myr M3-M5 dwarfs found 
bv ILuhmanl ([l999|) is 3415-3125 K. 

Such disagreement at these low effective temperatures 
is not unusual, given the increasing complexity of stel- 
lar spectra at <3000 K. The problem is even more ag- 
gravated at young ages, when the lower surface gravi- 
ties of the objects further affect their photospheric ap- 
pearance. Because of its specific pertinence to young 
M dwarfs, when considering the possibility below that 
ScoPMS 214"B" is a member o f Upper Scorpiu s, we will 
adopt the temperature scale of lLuhmanl (|1999f ). 

We proceed by examining two probable scenarios: (1) 
a "young" (5 Myr) ScoPMS 214"B" that is a mem- 
ber of Upper Scorpius, probably as a companion to 
ScoPMS 214, with 3125 < T cS < 3415 K, or (2) a "field- 
aged" (1-10 Gyr) ScoPMS 214"B" that is simply seen 
in projection along the line of sight toward ScoPMS 214, 
with 2700 K < T cff < 3250 K. 

6.3.2. Is ScoPMS 214 "B" a Member of Upper Scorpius? 

To decide which of the above two scenarios is valid, 
and by extension, whether ScoPMS 214A and "B" form a 
physical pair, we compare the locations of ScoPMS 214A 
and "B" on an HR di agram with res p ect to the NextGen 
model isochrones of Baraff e~et al.l ([19981 ). Mirroring 
the approach of iLuhmanT l 1999ft . we expect that if 



Adaptive Optics Survey of Young Solar Analogs 



13 



ScoPMS 214A and "B" were bound and hence co-eval, 
they should lie on the same isochrone. Since the temper- 
ature of ScoPMS 214 "B" is ~3000 K regardless of the 
considered scenario, the use of the dust-free NextGen 
models is justified. Indeed, t he more recent DUS TY 
models from the Lyon group (jChabrier et al.l I2000T ) do 
not extend above 3000 K, since dust is not expected to 
form in stellar photospheres at such high effective tem- 
peratures. 

The HR diagram analysis is illustrated in Figure [TTJ 
In the "young" ScoPMS 214 "B" scenario we have 
adopted the mean distance to Upper Scorpius mem- 
bers, 145 ± 40 pc, for both ScoPMS 214A and "B" . 
The bolometric luminosity of ScoPMS 214 "B" is then 
log L/Lq = —2.37 ± 0.24, where we ha ve used bolo- 
metric corre ctions for M3-M5 dwa rfs from lTinnev et ahl 
(|1993h and iLeggett et all (|1996l ). In this scenario, 
ScoPMS 214"B" lies on the 1 Gyr isochrone (i.e., on 
the main sequence), in disagreement with the positioning 
of ScoPMS 214A above the main sequence. Presuming 
that ScoPMS 214A is itself not an unresolved binary, the 
discrepancy indicates that the assumed distance range 
for ScoPMS 214 "B" is incorrect, and that probably it is 
not a member of the Upper Scorpius association. While 
ScoPMS 214A also lies slightly beneath an extrapolation 
of the 5 Myr isochrone, its position is not inconsistent 
with the adopted age for Upper Scorpius, especially given 
the physical extent (~40 pc core radius) of the associa- 
tion. 

In the "field-aged" ScoPMS 214 "B" case, the object 
is not a member of Upper Scorpius, and hence its helio- 
centric distance and bolometric luminosity are not con- 
strained. This case is presented by the shaded region in 
Figure QT] Given the range of luminosities at which the 
shaded region intersects the main sequence, the distance 
to ScoPMS 214 "B" is between 70-145 pc. 

Therefore, we conclude that ScoPMS 214"B" is proba- 
bly not a member of Upper Scorpius and hence prob- 
ably not a physical companion to ScoPMS 214. In- 
stead, it is likely to be a foreground M dwarf seen in 
projection against Upper Scorpius. We arrive at this 
conclusion despite the apparent agreement between the 
proper motions of ScoPMS 214 and ScoPMS 214"B" over 
nearly five years. The reason for the apparent agreement 
is the relatively small proper motion of ScoPMS 214 
(23 mas yr" 1 ), which tests the precision limits of our 
astrometric calibration even over a five-year period. On- 
going astrometric monitoring of this system and mea- 
surements of the individual radial velocities of the two 
components will allow us to discern the difference in their 
space motions. 

7. SURVEY INCOMPLETENESS AND SAMPLE BIASES 

Before addressing the frequency of wide substellar com- 
panions in our sample (§[SJ), we present a brief summary 
of the factors that affect the completeness of our survey 
(§ 17. ip . and discuss the various sample biases (§ 17. 2j) . 
The detailed completeness analysis is relegated to the 
Appendix. 

7.1. Survey Incompleteness 

The principal factors that influenced the completeness 
of our deep survey can be divided into three categories: 
(1) geometrical, defined by the inner and outer working 



angles (IWA and OWA) of the survey (0'/55 and 12'.'5, 
respectively) and by the distribution of heliocentric dis- 
tances of the sample targets; (2) observational, defined 
by the flux limits of the survey relative to the predicted 
brightness of substellar companions; and (3) orbital, de- 
fined by the fraction of orbital phase space observed. For 
the sake of simplicity in estimating the total survey in- 
completeness, we have assumed that the distributions of 
orbital semi-major axis and mass for substellar compan- 
ions are dN/d\oga oc a and dN/dM2 <x M 2 °, respec- 
tively. Other common distributions for these parameters 
are explored in § IA.4L and are found not to affect the final 
completeness estimate by more than a factor of 1.24. 

We find that the combined completeness of the deep 
survey to substellar companions in 28-1590 AU semi- 
major axes ranges from 64.8% at the 0.072 M Q hydrogen- 
burning mass limit to 47.0% at the 0.012 Mq deuterium- 
burning mass limit. The deep survey is severely in- 
complete (<30% completeness) to companions below 
0.012 Mq and maximally complete (64.9%) at and above 
0.090 M Q . Over the combined 0.012-0.072 Mq brown 
dwarf mass range, we estimate that the deep survey is 
complete to 62% of substellar companions with orbita l 
semi-major axes between 28 AU and 1590 AU (§§ IA.3.3T - 
IA.5[) . We combine this estimate with the observational 
results to obtain the underlying substellar companion fre- 
quency in § 18.11 

7.2. Sample Biases 

Our survey sample carries an important bias against 
visual binaries, inherited from the FEPS target selection 
policy. The FEPS sample excluded certain types of visual 
binaries to minimize photometric confusio n in Spitzer's 
F/5-3 0" beam at 3.6-70 /xm wavelengths ()Mever et al.l 
120061 ). In particular: 

1. all FEPS sample stars were required to have no 
projected companions closer than 5" in 2MASS, 
and 

2. stars older than 100 Myr were also required to have 
no projected 2MASS companions closer than 15", 
unless the companions were both bluer in J — K$ 
and fainter at K$ by > 3 mag. 6 

The above two criteria create a non-trivial bias against 
stellar-mass companions in our AO sample. Because of 
the seeing-limited dynamic range of 2 MASS (~ 4.5 ma g 
at 5", - 2.5 mag at 3"; see Fig. 11 in lCutri et aI1l2003h . 
criterion [T] excludes near-equal magnitude (i.e., near- 
equal mass) stellar companions. Criterion [2] then further 
excludes fainter, red (and hence, lower-mass) compan- 
ions, although only around the >100 Myr-old stars. 

Therefore, any analysis of the stellar multiplicity in 
our survey would tend to underestimate the true stel- 
lar companion rate. In particular, if we adopt the me- 
dian distance for the complete AO sample (Table [3]) and 
the orbital period distribution for solar mass binaries 
from lDuquennov &: Mavorl (|1991l ). we find that the above 
FEPS selection criteria have probably excluded ^25% of 
stellar binaries, mostly near-equal mass systems. We do 
not address this bias further. We only note in § 19.11 that 

6 For reasons of generating a statistically significant sample size, 
<100 Myr-old stars were allowed to violate this criterion in FEPS. 
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it has a systematic effect on our estimate of the CMF, 
in the sense that we have underestimated the relative 
frequency of near-equal mass binaries. 

An additional bias against binary stars, relevant only 
to the deep portion of our AO survey, is incurred by 
our on-the-fly selection against AKs < 4 mag projected 
companions at 0'.'8 — 13'.'0 from our deep-sample corona- 
graphic targets (criterion [2] in § 12. ip . However, by keep- 
ing track of which stars were delegated to the shallow 
sample in this manner, we have accounted for this bias 
in our analysis of the CMF in § 19.11 

Finally, our AO sample also carries a slight bias against 
substellar secondaries because of the second FEPS selec- 
tion criterion above. This bias affects only 100-500 Myr- 
old targets in the deep sample with well-separated (>5") 
massive brown dwarf secondaries. Fortunately, because 
of the shallow depth of 2MASS (K s < 15 mag) and its 
limited dynamic range (AKs ^ 6 mag) within our 12'.'5 
AO survey radius, the effect of this bias is negligible. 
Based on the range of assumed semi-major axis distri- 
butions for substellar companions considered in § IA.41 
we find that this criterion would have excluded <0.5% 
of detectable >60 Mj up substellar companions. Over the 
entire substellar companion mass sensitivity range of our 
survey (13-75 Mj up ) the effect of this bias is negligible 
(<0.1%). We will therefore ignore it in the rest of the 
discussion. 

8. THE FREQUENCY OF WIDE SUBSTELLAR 
COMPANIONS 

Throughout the remainder of this paper we will use the 
general terms "substellar companion" and "brown dwarf 
companion" to refer to a 0.012-0.072 M Q (13-75 M Jup ) 
brown dwarf secondary in a 28-1590 AU orbit around a 
young Sun-like star, unless otherwise noted. 

8.1. Results from the Present Survey 

Having discovered two bona fide brown dwarf compan- 
ions among the 100 stars in the deep sample, we esti- 
mate the range of true substellar companion fractions 
that these detections represent. We do so by following 
a Bayesian approach to derive confidence ranges for the 
implied frequency of detectable substellar companions, 
and by applying the incompleteness correction derived 
in§[X5] 

Strictly speaking, the probability of obtaining x suc- 
cessful outcomes (e.g., brown dwarf companion detec- 
tions) from a number of repetitions of an experiment 
with a binary outcome is governed by a binomial dis- 
tribution. In practice, the large number of experiments 
(100) and the small number of successful outcomes in 
our case (a; = 2) mean that the probability of detecting 
x brown dwarfs given an expected mean rate p, is well 
approximated by a Poisson probability distribution: 



P(x\ii) 



H e 



(2) 



We are interested in finding what is the probability distri- 
bution for the actual mean rate li given x detections, i.e., 
we seek the probability density function (p.d.f.) P(/i\x). 

The result follows from B ayes' Theorem 
(|Rainwater & Wdll947l : |Papoulislfl98l : 

P (M|l) = Mri (3) 



where the P's denote "probability distributions" rather 
than identical functional forms. P(li) is the "prior" and 
summarizes our expectation of the state of nature prior 
to the observations. P(x\/i) is the "likelihood" that 
x positive outcomes are observed given a mean of li. 
P(fi\x), the distribution of interest, is the "posterior" 
probability that the state of nature is fi, given x positive 
outcomes. P(x) is a normalization factor equal to the 
sum of all probable outcomes P(x\ii), given the distribu- 
tion of the prior P(ii): 



P(x) 



P(x\ii)P(ji)dn. 



(4) 



We assume no previous knowledge of the state of na- 
ture, and adopt a prior that minimizes the introduction 
of subjective information, imposing only a condition of 
nonnegativity: P(/i) = 1 for /i > 0, P([i) = for /i < 0. 
That is, we assume that all positive substellar compan- 
ion detection rates are equally probable. Inserting Equa- 
tion [2] into Equations [3] and [4] we obtain 



P(li\x) = P(x\li) = 



fi x e M 



(5) 



That is, the p.d.f. of fj, is a Gamma distribution that 
peaks at the observed detection rate x (Fig. [T2|) . Due 
to the asymmetry of the Gamma distribution, the mean 
value (li) is higher than the most likely value llml- (a 4 ) = 
x + 1 = 3 > /jml- 

We determine the confidence interval [///, llh] of the fre- 
quency of substellar companions li at a desired confidence 
level CL by integrating P(li\x) between m and Lih- We 
set the lower and upper bounds lli and ll u of the confi- 
dence interval CL so that (see Fig. [T2|) 



and 



P(n'\x)dLi' = CL 



P(Hi\x) = P(li u \x). 



(6) 



(7) 



Equations ([6]) and |(7j) define the minimum siz e confidence 
interv al \pi,fJ, u ] for confidence level CL (Kraft et al. 
Il991f h The system of equations can not be inverted ana- 
lytically, and has to be solved for \i\ and \i u numerically. 
We do so for the equivalent to the 1,2, and 3 Gaussian 
sigma (68.2%, 95.4%, and 99.7%) confidence intervals. 
The respective confidence ranges for /i are 0.9-3.9, 0.3- 
6.5, and 0.07-9.9 detectable brown dwarf companions for 
a survey of 100 stars. 

Having thus addressed the statistical uncertainties as- 
sociated with the small number of companion detections, 
we now apply the estimated survey completeness correc- 
tion (62%) to /iml and to the confidence interval limits 
of fi. We find that the most likely rate of occurrence of 
brown dwarf companions in 28-1590 AU orbits around 
3-500 Myr-old F5-K5 stars is /j M l = 2%/0.62 = 3.2%. 
The confidence intervals on this estimate are 1.5-6.3% 
at the la level, 0.5-10.5% at the 2a level, and 0.1- 
16.0% at the 3a level, and ar e not a strong function 
of the prior (jKraft et al.lH99l . The mean frequency, 
3%/0.62 = 4.8%, is higher than the most likely value, 
but the exact value of the mean frequency is depen- 
dent on the Bayesian prior. The higher mean frequency 
of wide b rown dw a rf com panions (6.8%) that we re- 
ported in iMetchevI (|2006h was due to the inclusion of 
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ScoPMS 214 "B" as a substellar companion. We have 
now shown that ScoPMS 214"B" is most probably an 
unrelated foreground star seen in projection along the 
line of sight towards ScoPMS 214 (§E|. 

Our results for the frequency of substellar compan- 
ions are built upon simple assumptions for the semi- 
major axis and mass distributions of substellar secon- 
daries (§ 17.11 for greater detail, see § IA.2[) . However, 
our conclusions do not depend strongly on these assump- 
tions. As we show in § IA.4[ when either or both distri- 
butions are varied within empirically reasonable limits, 
the substellar companion frequency remains unchanged 
to within a factor of 1.24. If, as we argue in § 19. 2\ the or- 
bital period distribution of substellar companions is the 
same as for stellar companions, our frequency estimate 
is accurate to within a factor of 1.06. 

8.2. Comparison to Previous Companion Searches 
8.2.1. Radial Velocity Surveys 

Precision radial velocity surveys for extrasolar plan- 
ets have revealed that brown dwarf secondaries are un- 
usually rare (< 0.5%) in 0-3 AU orbits from G and K 
stars: a phenomenon t ermed "the brown dwarf desert" 
(|Marcv fc Butler|[200l . The dearth of brown dwarfs in 
radial velocity surveys is evident with respect to the ob- 
serve d 0-3 AU frequencies of both extra-solar p lanets (5- 
15%; IMarcv fc Butler|[2000HFischer et al.ll2002fl a nd stel- 
lar secondaries (ll%; lDuauennov fc Mavoilll991h around 
Sun-like stars. That is, brown dwarfs are « 20 times 
rarer than planets and stellar companions in 0-3 AU or- 
bits. 

We found that 3. 2+2'^% (2<r confidence interval) of 
young Sun-like stars have 0.012-0.072 Mq companions 
with semi- major axes between 28 and 1590 AU (§ 18.1ft - 
The much wider orbits probed in the present survey pre- 
vent a direct parallel with the radial velocity results. 
Nevertheless, at face value the evidence indicates that 
the frequency of wide brown dwarf companions to Sun- 
like stars is, on average, a factor of ~ 3 smaller than that 
of 0-3 AU extrasolar planets, and a factor of ~ 6 greater 
than the frequency of 0-3 AU brown dwarfs. The dif- 
ference with the exoplanet frequency is not statistically 
significant. The frequencies of 28-1590 AU and 0-3 AU 
brown dwarfs differ at the 98.6% significance level. That 
is, wide brown dwarf companions to Sun-like stars are 
roughly comparable in frequency to radial velocity ex- 
trasolar planets, and are probably more common than 
radial velocity brown dwarfs. 

8.2.2. Wide Stellar Companions 

Based on the iDuquennov fc Mayor! (|199lD orbital pe- 
riod distribution and multiplicity of Sun-like stars, the 
frequency of 28-1590 AU stellar companions is «24%. 
Our estimated frequency of brown dwarfs is a factor of 
~ 8 smaller, and significantly (at the 1 — 10~ 8 level) so. 
Therefore, brown dwarf secondaries are indeed less com- 
mon than stellar secondaries in the 28-1590 AU orbital 
range. 

8.2.3. Other Direct Imaging Surveys for Substellar 
Companions 

A large number of direct imaging surveys have been 
completed to date, covering a wide range in primary mass 



and in sensitivity to substellar companions. Despite the 
disparate characteristics of these surveys, there are now 
enough data to analyze the ensemble of the results. 

We compare our AO survey to all previously pub- 
lished direct imaging surveys for substellar companions 
to >0.2 Mq primaries. We include only surveys target- 
ing >15 stars that also contain at least a cursory ref- 
erence to the parent sample statistics and to the sub- 
stellar companion discovery rates. All such surveys, to 
our present knowledge, are listed in Table [T3J Addi- 
tional direct imaging surveys certainly exist. However, 
any published results from these have tended to report 
only individual detections. To this list of direct imaging 
surveys we have also ad ded the radial- velocity results of 
iMarcv fc Butlerl (2000) for comparison. For each pub- 
lished survey, Table [TH lists the number, median spectral 
type, age, primary mass, and heliocentric distance of the 
sample stars. For most surveys, these values have been 
inferred from the description or listing of the sample in 
the referenced publication. For some surveys, however, 
these parameters have been inferred based on the stated 
focus of the survey (e.g., Sun-like stars), or where ap- 
propriate, based on the properties of stars in the solar 
neighborhood. Table [T^] also lists the maximum probed 
projected separation, the sensitivity to companion mass, 
the number of detected brown dwarf companions, and 
the rate of detection of brown dwarf companions. 

Although an incompleteness analysis is crucial for the 
correct interpretation of survey results, the majority of 
published results do not contain such. Therefore, any 
comparison among surveys has to be based solely on the 
mean or median survey sample statistics and sensitivi- 
ties. Taking the ensemble statistics of all direct imaging 
surveys for substellar companions at face value, without 
accounting for their varying degrees of incompleteness, 
we find that the mean detection rate is 1.0 substellar 
companions per 100 stars. Given the very low number of 
detections per survey (0-2), the results from all imaging 
companion surveys are fully consistent with each other. 

We have plotted the substellar companion detection 
rates of all surveys on a primary mass versus outer 
probed separation diagram in Figure 1131 The outer 
probed separation is defined simply as the product of 
the survey angular radius (generally, the half-width of 
the imaging detector) and the median heliocentric dis- 
tance for the survey sample. The diagram reveals that 
the surveys with the highest detection rates of substel- 
lar companions reside in a distinct locus in the upper 
right quadrant of the diagram, delimited by the dotted 
line. All surveys outside of this region have detection 
rates < 0.6%, whereas all surveys within the region have 
generally higher, 0.5-5.0% detection rates. This fact 
transcends survey sensitivity considerations. Some of 
the most sensitive companion surveys with the smallest 
likely degree s of incompleteness, such as the radial veloc- 
ity survey of IMarcv fc Butlerl (|2000l ) an d the simultane- 
ous di fferen tial imaging (SDI) surveys of lMasciadri et all 
(120051) and iBiller etaLl (|2007f ). lie outside of the dotted 
region and have detection rates well below 1%. That is, 
unless all of these highly sensitive surveys did not de- 
tect brown dwarf companions through some improbable 
happenstance, a significant population of brown dwarf 
companions apparently exists at >150 AU separations 
from >0.7 Mq stars. Brown dwarf companions appear 
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to be less frequent both at smaller orbital separations 
from Sun-like stars, and at wide separations from lower- 
mass stars. The dearth of brown dwarf companions to 
0.2-0.6 M© stars is likely due to a combination of the 
lower multiplicity rate of low mass stars and the ten- 
dency of low mass binaries to exist p redominantly in 
close-in near-equ al mass systems (e.g., iBurgasser et al.l 
l2007t lAllenll2007l and references therein). The surveys 
with the highest detection rates are those targeting very 
wide companions to ~ 1 Mq stars. 

It is important to re-iterate again that none of the de- 
tection rates for any of the surveys in Figure [13] have 
been corrected for systematic or incompleteness effects. 
In particular, there is a strong bias against the detection 
of substellar companions in narrow orbits in all direct 
imaging surveys because of contrast limitations. In ad- 
dition, the position of each survey along the abscissa is 
based on the median outer probed separation, whereas 
most companions are detected at smaller projected sep- 
arations. Therefore, the increase in the frequency of 
substellar companions to >0.7 Mq stars probably be- 
gins well within 150 AU. In the § 19.21 we argue that the 
peak of the brown dwarf companion projected separa- 
tion distribution ma y in fact occur at ^25 AU, a s would 
be expected from the lDuauennov fc Mavorl (|1991h binary 
orbital period distribution. 

9. DISCUSSION 

9.1. The Sub-Stellar and Stellar Companion Mass 
Function 

The salient characteristic of the present imaging sur- 
vey is its high sensitivity to low-mass (M2 < 0.1 Mq) 
companions to solar analogs, i.e., to systems with mass 
ratios q < 0.1. We found only seven such companions 
among 74 binary and one triple systems: the two brown 
dwarfs HD 49197B and HD 203030B, and the 0.08- 
0.14 Mq stars HD 9472B, HE 373B, RX J0329.1+0118B, 
HD 31950B, and PZ99 J161329.3-231106B (Table EJ. 
A naive expectation from the MF of isolated objects 
tiKroupal l200l: iChabrierl l200lh would require approxi- 
mately as many < 0.1 Mq companions as there are 
> 0.1 Mq companions. Therefore, it appears that there 
is a dearth of widely-separated both substellar and low- 
mass stellar companions to Sun-like stars. 

To assess the reality and magnitude of this effect we 
need a uniform survey of a well-characterized sample of 
binaries. Unfortunately, our full survey sample is not ad- 
equate for such an analysis because the imaging depths 
of the deep and the shallow sub-surveys are vastly differ- 
ent, and because the sample is subjected to the combined 
effect of three different biases against binary stars (§ 17. 2j) . 
We could, in principle, focus only on the deep survey of 
100 young stars, for which we have a well-characterized 
completeness estimate. However, doing so would not 
avoid any of the binarity biases. Furthermore, the deep 
survey sample contains only 19 of all 75 binaries and 
triples, only six of which are in the 0'.'55-12'.'5 angular 
separation range, for which we estimated incompleteness 
(§ 17. ip . This number is too low for a meaningful analysis 
of the CMF. 

Nevertheless, we can un-do some of the binarity biases 
and recover certain rejected stellar secondary compan- 
ions by re-visiting how binaries were removed from the 
deep sample during survey observations. We discuss this 



procedure and reconstruct a sample that is minimally 
biased against binaries in the following. 

9.1.1. Defining a Minimally Biased Sample 

We construct a less biased, larger sample of stars by 
adding to our 100-star deep sample all other <500 Myr- 
old stars that were initially selected to be in the deep 
sample but for which no coronagraphic exposures were 
taken because of the discovery of a close-in (0'.'8-13'.'0) 
AKs < 4 mag companion (see § 12. ip . Since we did not 
inherit this bias from the larger FEPS program sample, 
but rather imposed the criteria ourselves, we knew the 
parent sample and were able to un-do the bias exactly. 
The resulting "augmented" deep (AD) sample is mini- 
mally biased against binaries to the extent to which we 
controlled the sample generation. 

There are 28 binaries excluded in this manner, that 
contribute to a total of 128 young stars in the AD sam- 
ple. Among these are a total of 46 binaries and one triple, 
of which 30 systems (including the triple) have compan- 
ions between 0'.'55 and 12'.'5 from the primary. Members 
of the AD sample are distinguished in the last column 
of Table [TTJ where the 30 members with 0'.'55-12'.'5 com- 
panions are marked with "AD30" . 

We assume that the young binaries added from the 
shallow sample do not have additional fainter tertiary 
companions between 0'.'55-12'.'5 that would have been de- 
tectable had we exposed them to the depth of the longer 
coronagraphic images. Given the ~10% ratio of double 
to trip le systems in the study of iDuquennov fc Mayorl 
(|1991h . and the fact that the 28-1590 AU orbital range 
(1=3 10 4 - 7 — 10 7 - 3 days at 1 Mq total mass) includes ap- 
proximat ely 42% of all companions (0- 10 -day periods) 
probed in IDuquennov fc Mavorl (|1991h . we would expect 
that « 0.10 x 0.42 4% of systems have a tertiary 
component in a 28-1590 AU orbit. In comparison, the 
1 : 46 ss 2% ratio of triples to binaries in our sub-sample 
indicates that such an assumption is not unreasonable: 
on average, we may have missed one low-mass (possi- 
bly substellar) tertiary component. Therefore, we have 
potentially suffered only a small loss in completeness by 
including stars for which we do not have deep corona- 
graphic exposures. 

The AD and the AD30 samples remain biased against 
binaries, although mostly against near-equal mass sys- 
tems (§ 17. 2p . Because we have placed an upper age limit 
of 500 Myr for membership in these samples and because 
of the logarithmically uniform distribution of stellar ages 
in the parent FEPS sample, the bias against lower mass 
stellar secondaries (FEPS binarity criterion [2l § 17. 2p af- 
fects less than a quarter of the stars in the AD sample: 
only those that are 100-500 Myr old. 

The detectability of the AD30 secondaries within the 
greater AD sample is subject to the same set of target 
selection criteria and to the same geometrical, observa- 
tional, and orbital incompleteness factors (see § IA.3P as 
for the deep survey. Therefore we can estimate the com- 
pleteness of the AD sample to the detection of secon- 
daries in 28-1590 AU semi-major axes in the same man- 
ner as done for the deep survey. 

The completeness to 0.012-0.072 Mq substellar com- 
panions in 28-1590 AU semi-major axes in the deep sur- 
vey ranges from 47.0% to 64.8%, depending on compan- 
ion mass (see § IA.3p . For masses >0.090 Mq the deep 
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survey is maximally complete at 64.9%. The integrated 
completeness to 0.01-1.0 M© companions is ps 64% (cf., 
62% integrated completeness to 0.012-0.072 M© substel- 
lar companions; § 17. ip . Given the 30 0'.'55-12'.'5 bina- 
ries in the AD30 sample, we would expect a total of 
30/0.64 « 47 ± 9 binaries with 28-1590 AU semi-major 
axes in the 128-star AD sample, where the error on that 
estimate is propagated as -\/30/0.64. (By pure coinci- 
dence, this is exactly how many multiple systems (47) 
are present in the AD sample.) The incompleteness- 
corrected frequency of 0.01-1.0 M© companions in 28- 
1590 AU orbits in the AD sample is thus 47/128 = 
37 ± 7%. This is somewhat higher than the 24% in- 
tegrated ove r the corresponding 1 4 7 -10 73 -day orbital 
period from iDuquennov fc Mavorl (jl991) . Despite the 
bias against binaries, the higher multiplicity fraction of 
stars in our survey is not unexpected because of our su- 
perior sensitivity to very low mass companions and our 
focus on young st ars, which tend to more often be found 
in multiples (e.g-. lGhez et al.lll993l fl997h . 

9.1.2. The Distribution of Companion Masses 

In their G dwa rf multiplicity study, 
IDuquennov fc Mavorl || 19911 ) found that the MFs of 
isolated field stars and of 0.1-1.0 M© stellar compan- 
ions to solar-mass primaries were indistinguishable. 
We now re-visit this conclusion in light of our more 
sensitive imaging data and in the context of more recent 
determinations of the field MF. 

The mass ratio distribution for our selection of 30 
young binaries in the AD30 sample is shown in Figure [T4l 
The distribution is fit well by a power law of the form 

1o z(tT^-) =Slogq + b, (8) 
\d\ogqJ 

equivalent to dN/d\ogq oc q s , or dN/dq cx q 5 ^ 1 = q 13 . 
The best- fit value for the power law index is 8 = 0.61, 
or equivalently, (3 — 0.39 « 0.4. The reduced x 2 of the 
fit is adequate, 1.5, and given only three of degrees of 
freedom a higher-order functional fit is not warranted. 
The x 2 contours of (3 and b indicate that the parameters 
are correlated. By integrating over all possible values 
for b, we find that the 68% (one Gaussian a) and 95% 
confidence intervals for (3 are —0.75 < j3 < —0.03 and 
-0.93 < (3 < 0.14, respectively. 

We compare this mass ratio dist ribution to the k nown 
MF of isolated field objects from IChabrien (|2003f) . Be- 
cause the masses of the primary stars in our sample are 
distributed closely around 1 M© (Fig. []J>), we can di- 
rectly compare the distribution of the (unitless) mass ra- 
tios to the field MF (in units of M@). That is, the mass 
ratio distribution of our sample is essentially equivalent 
to the CMF in units of M© since q = M 2 jM 1 w M 2 /M©. 
The power law index (3 of the CMF is analogous to the 
linear slope a (Salpeter value —2.35) of the field MF. 
As is evident from Figure [14j the CMF of our sample of 
young binaries is very different (reduced x 2 = 7.6) from 
the MF of field objects. 

A potentially more sensitive comparison between the 
observed mass ratio distribution and any model MFs 
(e.g., log- normal, power law) could be obtained using 
a Kolmogorov-Smirnov (K-S) test. We do not perform 
Monte Carlo simulations to degrade the MF models to 



match the observed data, as would be necessary in the 
rigorous sense, but instead compare the models to the 
incompleteness-corrected data. Although the K-S test is 
not strictly applicable with such an approach, the results 
from the test are nevertheless illustrative. Thus, a one- 
sample Kolmogorov-Smirnov test finds only a 2 x 10~ 8 
probability that the o bserved CMF or iginates from the 
log-normal field MF of lChabrierl (f2003h . The K-S prob- 
ability that the fitted power law in Equation [8] with 
(3 = —0.39 is the correct parent CMF is 7%. Ostensi- 
bly the best agreement (58% K-S probability) with the 
incompleteness-corrected data is reached by a log-normal 
mass ratio distribution with a mean and standard devi- 
ation of 0.39. A similar log- normal CMF was inferred 
independently bv lKraus et al.l ([2008) in their analysis of 
resolved binaries in Upper Scorpius. However, we note 
that in our case the difference between the probabilities 
of the power-law and log-normal parent CMFs (7% ver- 
sus 58%) is not statistically significant in the context 
of the K-S test. Therefore, given the already adequate 
reduced x 2 °f the power law fit from Equation [51 we dis- 
regard the potentially better, but statistically less well 
motivated, agreement with the data of the higher-order 
log-normal parameterization (three free parameters), in 
favor of the lower-order (two free parameters) power law. 

A value near zero for our CMF power law exponent 
(3 is consistent with the MF of < 0.1 M r .) objects in 
the fi eld (-1.0 < a < 0.6; iChabrieri |200H iKroupal 
120021: lAllen et al.l HOfll and in young stellar associa- 
tions (-1 < a < 0: iHillenbrand fc Carpenten 120001 
ISlesnick et alll2004l lLhhmanll2004l )." However, the mono- 
tonic rise of the CMF throughout the entire 0.01-1.0 M© 
mass range and the lack of a turnover near 0.1 M© dis- 
agree with MF determinations for stars in the 0.1 -1.0 M Pl 
interv al, where a ranges between —0.5 and —2 ()Kroupal 
|2002| ). That is, in the stellar mass regime, the CMF and 
the MF of isolated stars are distinctly different. 

We should note that the results from our companion 
survey may not be ideally suited for determining the 
CMF of both brown dwarf and stellar companions. In- 
deed, we recall that our AD sample is biased against 
various types of visual binaries IT. 2|) . However, as we 
discussed in § 19.1. 1[ the bias against binarity in the AD 
sample is mostly against near-equal mass systems, the 
secondaries in which would populate the highest mass 
ratio bin in Figure 1141 That is, the power-law index (3 
of the CMF would only further increase in value if the 
bias against near equal-mass binaries in our survey sam- 
ple were taken into account, and the CMF would become 
even more disparate from the field MF. 

Our conclusion counters the established view that the 
MF of 0.1-1.0 M© binary components is indistinguish- 
able from the M F of isolated objects. In a rriving at 
the original result. IDuquennov fe Mavo~rl ()1991l ) had com- 
pared the 0.1 < q < 1.0 binary mass ratio distribution of 
their sample s t ars to an earlier form of the field MF from 
iKroupa et all (|1990f ). Since the mass ratio distribution 
of q > 0.1 binaries in ou r sample is consistent with that 
of IDuquennov fc Mavorl (see Fig. [T4"b). the difference in 
the results stems from our superior sensitivity to q < 0.1 
binaries and from the recently improved knowledge of the 
MF of low-mass (< 0.2 M©) stars in the field. 

Similar conclusions were reache d independently by 
IShatskv fe Tokovm in (2002) and bv iKouwenhoven et al.1 
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(|2005f ) in their direct imaging studies of the vi- 
sual multiplicity of intermediate mass (2-20 Mq) B 
and A stars. These two surveys found that the 
mass ratio distribution of 45-900 AU intermediate 
mass binaries fo llows an f(q) oc q@ p ower law, 
where is -0.5 (IShatskv fc Tokovininl 120021) or -0.33 
(|Kouwenhoven et alj l2005f T Our determination that 
= —0.39 ± 0.36 for companions to solar mass 
st ars indicates that the sha pe of the CMF found 
by IS hat sky fc Tokovininl (|2002f ) and iKouwenhoven et al.l 
(|2005f ) is not specific to intermediate mass stars. Consid- 
ered together, these three sets of results provide a strong 
evidence for a significant difference between the MFs of 
wide secondaries and of isolated field objects. That is, 
the mass ratio distribution of 28-1590 AU binaries is in- 
consistent with random pairing of stars drawn from the 
IMF over a vast range of primary and companion masses. 
We discuss the implications of this conclusion on shaping 
the dearth of brown dwarf secondaries to stars below. 

9.2. The Brown Dwarf Desert as a Result of Binary 
Star Formation 



The inferred 0.01-1.0 M CMF (§ BXg| naturally ex- 
plains the scarcity of wide brown dwarf companions with- 
out the need to invoke formation or evolutionary sce- 
narios specific to substellar companions. The functional 
form of the wide-binary CMF is also consis tent with re- 
sults from radial velocity studies. Thus, Maze h~et al.l 
( 2003) found that the CMF of K-dwarf binaries in 0- 
4 AU orbits is also a rising function of mass over the 0.07- 
0.7 Mq range. Their data are consistent with a power- 
law index of « —0.4, in full agreement with the —0.3 < 
< —0.5 values fo r 28-1 590 AU binaries found by 
Shats kv fc Tokovininl (|2002f) . IKouwenhoven et alj ([20051 



2007D, and here. 



It may be argued perhaps that, given the disparate 
sensitivity systematics and statistical treatments in these 
diverse samples, such an overall agreement is merely co- 
incidental. Indeed, differences in the mass ratio distribu- 
tions of short- vs. long-period binaries within single uni- 
form samples have been previously suggested, with divid- 
ing p eriods of 1000 days (~ 2 A U; IDuquennoy fc Mavorl 
I1990D or 50 days (~ 0.3 AU; Ha lbwachs et all 1 2003). 
Howev er, s ubsequent analyses by IDuquennoy fc Mavorl 
(|1991f) and lMazeh et all (|2003f ) have shown that the evi- 
dence for such discontinuities was inconclusive because of 
relatively small number statistics. At the same time, the 
combined set of direct imaging and spectroscopic data 
referenced here point to an approximately uniform func- 
tional form for the CMF over 1.5 orders of magnitude in 
primary mass (0.6-20 Mq), 3.3 orders of magnitude in 
companion mass (0.01-20 Mq), and 4.7 orders of mag- 
nitude in physical separation (0.03-1590 AU). That is, 
we see strong evidence for a universally uniform shape of 
the CMF. 

Given such universality, it is interesting to consider 
whether the CMF can explain the very low frequency 
of brown dwarfs not only in direct imaging, but also 
in radial velocity surveys. Because stellar and substel- 
lar companions to Sun-like stars appear to be derived 
fro m the same CMF (§19.1.21). we can presume that 
the IDuquennoy fc Mavorl (|1991[ ) period distribution of 
>0.1 Mq stellar secondaries also holds for substellar com- 
panions. Based on this period distribution, the frac- 



tion of all secondary companions in 0-3 AU orbits is 
«22%. Brown dwarfs account for < 0.5%/22% = 2.3% 
of these. For comparison, brown dwarfs account for 
- 3.2%/42% = 7.6% of all secondaries in 28-1590 AU 
periods. In the context of our inferred power-law CMF, 
we find that the value of the index would need to be 
as high as 0.2 to reproduce the ^3 times smaller rel- 
ative frequency of radial velocity brown dwarfs. This 
does not agree well with our 95% confidence limits on 
(-0.93 < < 0.14; §HX1J). However, we also noted that 
our estimate for is systematically underestimated be- 
cause of the bias against near-equal mass binaries in our 
AO sample. It is therefore conceivable that the radial 
velocity brown dwarf desert around G stars represents 
just the low-mass, narrow-orbit end of a CMF that spans 
3.3 dex in secondary mass and 4.7 dex in orbital semi- 
major axis. The problem that would need to be addressed 
then is not why brown dwarf companions specifically are 
so rare, but why the CMF differs so significantly from 
the MF of isolated substellar and stellar objects over all 
orbital ranges. 

In such a universal CMF scenario, brown dwarfs would 
be expected to peak in frequency at semi-major axes de- 
termined by the binary period distribution: at 31 AU 
from solar mass stars, or at projected separations of 
pa 31/1.26 = 25 AU (see ? IA.2I for explanation of factor 
of 1.26). At first glance, this is not consistent with the 
diagram of survey detection rates in Figure 1131 where 
we found that (prior to correction for survey incom- 
pleteness) the highest detection rates occurred in sur- 
veys probing projected separations > 150 AU. However, 
we pointed out that Figure [T3l compares only the median 
outer projected separations probed by the various sur- 
veys, whereas most companions tend to be discovered at 
smaller projected separations f§ 18.2.3)) . In addition, we 
need to consider that projected separations of 25 AU are 
usually well within the contrast-limited regime of existing 
direct imaging surveys of young nearby (50-100 pc) stars. 
Our own survey is less than 40% complete to objects at 
the hydrogen-burning mass limit in 31 AU semi-major 
axis orbits (see § IA.3|> . That is, a number of ~ 30 AU 
brown dwarfs around solar-mass stars may have simply 
been missed in direct imaging surveys because of insuffi- 
cient imaging contrast. 

Unfortunately, neither of the two most sensitive di- 
rect imaging surveys that probe well within 150 AU 
(|Masciadri et al.ll2005l : iBiller et al.ll2007D detect any sub- 
stellar companions. However, their sample sizes are not 
large (54 and 28, respectively), and the null detection 
rates do not place significant constraints on the univer- 
sal CMF hypothesis. Conversely, the recent discovery of 
several pr obable radia l veloc ity brown dwarfs in > 3 AU 
orbits by iPatel et all (|2007l ) lends support to the idea 
that brown dwarfs are more common at wider separa- 
tions, as would be inferred by extrapolation from the 
IDuquennoy fc Mavorl (|199lD orbital period distribution 
for higher-mass, stellar companions. 

Finally, the «0.012 M s («13 M Jup ) deuterium- 
burning mass, above which we limit our analysis, does 
not necessarily mark the bottom of the MF of isolated ob- 
jects. Based on results from three-di mensional smoothed 
part icle hydrodynam i c sim ulations, IBate et all (2002) 
and IBate fc Bonne ll (2005) estimate that the opacity 
limit for gravo-turbulent fragmentation may be as low 
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as 3-10 Mj up . Adopting 3 Mj up as the limit and extrap- 
olating the inferred CMF to <13 Mj up masses, we find 
that sub-deuterium-burning "planetary-mass" compan- 
ions, if able to form through gravo-turbulent fragmenta- 
tion, exist in >30 AU orbits around only <1% of Sun-like 
stars. 

10. CONCLUSION 

We have presented the complete results from a direct 
imaging survey for substellar and stellar companions to 
266 Sun-like stars performed with the Palomar and Keck 
AO systems. We discovered two brown dwarf compan- 
ions in a sub-sample of 100 3-500 Myr-old stars imaged 
in deep cor onagraphic observations. Both were already 
published in lMetchev fc Hillenbrand] (|2004l2006l ). In ad- 
dition, we discovered 24 new stellar companions to the 
stars in the broader sample, five of which are in very low 
mass ratio q ~ 0.1 systems. 

Following a detailed consideration of the completeness 
of our survey, we found that the frequency of 0.012- 
0.072 Mq brown dwarf companions in 28-1590 AU or- 
bits around 3-500 Myr-old Sun-like stars is 3.2±l\% (2a 
limits). This frequency is marginally higher than the 
frequency of 0-3 AU radial velocity brown dwarfs, and 
is significantly lower than the frequency of stellar com- 
panions in 28-1590 AU orbits. The frequency of wide 
substellar companions is consistent with the frequency of 
extrasolar giant planets in 0-3 AU orbits. A comparison 
with other direct imaging surveys shows that substellar 
companions are most commonly detected at > 150 AU 
projected separations from > 0.7 Mq stars. However, 
because of bias against the direct imaging of faint close- 
in companions, brown dwarf secondaries are likely also 
common at smaller projected separations. 

Considering the two detected brown dwarf companions 
as an integral part of the broader spectrum of stellar and 
substellar companions found in our survey, we infer that 
the mass ratio distribution of 28-1590 AU binaries, and 
hence, the MF of 28-1590 AU secondary companions to 
solar-mass primaries, follows a dN/dM 2 (x M| power 
law, with (3 — —0.39 ±0.36 (la limits). This distribution 
differs significantly from the MF of isolated objects in the 
field and in young stellar associations, and is inconsis- 
tent with random pairing of individual stars with masses 
drawn from the IMF. In this context, the observed defi- 
ciency of substellar relative to stellar companions at wide 
separations arises as a natural consequence of the shape 
of the CMF, and does not require explanation through 
formation or evolutionary scenarios specific to the sub- 
stellar or low-mass stellar regime. 



Comparing our CMF analysis to results from other di- 
rect imaging and radial velocity surveys for stellar and 
substellar companions, we find tentative evidence for uni- 
versal behavior of the CMF across the entire 0-1590 AU 
orbital semi-major axis and the entire 0.01-20 Mq com- 
panion mass range. Such a universal CMF is not incon- 
sistent with the marked dearth of brown dwarfs in the 
radial velocity brown dwarf desert around Sun-like stars. 
That is, the properties of brown dwarf companions at 
any orbital separation are conceivably an extension of the 
properties of stellar secondaries. Hence, we predict that 
the peak in semi-major axes of brown dwarf companions 
to solar-mass stars occurs at « 30 AU. Extrapolating the 
inferred CMF to masses below the deuterium burning 
limit, we find that if 0.003-0.012 Mq "planetary-mass" 
secondaries can form through gravo-turbulent fragmen- 
tation, they should exist in > 30 AU orbits only around 
less than 1% of Sun-like stars. 
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APPENDIX 

INCOMPLETENESS OF THE DEEP SURVEY 

Here we examine the factors affecting the sensitivity of the deep survey to substellar companions (§ IA.1[) . and, based 
on several assumptions about the semi-major axis and mass distributions of wide substellar companions (§ IA.2|) . we 
estimate the completeness of the surv ey (§IA.3j) . We find that variations in the parameters of the se mi-m ajor axis and 
mass distributions have little effect f§ IA.4|) on the final completeness estimate. This final estimate (§ IA.5|) is used in §[8] 
in combination with the observational results from our survey to obtain the actual frequency of substellar companions. 

Factors Affecting Survey Completeness 

Several factors need to be taken into account when estimating the detectability of substellar companions to our stars. 
These include: (i) possible sample bias against stars harboring substellar secondaries, (ii) choice of substellar cooling 



20 



Metchev & Hillenbrand 



models, (iii) observational constraints (i.e., survey radius, imaging contrast, and depth), and (iv) physical parameters 
of the stellar/substellar systems (flux ratio, age, heliocentric distance, orbit). 

As discussed in § 17. 2[ the deep sample is largely unbiased toward substellar companions, i.e., factor (i) can be 
ignored. For the bas i s substellar cooling models (ii) we rely on the DUSTY and COND models of the Lyon group 
(jChabrier et al . 2000; Baraffc ct al. 2003). These have been used, either alone or in parallel with the models of the 
Arizona group ([Burrows et al.lll997l L in all other studies of substellar multiplicity. Our choice therefore ensures that 
our results will be comparable with the existing work on the subject. The remaining factors (iii and iv) motivate the 
rest of the discussion here. 

Assumptions 

We will base our incompleteness analysis on three assumptions: (1) that the distribution of semi-major axes a of 
substellar companions to stars is flat per unit logarithmic interval of semi-major axis, dN/dloga oc a (or equivalently, 
dN/da oc a" 1 ) between 10 AU and 2500 AU, (2) that this implies a logarithmically flat distribution in projected 
separations p: dN/dlogp oc p° (i.e., dN/dp oc p _1 ), and (3) that the mass function of substellar companions is flat per 

linear mass interval (dN/dM2 oc Mf = Mj) between O.OIMq and 0.072M©. These assumptions, albeit simplistic, have 
some physical basis into what is presently known about binary systems and brown dwarfs. We outline the justification 
for each of them in the following. 

Assumption (1). — Adopting a total (stellar+substellar) system mass of 1 Mq, the 10-2500 AU range of projected 
separations corresponds approximately to orbital periods of 10 4 — 10 7 5 days. This range straddles the peak (at 
P = 10 4 8 days; o = 31 AU), and f alls along the long-period slope of the Gaussian period distribution of G-dwarf 
binaries (|Ducmennov fc Mavor|[i~99"l . If we were to assume a similar formation scenario for brown dwarfs and stars, 
brown dwarf secondaries would also be expected to fall in frequency beyond ^30 AU separations. However, our limited 
amount of knowledge about brown dwarf companions sug gests the opposite : brown dwarf secondaries may appear as 
common as stellar seco ndaries at >1000 AU separations (|Gizis et al.ll2001[ ), whereas a brown dwarf desert exists at 
<3 AU semi-major axes (jMarcv fc Butlerll2 000: Ma zeh et al.ll2003l ). A smattering of brown dwarfs have been discovered 
in between. A logarithmically flat distribution of semi-major axes for substellar companions, dN/d\oga oc a , or 
equivalently dN/da oc a -1 , represents a middle ground between the known distribution of stellar binary orbits and 
the possible orbital distribution of known brown dwarf companions. The assumption is also attractive because of 
its conceptual and computational simplicity. As we discuss in § IA.4| varying the linear exponent on the semi-major 
axis distribution between and —1, or adopting a log- normal semi-major axis distribution as motivated by the 
iDuauennov fc Mayor! (|1991h binary period distribution, changes the overall completeness estimate by a factor of 

< 1.20. 

Assumption (2). — For a random distribution of orbital inclinations i on the sky, the true and apparent physical 
separations are related by a constant multiplicative factor: the mean value of sini. However, a complication is 
introduced when relating the projected separation to the true semi-major axis because of the need to consider orbital 
eccentricity. Because an object spends a larger fraction of its orbital period near the apocenter than near the pericenter 
of its orbit, the ratio of the semi-major axis to t he apparent separation will tend to values >1. Analytical treatment 
of the problem (|Couteaulll960l : Ivan A lbada 1968) shows that this happens in a n eccentricity-d e pendent manner. Yet, 
when considering the eccentr icity distributions of observed binar y populations (|Kuiped ll935a.b: Duaucnno v fc Mavoil 
H99lHFischer fc Marcvlll992f K both analytical (jvan Albadalll968| ) and empirical Monte Carlo ([Fischer fc Marcvlll992ft 
approaches yield the same identical result: (logo) w (logp) + 0.1. That is, the true semi-major axis and the measured 
projected separation are, on average, related by a multiplicative factor of 1.26, such that (a) = l.26(p). Given 
assumption (1), this then confirms the appropriateness of the current assumption that dN / d\og p oc p°. Furthermore, 
it allows us to relate the projected separations of an ensemble of visual companions to their expected semi-major axes 
in a mean statistical sense. 

Assumption (3). — The assumption for a linearly flat substellar mass distribution (dN/dM2 oc ; [3 — 0) parallels 
results from spectroscopic studies of the in itial mass function (IMF) of low-mass o bjects in star-forming regions 
(jBriceno et al.l l2002l : iLuhman et~aD l2003al fbl: ISlesnick et "all 12004 ILuhman et al.ll2004f ). which are broadly consistent 
with a ~ (wh ere a is the e xpon ent in dN/dM oc M a ). Independently, in a recent analysis of the field substellar mass 
function (MF). [Allen et all (|2005h find a = 0.3±0.6, also consistent with zero. Therefore, assuming that the substellar 
MFs in young stellar associations and in the field are representative of the MF of wide substellar companions, we adopt 
a linearly flat dN/dM2 oc M° CMF for our analysis. This is consistent with our subsequent fit to the CMF in § 19.1.21 
where we determine that (3 is in fact — 0.39± 0.36 over the entire 0.01-1.0 Mq substellar and stellar companion mass 
range. 

The latter result may seem circuitous, since the derivation that (3 is near zero is in fact dependent on the initial 
assumption that (3 is zero. Nevertheless, we find that the initial guess for the CMF exponent is largely unimportant. 
As we discuss in § IA.4| initial values for (3 ranging between — 1 and 1 change the overall completeness estimate by 

< 1.08, and as a result have negligible effect on the final value for (3. 

Incompleteness Analysis 

Adopting the preceding set of assumptions, we now return to the discussion of the remaining factors affecting survey 
incompleteness: factors (iii) and (iv) from § IA.1I We address the individual factors in three incremental steps, as 
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pertinent to: geometrical incompleteness, denned solely by the IWA and OWA of the survey and by the distribution 
of stella heliocentric distances; observational incompleteness, defined by the flux limits of the survey and by the 
predicted brightness of substellar companions; and orbital incompleteness, defined by the fraction of orbital phase 
space observed. These are the same incompleteness categories as already mentioned in § 17.11 Throughout, we adopt 
the aperture-normalized r.m.s. detection limits determined for each star in § 14.31 and assume that the primary ages 
and distances are fixed at their mean values listed in Table [T] 

Geometrical Incompleteness 

In deciding the range of projected separations that the study is most sensitive to, we consider the full range of 
separations that have been explored between the IWA and OWA of the deep survey. For the IWA we adopt 0'.'55, 
i.e., approximately one half width of the O'.'l PALAO i^s-band PSF wider than the 0'.'49 radius of the PHARO 
coronagraph. For the OWA, we adopt 12'.'5, which is 0"3 less than the half-width of the PHARO FOV. Figure Hoi 
shows the fraction of the deep sample stars (solid line) around which successive 1 AU intervals are probed as a function 
of projected separation. It is immediately obvious that only a very narrow range of orbital separations, between 105 AU 
and 125 AU, is probed around 100% of the stars. All other projected separations carry with them some degree of 
incompleteness that needs to be taken into account. From a purely geometrical standpoint, i.e., ignoring imaging 
sensitivity, the limitations imposed by the choice of IWA and OWA amount to a factor of 1.96 incompleteness (for a 
dN/dloga oc a semi- major axis distribution) between 6 AU and 2375 AU: the projected separation range contained 
between the IWA for the nearest star and the OWA for the farthest star in the deep sample. That is, provided that 
substellar companions are detectable regardless of their brightness anywhere between 0'.'55 and 12'.'5 from each star, 
and provided that their distribution of semi-major axes a is logarithmically flat, only about half of the companions 
residing in the 6-2375 AU projected separation range would be detected. 

As is evident from Figure I15[ such a wide range of orbital separations includes regions probed around only a small 
fraction of the stars. Consideration of the full 6-2375 AU range will thus induce a poorly substantiated extrapolation 
of the companion frequency. Instead, we choose to limit the analysis to projected separations explored around at least 
one-third (i.e., 33) of the stars in the deep sample. The corresponding narrower range, 22-1262 AU, is delimited by 
the dashed lines in Figure [T"5l The region has a geometrical incompleteness factor of 1.40 (compared to 1.96 for the 
full 6-2375 AU range above). That is, 1/1.40 = 71.4% of all companions with projected separations between 22 and 
1262 AU should be recovered in our deep survey, if they are sufficiently bright. 

Observational Incompleteness 

Following an approach analogous to the one described in the preceding discussion, we infer the projected separation 
range over which our survey is sensitive to a companion of a given mass. That is, we now take into account that 
not all companions are sufficiently bright to be detected at all probed projected separations. Rather their visibility is 
determined by their expected brightness and the attained imaging contrast. 

Because mass is not an observable, we use the absolute if-band magnitude of a substellar object as a proxy for 
its mass, and employ the Lyon suite of theoretical models to convert between absolute magnitude and mass at the 
assumed stellar age. 

We calculate the observational incompleteness of the deep survey for a grid of 11 discrete companion masses (0.005, 
0.010, 0.012, 0.015, 0.020, 0.030, 0.040, 0.050, 0.060, 072, and 0.090 Mq) and over the entire 3-500 Myr age range 
of our deep sample. We use the DUSTY models from IChabrier et~a l. (2000) when the predic ted companion effect ive 
temperature is above 1400 K (i.e., for spectral types L or earlier), and the COND models from lBaraffe et all (|2003f ) at 
lower effective temperatures (spectral type T) . We compare the estimated companion fluxes at the age of each of our 
sample stars to the corresponding flux limits for each star (see Table |6|), and obtain a minimum projected separation 
at which a companion of a given mass would be visible around each star. Thus, summing over all stars in the deep 
sample, we estimate the observational incompleteness of the entire deep survey to companions of this mass. 

The observational completeness estimates for each of the discrete set of 0.005-0.090 Mq companion masses are shown 
by the filled circles in Figure \Wh . The geometrical completeness limit (i.e., if companion brightness were not a limiting 
factor) is shown by the horizontal continuous line. Figure [T6k demonstrates that the deep survey is nearly as complete 
as is theoretically possible to stellar-mass companions at angular separations between the IWA and OWA, since the 
observational completeness reaches the geometrical limit at 0.090 Mq, just above the minimum hydrogen-burning 
mass. Figure \Wk also illustrates that the observational completeness of the deep survey is >50% for all substellar 
objects above the deuterium-burning limit. The survey completeness drops rapidly below the d euterium-burning lim it 
because of the significantly fainter luminosities expected of non-deuterium fusing objects (e.g.. iBurrows et al.l l2001). 

With the aim to minimize our incompleteness correction, we limit our analysis to substellar companions in the 
0.012-0.072 Mq range, i.e., between the deuterium- and hydrogen-burning mass limits. The sum of the geometrical + 
observational completeness in thi s ma ss range is between 53.0% and 71.3%. Adopting a dN/dM2 oc (i.e., (3 = 0) 
MF for substellar companions (§ I A. 2[) . we find that the observational survey is 68.2% complete to 0.012-0.072 M© 
substellar companions at projected separations of 22-1262 AU from their host stars. 

Orbital Incompleteness 

The analysis so far has dealt only with the projected separation of substellar companions. We now consider the 
effect of realistic orbital semi-major axes, inclinations, and eccentricities. 
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We first adopt the multiplicative factor of 1.26 to relate the projected separation p to the true semi-major axis a: 
(a) = 1.26(/?) (see § IA.2j) . That is, the orbital semi-major axes probed by the survey are on average a factor of 1.26 
further from the star, at 28-1590 AU, than the range of probed projected separations. 

The multiplicative transformation from (p) to (a) does not exhaust the discussion of orbital incompleteness. Because 
companions on orbits with non-zero inclinations and eccentricities spend most of their time at projected separations 
p =/= a 1 1.26, they may still be missed in the survey. The most likely scenarios in which this can occur are for companions 
on highly inclined and/or eccentric orbits. 

With a small number of positive substellar companion detections, orbital incompleteness issues are best addressed 
through Monte Carlo si mulations. Suc h have been performed for a wide r ange of realistic orbital inclinations and ec- 
centricities in a study bv lBrownl (|2004h . the results of which we adopt here. iBrownl 's work investigates the detectability 
of populations of habitable extra-solar terrestrial planets with a range of orbital distributions by the Terrestrial Planet 
Finder- Coronagraph (TPF-C). Although the angular scales and the levels of imaging contrast between the present 
coronagraphic survey and the design specifications for TPF-C are vastly different ( TPF-C projections call for a factor 
of s» 2.5 smaller IWA and ~ 10 6 higher contrast), the problem is conceptually the same: t o determine th e completeness 
to orbits with a certain semi-major axis, given an opaque coronagraph of a fixed radius. iBrownl (|2004f ) parameterizes 
this problem in terms of the ratio Q (which he defines as a) of the semi-major axis to the obscuration radius, so his 
results are universally scalable. H is analysis do es not include treatment of imaging contrast or limiting flux (these 
are addressed in a follow-up work: IBrownl l200lf ). which makes it suit able to apply to results that have already been 
c orrected for th ese effects, as we have already done for our survey in § IA.3.21 

IBrownl (|2004l ) finds that the detectability of orbiting companions in a single-visit observation, what he terms the 
"single visit obscurational completeness" (SVOC), is a s trong functio n of £ between £ = 1 and 2. The SVOC varies 
between «30% at C = 1 and «85% at < = 1.9 (Fig. 3 in IBrownl [200l . Higher SVOC, at the 95% and 99% levels, is 
achieved only for £ = 3.2 and 7.1, respectively, i.e., far from the coronaghraphic edge. The result is largely independent 
(< 10% variation) of the assumed orbital eccentricity e for < e < 0.35. 

We adopt the results o f IBrownl 's analysis and use the SVOC values for a representative orbital eccentricity of 0.35 
(Table 4 in [Br own 200 4]) — a value n ear the peak of the eccentricity distribution of G-dwarf binaries with > 10 3 day 
periods (jDuquennov fc Mavorlll99il ). We calculate the SVOC on the deepest image of each sample star, for each of 
the discrete candidate companion masses in the 0.005-0.090 Mq range considered in § IA.3.21 We define the minimum 
projected separation at which a companion of a given mass is detectable as the effective obsc uration radius for that 
companion mass. The results from the combined treatment of the observational completeness (§ IA.3.2|) and the SVOC 
are shown in FigurefTBk by filled triangles and in FigurefTBb with the dotted lines. The long-dashed lines in FigurclTGk.b 
delimit the maximum attainable SVOC, that is, when the companion brightness is not a limiting factor. Figure ITBk 
shows that the completeness to > 0.072 M Q objects is very near (64.8%) the SVOC limit (64.9%). That is, the deep 
survey is almost maximally complete to stellar companions. The survey is only 47.0% complete to companions at the 
low end of the brown dwarf mass range at 0.012 M . 

The additional consideration of orbital incompleteness does not affect significantly the overall incompleteness of the 
survey within the posited 22-1262 AU projected separation (28-1 590 AU semi-major axis) range. Given the assumed 
companion mass and orbital semi- major axis distributions (§ IA.2|) . the overall (geometrical + observational + orbital) 
completeness becomes 62.0%. 

We note that the consideration of the SVOC, as defined bv lBrownl (I2004I ). does not address all possibilities for orbital 
incompleteness. Other than being obscured by the coronagraph or lost in the glare of its host star, a companion on a 
highly-eccentric orbit may fall outside the OWA, even if its semi-major axis was in the explored range. This additional 
factor, among possible other sources of orbital incompleteness, is not taken into account here. However, judging by the 
small decrease (68.2% — 62.0% = 6.2%) in the overall incompleteness correction induced by the consideration of the 
SVOC, it is unlikely that inclusion of the remaining factors affecting orbital incompleteness will decrease the overall 
survey completeness below 50%. 

Effect of Variations in the Assumed Companion Semi-major Axis and Period Distributions 

The above final completeness estimate is based on the assumptions for the semi-major axis and companion mass 
distributions adopted in § IA.2I These assumptions are merely guesses, and in reality the companion orbital and mass 
distributions may take different forms. Indeed, in § 19.21 we argue that the orbital period distribution of substellar and 
stellar mass companions are probably the same, while in § 19.1. 2l we conclude that the MFs of companions and isolated 
objects are different. Both of these results are at odds with the corresponding assumptions. It is conceivable that 
other initial guesses for the orbital and mass distributions of the companions may lead to different conclusions. 

We therefore analyzed the completeness of the survey to substellar companions under a broader set of functional forms 
for the companion semi-major axis and mass distributions. For the semi-major axis distributio n we also considered: 
(1) th e equivalent of the log-normal orbital period distribution for sun-like binary stars from iDuauennov fc Mayor! 
(|1991[ ) under the assumption that the to tal system mass is 1 M@, (2) the extrasolar planet period dN/d\ogP oc P - 26 
distribution from lCumming et al.l (|2008h . which converts to dN/da oc a~ 61 for solar- mass primaries, and (3) a linearly 
flat dN/da cx a distribution. For the CMF exponent (3 we tested values in the —1 to 1 range. 

The estimates for the completeness to substellar companions for the three different semi-major axis distributions 
with (3 fixe d at zero were 62.7%, 59.9%, and 51.9% respectively, all within a factor of 1.20 of the one already obtained 
in § IA.3.31 In particular, we note that the assumption of either the star-like log- normal or the planet-like period 
distribution altered the completeness estimate very little (by a factor of <1.04). If we set the CMF index (3 to either 
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— 1 or 1 but held the assumed semi-major axis distribution fixed at dN/da oc a" 1 , the completeness became 58.1% 
or 64.9%, respectively. If we allowed both of the companion orbital and mass distributions to vary, the completeness 
estimates ranged from 50.2%-66.0%. 

Overall, we found that the inferred frequency of wide substellar companions to young solar analogs in § 18.11 would 
be affected by a factor of < 1. 24. In the likely case that t he orbital period distributions of substellar and stellar 
companions are the same, as in iDuquennov fc Mayor] (|1991h . our inferred frequency would be accurate to a within 
factor of 1.06. 

Such small changes to the incompleteness estimate of our survey affect the resultant CMF power law index (3 only 
minimally. Because the relative changes in the completeness-corrected numbers per mass bin of the CMF are much 
smaller than the observed trend, the variations in the fitted value for (3 are well within the derived la range. 

Summary of Incompleteness Analysis and Further Considerations 

We adopt 62% as the final estimate for the completeness to substellar companions in our deep survey. That is, 
given two detected brown dwarf companions with semi-major axes in the 28-1590 AU range, on average 0.62 -1 = 1.6 
more companions with semi-major axes in the same range have been missed. This estimate is based on the combined 
consideration of the geometrical, observational, and orbital incompleteness factors described in § IA.3I 

In closing, we recall that because the physical association status of a large fraction (31.4%) of candidate companions 
discovered in the survey remains undecided (§ I5.2|) , it is possible that more bona fide substellar companions may be 
confirmed in this data set in the future. This is not very likely, given that the vast majority of the undecided candidates 
are faint, reside in relatively high-density fields, and are at wide angular separations from their candidate primaries 
(Fig. [5]), i.e., they have very high probabilities of being background stars. Because of the presently unknown and likely 
unimportant nature of the additional candidate companions, and for the sake of preserving statistical rigor, we have 
assumed that none of the remaining candidates are bona fide brown dwarfs, and that the derived value of 62% provides 
an accurate estimate of the completeness of our deep survey. 
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TABLE 1 
Deep Sample 



a 5 n a cos S fij d K$ 



Star 


(J2000.0) 


(J2000.0) 


(mas yr 1 ) 


(mas yr 1 ) 


(PC) 


(mag) 


Sp.T. 


Association 


log Age/ 


HD 377 


00:08:25.74 


+06:37:00.50 


85.2 ± 1.5 


-2.6 ± 1.4 


40.0 ± 2.0 


6.1 


G2V 




7.6 


HD 691 


00:11:22.44 


+30:26:58.52 


209.7+ 1.0 


35.5 ± 1.0 


34.0 + 1.0 


6.2 


K0V 




8.5 


HD 984 


00:14:10.25 


-07:11:56.92 


104.9 ± 1.3 


-67.6 ± 1.2 


46.0 ± 2.0 


6.1 


F7V 




7.6 


HD 1405 


00:18:20.78 


+30:57:23.76 


141.5 ± 2.2 


-177.0 ± 2.1 


29.0 ± 10.0 


6.4 


K2V 




8.0 


QT And 


00:41:17.32 


+34:25:16.77 


44.8 + 0.7 


-36.2 + 0.8 


50.0 ± 25.0 


7.4 


K4 




7.8 


HD 7661 


01:16:24.19 


-12:05:49.33 


134.8 ± 1.1 


-5.7 ± 1.1 


27.0 ± 1.0 


5.7 


K0V 




8.6 


HIP 6276 


01:20:32.27 


-11:28:03.74 


116.0 ± 1.1 


-140.2 ± 1.1 


35.0 ± 1.0 


6.5 


GO 




8.5 


HD 8907 


01:28:34.35 


+42:16:03.70 


51.7 ± 1.0 


-99.2 ± 1.1 


34.0 ± 1.0 


5.4 


F8 




8.3 


HD 12039 


01:57:48.98 


-21:54:05.32 


102.4 ± 1.2 


-48.0 ± 1.1 


42.0 + 2.0 


6.5 


G3/5V 




7.5 


HD 15526 


02:29:35.03 


-12:24:08.56 


42.1 ± 1.3 


-12.2 ± 1.1 


106.0 ± 26.0 


8.0 


G5/6V 




7.6 


1RXS J025216.9+361658 


02:52:17.59 


+36:16:48.14 


53.4 ± 1.3 


-40.1 ± 0.7 


100.0 ± 50.0 


7.6 


K2IV 




7.8 


HD 17925 


02:52:32.14 


-12:46:11.18 


397.3 ± 1.2 


-189.9 ± 1.3 


10.0 ± 0.1 


4.1 


K1V 




7.9 


1RXS J025751. 8+115759 


02:57:51.68 


+11:58:05.83 


31.4 ± 1.2 


-28.4 ± 1.2 


118.0+ 16.0 


8.5 


G7V 




7.8 


RX J0258.4+2947 


02:58:28.77 


+29:47:53.80 


17.4 ± 1.2 


-40.0 ± 0.6 


100.0 ± 50.0 


9.1 


K0IV 




8.0 


1RXS J030759.1+302032 


03:07:59.20 


+30:20:26.05 


31.2 ± 0.6 


-66.6 ± 0.7 


75.0 + 37.5 


7.4 


G5IV 




8.3 


HD 19668 


03:09:42.28 


-09:34:46.46 


88.0+ 1.2 


-113.3 ± 1.1 


40.0 ± 2.0 


6.7 


G8/K0V 




8.4 


IE 0307.4+1424 


03:10:12.55 


+14:36:02.90 


-4.0 ± 1.2 


-25.3 ± 1.2 


160.0 ± 80.0 


8.8 


G6V 




7.8 


V525 Per 


03:19:02.76 


+48:10:59.61 


16.4 ± 4.0 


-23.6 ± 1.4 


190.0 ± 11.0 


9.4 


K2 


aPer 


7.9 


1RXS J031907.4+393418 


03:19:07.61 


+39:34:10.50 


27.3 + 0.9 


-25.3 ± 1.7 


100.0 + 50.0 


9.5 


K0V 




7.8 


HE 622 


03:24:49.71 


+48:52:18.33 


22.3 ± 0.9 


-26.3 ± 0.7 


190.0 ± 11.0 


9.6 


G7 


aPer 


7.9 


IE 0324.1-2012 


03:26:22.05 


-20:01:48.81 


25.0 ± 1.6 


7.4 ± 1.6 


160.0 ± 80.0 


8.9 


G4V 




7.8 


RX J0329.1+0118 


03:29:08.06 


+01:18:05.66 


4.4 ± 1.3 


-4.5 ± 1.3 


100.0 ± 50.0 


9.2 


G0(IV) 




7.8 


HE 1101 


03:35:08.75 


+49:44:39.59 


20.9 ± 1.3 


-28.5 ± 0.9 


190.0 ± 11.0 


9.3 


G5 


aPer 


7.9 


HD 22179 


03:35:29.91 


+31:13:37.45 


42.6 ± 0.6 


-46.0 ± 0.7 


140.0 ± 70.0 


7.4 


G5IV 




7.8 


HD 23208 


03:42:39.80 


-20:32:43.80 


3.8 ± 1.5 


24.1 ± 1.0 


57.5 ± 4.7 


7.2 


G8V 




6.7 


HII 120 


03:43:31.95 


+23:40:26.61 


18.0 ± 0.7 


-46.8 ± 0.6 


133.0 ± 6.0 


9.1 


G5 


Pleiades 


8.1 


HII 2147 


03:49:06.11 


+23:46:52.49 


15.9 ± 0.9 


-43.8 ± 0.8 


133.0 ± 6.0 


8.6 


G7IV 


Pleiades 


8.1 


1RXS J035028. 0+163121 


03:50:28.40 


+16:31:15.19 


26.2 ± 1.3 


-23.4 ± 2.1 


138.0 ± 21.0 


8.6 


G5IV 




7.8 


RX J0354.4+0535 


03:54:21.31 


+05:35:40.77 


-1.4 ± 1.3 


-7.6 ± 1.3 


100.0 ± 50.0 


8.7 


G2(V) 




8.3 


Pels 191 


03:54:25.23 


+24:21:36.38 


17.1 ± 0.7 


-46.8 ± 0.8 


133.0 ± 6.0 


9.1 


G5IV 


Pleiades 


8.1 


RX J0357.3+1258 


03:57:21.39 


+12:58:16.83 


22.7 ± 1.8 


-21.9 ± 1.5 


149.0 ± 23.0 


9.0 


GO 




7.8 


HD 285751 


04:23:41.33 


+15:37:54.87 


8.2 ± 1.7 


-15.8 ± 1.4 


150.0 ± 75.0 


8.8 


K2(V) 




6.8 


RX J0442. 5+0906 


04:42:32.09 


+09:06:00.86 


28.9 ± 2.4 


-22.3 ± 2.0 


119.0 ± 21.0 


9.1 


G5(V) 




7.8 


HD 286179 


04:57:00.65 


+15:17:53.09 


-1.8 ± 1.5 


-17.3 ± 1.4 


140.0 ± 70.0 


8.5 


G3(V) 




7.3 


HD 31950 


05:00:24.31 


+15:05:25.28 


0.3+ 1.1 


-15.2 ± 1.1 


100.0 + 50.0 


8.4 


F8 




7.8 


HD 35850 


05:27:04.77 


-11:54:03.38 


17.5 ± 0.7 


-49.8 ± 0.8 


27.0 ± 1.0 


4.9 


F7/8V 




7.5 


1RXS J053650.0+133756 


05:36:50.06 


+13:37:56.22 


4.9 ± 1.3 


-108.8 ± 1.2 


56.0 ± 28.0 


8.1 


K0V 




8.3 


HD 245567 


05:37:18.44 


+13:34:52.52 


7.5 ± 0.9 


-33.2 ± 0.9 


119.0 ± 21.0 


7.6 


GOV 




6.6 


SAO 150676 


05:40:20.74 


-19:40:10.85 


19.2 ± 1.2 


-12.9 ± 1.2 


78.0 ± 30.0 


7.5 


G2V 




7.8 


HD 38949 


05:48:20.06 


-24:27:50.04 


-29.8+ 1.1 


-37.8 + 1.2 


43.0 + 2.0 


6.4 


G1V 




8.4 


HD 43989 


06:19:08.05 


-03:26:20.39 


10.6 ± 0.9 


-43.7+ 1.0 


50.0 ± 2.0 


6.6 


GOV 




7.8 


HD 49197 


06:49:21.34 


+43:45:32.87 


-37.6 ± 0.6 


-50.9 + 0.6 


45.0 + 2.0 


6.1 


F5 




8.7 


RE J0723+20 


07:23:43.58 


+20:24:58.64 


-66.2 ± 1.8 


-230.2 ± 2.6 


24.0 ± 12.0 


6.9 


K3(V) 




8.1 


HD 60737 


07:38:16.44 


+47:44:55.34 


-14.2 ± 1.0 


-165.0 ± 1.0 


38.0 ± 2.0 


6.3 


GO 




8.2 


HD 70573 


08:22:49.95 


+01:51:33.58 


-49.1 ± 1.1 


-49.7 ± 1.1 


46.0 ± 23.0 


7.2 


G1/2V 




8.0 


HD 70516 


08:24:15.66 


+44:56:58.92 


-63.1 + 0.9 


-178.4 + 1.0 


37.0 + 3.0 


6.1 


GO 




7.9 


HD 72905 


08:39:11.62 


+65:01:15.14 


-28.9 ± 1.0 


88.5 ± 1.0 


14.0 ± 0.1 


4.2 


G1.5VB 




8.3 


HD 75393 


08:49:15.35 


-15:33:53.12 


35.8 ± 1.4 


-33.6 ± 1.2 


42.0 ± 1.0 


5.9 


F7V 




8.4 


HD 82558 


09:32:25.72 


-11:11:05.00 


-248.3 ± 1.2 


35.1 + 0.6 


18.3 + 0.3 


5.4 


K3V 




8.0 


HD 82443 


09:32:43.92 


+26:59:20.76 


-147.5 ± 0.9 


-246.3 ± 0.5 


17.7 + 0.3 


5.1 


K0V 




8.0 


SAO 178272 


09:59:08.42 


-22:39:34.57 


-62.8+ 1.4 


-15.6 ± 1.7 


58.0 ± 29.0 


7.4 


K2V 




8.0 


HD 90905 


10:29:42.23 


+01:29:27.82 


-150.4 ± 0.8 


-124.1 ± 0.8 


32.0 ± 1.0 


5.5 


G1V 




8.3 


HD 91782 


10:36:47.84 


+47:43:12.42 


-71.4 ± 0.6 


-81.7 ± 0.7 


56.0 ± 3.0 


6.8 


GO 




8.2 


HD 92855 


10:44:00.62 


+46:12:23.86 


-268.8+ 1.1 


-61.9 + 1.2 


36.0 + 1.0 


5.9 


F9V 




8.2 


HD 93528 


10:47:31.20 


-22:20:52.80 


-122.7+ 1.1 


-29.4 + 0.8 


34.9 ± 1.2 


6.5 


K0V 




8.0 


HD 95188 


10:59:48.28 


+25:17:23.65 


-126.3 ± 1.4 


1.7 + 1.3 


36.0 + 1.0 


6.6 


G8V 




8.4 


HD 101472 


11:40:36.59 


-08:24:20.32 


-20.0 ± 0.8 


-13.8 ± 0.8 


39.0 ± 2.0 


6.1 


F7V 




8.4 


BPM 87617 


11:47:45.73 


+12:54:03.31 


-71.5 ± 1.9 


-0.4 ± 1.8 


50.0 ± 25.0 


7.8 


K5Ve 




8.1 


HD 104576 


12:02:39.46 


-10:42:49.16 


32.7 ± 1.0 


-18.4 + 0.9 


49.0 + 3.0 


6.7 


G3V 




8.2 


HD 104860 


12:04:33.71 


+66:20:11.58 


-56.1 ± 1.4 


49.7+ 1.4 


48.0 ± 2.0 


6.5 


F8 




7.6 


HD 107146 


12:19:06.49 


+16:32:53.91 


-175.6 ± 0.9 


-149.5 ± 1.0 


29.0 ± 1.0 


5.5 


G2V 




8.0 


SAO 15880 


12:43:33.36 


+60:00:53.28 


-125.2 ± 1.4 


-66.4+1.5 


60.0 ± 20.0 


7.3 


K0 




8.0 


SAO 2085 


12:44:02.88 


+85:26:56.40 


-129.6 ± 0.8 


43.2 ± 0.9 


66.0 ± 20.0 


7.3 


G5 




8.2 


HD 111456 


12:48:39.46 


+60:19:11.40 


107.8 ± 3.1 


-30.6 ± 2.7 


24.2 ± 1.9 


4.6 


F5V 




8.5 


HD 132173 


14:58:30.51 


-28:42:34.15 


-99.9 ± 1.5 


-93.0 ± 1.7 


49.0 ± 2.0 


6.2 


GOV 




8.2 


HD 139813 


15:29:23.61 


+80:27:01.08 


-218.0 ± 1.2 


105.8 ± 1.2 


22.0 ± 0.3 


5.5 


G5 




8.3 


HD 139498 


15:39:24.40 


-27:10:21.87 


-21.8+ 1.5 


-28.1 ± 1.5 


127.0 + 10.0 


7.5 


G8(V) 


ScoCen 


7.2 


HD 142361 


15:54:59.86 


-23:47:18.26 


-29.3+ 1.1 


-38.8 ± 1.1 


101.0+14.0 


7.0 


G3V 


USco 


6.7 


HD 143006 


15:58:36.92 


-22:57:15.35 


-10.6+ 1.7 


-19.5 ± 1.3 


145.0 ± 40.0 


7.1 


G6/8 


USco 


6.7 


PZ99 J155847.8-175800 


15:58:47.73 


-17:57:59.58 


-14.8 + 3.5 


-18.4 + 2.8 


145.0 ± 40.0 


8.3 


K3 


USco 


6.7 


ScoPMS 21 


16:01:25.63 


-22:40:40.38 


-9.4 + 2.8 


-23.8 ± 1.7 


145.0 ± 40.0 


8.5 


K1IV 


USco 


6.7 


PZ99 J160158.2-200811 


16:01:58.22 


-20:08:12.0 


-6.8 + 2.1 


-21.7 + 2.3 


145.0 ± 40.0 


7.7 


G5 


USco 


6.7 


PZ99 J160302. 7-180605 


16:03:02.69 


-18:06:05.06 


-11.3 + 2.9 


-22.7 + 1.7 


145.0 ± 40.0 


8.7 


K4 


USco 


6.7 


ScoPMS 27 


16:04:47.76 


-19:30:23.12 


-14.0 + 2.3 


-20.1 + 3.1 


145.0 ± 40.0 


8.0 


K2IV 


USco 


6.7 


ScoPMS 52 


16:12:40.51 


-18:59:28.31 


-8.4 + 2.4 


-28.5 + 4.1 


145.0 ± 40.0 


7.5 


K0IV 


USco 


6.7 


PZ99 J161318. 6-221248 


16:13:18.59 


-22:12:48.96 


-9.1 ± 1.2 


-21.0 ± 1.4 


145.0 ± 40.0 


7.4 


G9 


USco 


6.7 


PZ99 J161402. 1-230101 


16:14:02.12 


-23:01:02.18 


-8.8+ 1.7 


-22.8 + 1.7 


145.0 ± 40.0 


8.6 


G4 


USco 


6.7 


PZ99 J161411.0-230536 


16:14:11.08 


-23:05:36.26 


-12.1 ± 1.6 


-23.8 ± 1.9 


145.0 ± 40.0 


7.5 


K0 


USco 


6.7 
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TABLE 2 
Shallow Sample 



a 5 n a cos S fij d K$ 1 



Star 


(J2000.0) 


(J2000.0) 


(mas yr 1 ) 


(mas yr 1 ) 


(PC) 


(mag) 


Sp.T. 


Association 


log Age/yr 


HD 224873 


00:01:23.66 


+39:36:38.12 


-28.7 + 0.6 


-43.3 ± 0.7 


49.0 ± 5.0 


6.7 


K0 




8.5 


HD 6963 


01:10:41.91 


+42:55:54.50 


-154.6 ± 0.9 


-198.5 + 0.9 


27.0 + 1.0 


5.9 


G7V 




9.0 


HD 8467 


01:24:28.00 


+39:03:43.55 


210.6 ± 1.8 


-26.6 + 1.0 


31.0 + 1.0 


6.6 


G5 




9.3 


HD 8941 


01:28:24.36 


+17:04:45.20 


118.3 ± 0.7 


-34.8 ± 0.7 


50.0 ± 2.0 


5.4 


F8IV-V 




9.2 


HD 9472 


01:33:19.03 


+23:58:32.19 


0.0 ± 1.0 


28.4 + 0.9 


33.0 ± 1.0 


6.0 


GO 




8.9 


RE J0137+18A 


01:37:39.41 


+18:35:33.16 


65.8 ± 1.9 


-46.0 ± 2.5 


64.0 ± 8.0 


6.7 


K3Vc 




6.8 


HD 11850 


01:56:47.27 


+23:03:04.09 


-83.8 ± 1.0 


-18.1 ± 1.0 


33.0 ± 1.0 


6.2 


G5 




8.8 


HD 13382 


02:11:23.15 


+21:22:38.39 


273.1 ± 0.8 


-12.6 ± 0.7 


33.0 ± 1.0 


5.8 


G5V 




8.7 


HD 13507 


02:12:55.00 


+40:40:06.00 


56.9 ± 1.3 


-99.2 ± 1.3 


26.0 ± 1.0 


5.6 


G5V 




8.9 


HD 13531 


02:13:13.35 


+40:30:27.34 


57.6 ± 1.0 


-96.4 ± 1.0 


26.0 ± 1.0 


5.7 


G7V 




8.7 


HD 13974 


02:17:03.23 


+34:13:27.32 


1153.8 ± 0.8 


-245.1 ± 0.8 


11.0 ± 0.1 


3.2 


GOV 




9.2 


1RXS J025223. 5+372914 


02:52:24.73 


+37:28:51.83 


22.5 ± 0.7 


-24.5 ± 1.0 


170.0 ± 85.0 


9.1 


G5IV 




8.3 


2RE J0255+474 


02:55:43.60 


+47:46:47.58 


79.8 ± 0.6 


-76.1 ± 0.7 


50.0 ± 25.0 


7.2 


K5Vc 




7.9 


HD 18940 


03:03:28.65 


+23:03:41.19 


111.4 ± 0.8 


-0.7 + 0.7 


34.0 ± 1.0 


5.5 


GO 




8.9 


HD 19019 


03:03:50.82 


+06:07:59.82 


231.8 ± 1.8 


50.7+1.7 


31.0 + 1.0 


5.6 


F8 




9.2 


HD 19632 


03:08:52.45 


-24:53:15.55 


226.7+ 1.3 


136.3 ± 1.3 


30.0 ± 1.0 


5.7 


G3/5V 




8.6 


vB 1 


03:17:26.39 


+07:39:20.90 


167.2 ± 1.3 


-6.4 ± 1.4 


43.1 ± 0.6 


6.0 


F8 


Hyades 


8.8 


HE 350 


03:17:36.93 


+48:50:08.50 


23.2 + 0.8 


-23.0 + 0.9 


190.0 ± 11.0 


9.3 


G2 


oPer 


7.9 


HE 373 


03:18:27.39 


+47:21:15.42 


29.0 + 0.7 


-26.8 ± 2.0 


190.0 ± 11.0 


9.4 


G8 


oPer 


7.9 


HE 389 


03:18:50.31 


+49:43:52.19 


22.5 + 0.9 


-23.9 + 0.7 


190.0 ± 11.0 


9.5 


GO 


oPer 


7.9 


HE 696 


03:26:19.36 


+49:13:32.54 


19.8 ± 0.7 


-25.0 ± 0.7 


190.0 ± 11.0 


9.7 


G3 


oPer 


7.9 


HE 699 


03:26:22.22 


+49:25:37.52 


22.4 ± 0.8 


-24.5 + 0.7 


190.0 ± 11.0 


9.4 


G3 


aPcr 


7.9 


HE 750 


03:27:37.79 


+48:59:28.78 


22.0 + 0.7 


-25.6 + 0.7 


190.0 + 11.0 


9.1 


F5 


oPer 


7.9 


HE 767 


03:27:55.02 


+49:45:37.16 


21.1 ± 0.6 


-26.0 + 0.6 


190.0 + 11.0 


9.2 


F6 


oPer 


7.9 


HE 848 


03:29:26.24 


+48:12:11.74 


22.2 ± 0.6 


-26.4 + 0.6 


190.0 ± 11.0 


8.5 


F9V 


aPer 


7.9 


HE 935 


03:31:28.99 


+48:59:28.37 


21.3 ± 0.9 


-26.6 ± 0.6 


190.0 ± 11.0 


8.5 


F9.5V 


oPer 


7.9 


HE 1234 


03:39:02.91 


+51:36:37.11 


21.4 ± 0.8 


-33.7 + 0.7 


190.0 ± 11.0 


8.9 


G4 


oPer 


7.9 


HD 22879 


03:40:22.08 


-03:13:00.86 


691.6 ± 1.1 


-212.8 ± 1.1 


24.0 ± 1.0 


5.2 


F7/8V 




9.3 


HII 102 


03:43:24.54 


+23:13:33.30 


17.1 ± 0.6 


-43.7 + 0.6 


133.0 + 6.0 


8.7 


G6 


Pleiades 


8.1 


HII 152 


03:43:37.73 


+23:32:09.59 


19.5 ± 0.7 


-46.9 ± 1.0 


133.0 + 6.0 


9.1 


G4 


Pleiades 


8.1 


1111 174 


03:43:48.33 


+25:00:15.83 


18.8 ± 1.1 


-47.0 ± 0.9 


133.0 + 6.0 


9.4 


Kl 


Pleiades 


8.1 


HII 173 


03:43:48.41 


+25:11:24.19 


20.4 ± 0.8 


-48.4 ± 0.7 


133.0 ± 6.0 


8.8 


K0 


Pleiades 


8.1 


HII 250 


03:44:04.24 


+24:59:23.40 


20.1 ± 1.0 


-49.4 ± 0.7 


133.0 ± 6.0 


9.1 


G3 


Pleiades 


8.1 


HII 314 


03:44:20.09 


+24:47:46.16 


18.2 ± 0.7 


-49.8 + 0.8 


133.0 + 6.0 


8.9 


G3 


Pleiades 


8.1 


1RXS J034423.3+281224 


03:44:24.25 


+28:12:23.07 


46.4 + 0.7 


-50.6 + 0.6 


49.0 ± 10.0 


7.2 


G7V 




7.8 


HII 514 


03:45:04.01 


+25:15:28.23 


17.3 ± 0.7 


-46.3 ± 0.6 


133.0 ± 6.0 


9.0 


G4 


Pleiades 


8.1 


HII 571 


03:45:15.35 


+25:17:22.11 


15.1 + 0.9 


-48.5 + 0.9 


133.0 + 6.0 


9.2 


G9 


Pleiades 


8.1 


HII 1015 


03:46:27.35 


+25:08:07.97 


18.6 ± 0.7 


-48.5 + 0.9 


133.0 + 6.0 


9.0 


Gl 


Pleiades 


8.1 


HII 1101 


03:46:38.78 


+24:57:34.61 


18.4 + 0.8 


-48.1 + 0.7 


133.0 + 6.0 


8.8 


G4 


Pleiades 


8.1 


HII 1182 


03:46:47.06 


+22:54:52.48 


18.4 + 0.6 


-45.6 + 0.7 


133.0 + 6.0 


8.9 


Gl 


Pleiades 


8.1 


HII 1200 


03:46:50.54 


+23:14:21.06 


17.3 ± 0.6 


-40.2 ± 0.7 


133.0 ± 6.0 


8.5 


F6 


Pleiades 


8.1 


HII 1776 


03:48:17.70 


+25:02:52.29 


19.0 ± 1.0 


-47.1 ± 1.0 


133.0 ± 6.0 


9.2 


G5 


Pleiades 


8.1 


HII 2106 


03:48:58.49 


+23:12:04.33 


16.5 ± 1.3 


-44.9 ± 1.1 


133.0 ± 6.0 


9.4 


Kl 


Pleiades 


8.1 


RX J0348.9+0110 


03:48:58.76 


+01:10:53.99 


35.1 ± 1.6 


-22.1 ± 1.2 


100.0 ± 50.0 


8.3 


K3(V)/E 




8.2 


HII 2278 


03:49:25.70 


+24:56:15.43 


18.4 ± 0.9 


-47.0 ± 0.8 


133.0 ± 6.0 


8.8 


K0 


Pleiades 


8.1 


HII 2506 


03:49:56.49 


+23:13:07.01 


17.6 ± 0.7 


-43.9 ± 0.6 


133.0 ± 6.0 


8.8 


F9 


Pleiades 


8.1 


HII 2644 


03:50:20.90 


+24:28:00.22 


19.8 ± 0.8 


-46.8 ± 0.9 


133.0 ± 6.0 


9.3 


G5 


Pleiades 


8.1 


HII 2786 


03:50:40.08 


+23:55:58.94 


17.6 ± 0.7 


-45.2 ± 1.0 


133.0 + 6.0 


8.9 


F9 


Pleiades 


8.1 


HII 2881 


03:50:54.32 


+23:50:05.52 


17.7 ± 0.7 


-46.9 ± 1.1 


133.0 + 6.0 


9.1 


K2 


Pleiades 


8.1 


HII 3097 


03:51:40.44 


+24:58:59.41 


17.5 ± 0.7 


-46.1 ± 1.0 


133.0 ± 6.0 


9.1 


G6 


Pleiades 


8.1 


HII 3179 


03:51:56.86 


+23:54:06.98 


19.2 ± 0.6 


-46.5 ± 0.8 


133.0 ± 6.0 


8.6 


F8 


Pleiades 


8.1 


HD 285281 


04:00:31.07 


+19:35:20.70 


2.7 ± 1.1 


-12.9 ± 1.2 


49.0 ± 11.0 


7.6 


Kl 




7.0 


HD 284135 


04:05:40.58 


+22:48:12.14 


6.0 ± 0.6 


-14.9 ± 0.6 


140.0 ± 70.0 


7.8 


G3(V) 




6.8 


HD 281691 


04:09:09.74 


+29:01:30.55 


19.9 ± 0.7 


-36.3 ± 1.0 


140.0 ± 70.0 


8.4 


K1(V) 




7.8 


HD 26182 


04:10:04.69 


+36:39:12.14 


23.8 ± 0.7 


-36.7 + 0.7 


100.0 ± 50.0 


7.8 


GOV 




7.8 


HD 284266 


04:15:22.92 


+20:44:16.93 


1.8 ± 1.0 


-13.6 ± 0.7 


140.0 ± 70.0 


8.6 


K0(V) 




7.3 


HD 26990 


04:16:16.50 


+07:09:34.15 


-85.6 ± 1.5 


-52.1 ± 1.5 


35.0 ± 2.0 


5.9 


G0(V) 




8.9 


HD 27466 


04:19:57.08 


-04:26:19.60 


-58.6 ± 1.2 


-37.0 ± 1.2 


36.0 + 1.0 


6.3 


G5V 




9.2 


vB 39 


04:22:44.74 


+16:47:27.56 


173.3 ± 11.5 


4.7 ± 10.2 


39.3 ± 3.5 


6.2 


G4V 


Hyades 


8.8 


vB 49 


04:24:12.78 


+16:22:44.22 


87.6 + 1.3 


-21.9 + 1.2 


57.5 + 1.0 


6.8 


GOV 


Hyades 


8.8 


vB 52 


04:24:28.33 


+16:53:10.32 


113.1 ± 1.4 


-23.3 ± 1.2 


44.8 ± 0.8 


6.3 


G2V 


Hyades 


8.8 


vB 176 


04:25:47.56 


+18:01:02.20 


102.6 ± 2.2 


-29.9 ± 3.2 


48.0 ± 1.0 


6.8 


K2V 


Hyades 


8.8 


vB 63 


04:26:24.61 


+16:51:11.84 


106.7 ± 1.3 


-24.5 ± 1.2 


46.9 ± 1.0 


6.4 


G1V 


Hyades 


8.8 


vB 64 


04:26:40.11 


+16:44:48.78 


107.0 ± 1.1 


-26.8 ± 1.1 


46.4 ± 0.9 


6.5 


G2+ 


Hyades 


8.8 


vB 66 


04:27:46.07 


+ 11:44:11.07 


110.1 ± 1.3 


-13.2 ± 1.2 


44.6 ± 0.9 


6.2 


F8 


Hyades 


8.8 


7^ 




il 7.1 7-07 QA 


110 1 i 1 1 


OQ Q _l_ 1 n 


AA fi + n8 


A 




Hyades 


0.0 


vB 79 


04:29:31.61 


+17:53:35.46 


106.7 + 1.1 


-31.4 + 1.1 


45.6 + 0.8 


7.1 


K0V 


Hyades 


8.8 


vB 180 


04:29:57.73 


+16:40:22.23 


106.2 ± 1.1 


-27.1 ± 1.1 


46.0 + 0.8 


7.1 


K1V 


Hyades 


8.8 


vB 88 


04:31:29.35 


+13:54:12.55 


90.0 + 1.2 


-16.0 + 1.2 


53.1 ± 1.3 


6.5 


F9V 


Hyades 


8.8 


1RXS J043243.2-152003 


04:32:43.51 


-15:20:11.39 


2.3 ± 1.1 


14.2 + 1.1 


140.0 ± 70.0 


8.6 


G4V 




6.6 


vB 91 


04:32:50.12 


+16:00:20.96 


103.2 ± 1.0 


-25.9 + 1.0 


45.9 + 0.6 


6.8 


G7 


Hyades 


8.8 


vB 92 


04:32:59.45 


+15:49:08.37 


99.1 ± 1.2 


-24.1 ± 1.2 


47.8 ± 0.8 


6.9 


G7 


Hyades 


8.8 


vB 93 


04:33:37.97 


+16:45:44.96 


99.0 + 1.1 


-22.9 + 1.2 


48.3 ± 0.7 


7.4 


G7 


Hyades 


8.8 


vB 96 


04:33:58.54 


+15:09:49.04 


101.9 ± 1.3 


-29.4 + 1.3 


45.4 + 0.8 


6.5 


G5 


Hyades 


8.8 


RX J0434.3+0226 


04:34:19.54 


+02:26:26.10 


18.0 ± 2.0 


-16.4 + 1.9 


161.0 + 24.0 


9.5 


K4e 




7.8 


vB 183 


04:34:32.18 


+15:49:39.23 


91.0 + 1.0 


-20.0 + 1.0 


51.7 + 0.8 


7.6 


G7 


Hyades 


8.8 


vB 97 


04:34:35.31 


+15:30:16.56 


98.1 ± 1.0 


-26.7+1.1 


47.2 + 0.9 


6.4 


F8:V: 


Hyades 


8.8 


vB 99 


04:36:05.27 


+15:41:02.60 


95.0 + 1.0 


-23.1 ± 1.2 


48.7 + 0.7 


7.4 


G7 


Hyades 


8.8 
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TABLE 3 

Median Sample Statistics 





log (Age/yr) 


Distance (pc) 


Spectral Type 


Sample 


range median 


range median 


range median 


Deep 


6.6-8.7 8.0 


10-190 50 


F5-K5 G5 


Shallow 


6.6-10.0 8.8 


11-190 45 


F5-K5 G7 


Complete 


6.6-10.0 8.3 


10-190 46 


F5-K5 G5 



TABLE 4 

Calibration Binaries and Assumed Parameters of Their Astrometric Orbits 



Binary 
(WDS) 


a 

(arcsec) 


P 

(years) 


To 
(year) 


e 


i 

(degrees) 


n 

(degrees) 


(degrees) 


grade 


Ref. 


09006+4147 


0.6472 + 0.0010 


21.776 + 0.017 


1993.725 ± 0.023 


0.1507 + 0.0008 


131.26 + 0.13 


204.39 + 0.19 


32.52 + 0.36 


1 


1 


16147+3352 


5.927 


888.989 


1826.949 


0.7605 


31.795 


16.889 


72.201 


i 


2 


18055+0230 


4.5540 ± 0.0052 


88.38 + 0.02 


1895.94 + 0.02 


0.4992 ± 0.0004 


121.16 + 0.08 


302.12 + 0.10 


14.0 + 0.1 


1 


3 


20467+1607 


10.22 


3249 


2305 


0.88 


148.78 


88.06 


331.16 


4 


4 
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TABLE 5 

Extinction Due to PHARO and NIRC2 Optics 





AJ 


AH 


AK S 


Transmissive Optic 


(mag) 


(mag) 


(mag) 


PHARO ND 1% filter 


4.753 ± 0.039 


4.424 + 0.033 


4.197 + 0.024 


NIRC2 1" coronagraph 


8.36 + 0.28 


7.78 + 0.15 


7.10 + 0.17 


NIRC2 2" coronagraph 


9.26 ± 0.09 


7.79 + 0.22 


7.07 + 0.22 
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TABLE 6 

Deep Sample Observations and Sensitivity Limits 





Date 






Limiting A"c-Band Magnitude at b 


Star 


(UT) 


Observatory 


Optic a 


1" 


2" 


5" 


HD 377 


2002 Aug 28 


Palomar 


corona 










2003 Nov 09 


Keck 


corona 


16.7 


18.3 


20.0 




2004 Oct 07 


Keck 


corona 








HD 691 


2002 Aug 31 


Palomar 


corona 










2002 Nov 18 


Palomar 


corona 


15.1 


18.3 


20.3 


HD 984 


2002 Aug 29 


Palomar 


corona 


13.1 


16.9 


18.5 




2003 Sep 20 


Palomar 


corona 










2004 Oct 04 


Palomar 


corona 








HD 1405 


2004 Jun 06 


Keck 


corona 


16.3 


18.3 


19.2 


QT And 


2002 Aug 29 


Palomar 


corona 










2003 Nov 10 


Keck 


corona 


17.6 


19.1 


19.9 




2004 Oct 07 


Keck 


corona 








HD 7661 


2002 Aug 30 


Palomar 


corona 










2002 Nov 09 


Keck 


corona 


16.5 


18.0 


20.0 


HIP 6276 


2002 Aug 30 


Palomar 


corona 


14.5 


18.9 


20.7 




2004 Oct 08 


Keck 


corona 








HD 8907 


2002 Aug 27 


Palomar 


corona 


13.5 


17.3 


19.3 


HD 12039 


2002 Aug 28 


Palomar 


corona 


14.0 


18.5 


20.0 


HD 15526 


2003 Sep 20 


Palomar 


corona 


15.5 


19.4 


20.3 


1RXS J025216.9+361658 


2002 Nov 18 


Palomar 


corona 










2003 Nov 10 


Keck 


corona 


17.8 


19.3 


20.3 




2004 Feb 07 


Palomar 


corona 








HD 17925 


2003 Sep 21 


Palomar 


corona 


12.5 


15.3 


18.6 


1RXS J025751. 8+115759 


2002 Nov 16 


Palomar 


corona 


16.4 


19.8 


20.5 


RX J0258.4+2947 


2002 Feb 28 


Palomar 


corona 


17.7 


19.5 


19.6 


1RXS J030759. 1+302032 


2002 Nov 18 


Palomar 


corona 


15.1 


18.2 


19.8 


HD 19668 


2002 Aug 27 


Palomar 


corona 










2003 Nov 09 


Keck 


corona 


17.0 


18.5 


20.0 


IE 0307.4+1424 


2003 Jan 12 


Palomar 


corona 


15.6 


18.9 


20.5 


V525 Per 


2003 Sep 20 


Palomar 












2004 Oct 03 


Palomar 












2004 Oct 08 


Keck 


corona 


18.7 


20.3 


20.4 


1RXS J031907.4+393418 


2002 Aug 29 


Palomar 


corona 










2003 Nov 10 


Keck 


corona 


19.3 


20.5 


20.7 




2004 Oct 08 


Keck 


corona 








HE 622 


2003 Sep 20 


Palomar 












2004 Oct 08 


Keck 


corona 


19.2 


20.8 


20.9 


IE 0324.1-2012 


2003 Jan 12 


Palomar 


corona 










2004 Feb 07 


Palomar 


corona 










2004 Oct 08 


Keck 


corona 


18.1 


19.5 


19.8 


RX J0329. 1+0118 


2003 Sep 21 


Palomar 


corona 


16.6 


19.7 


20.1 


HE 1101 


2003 Sep 20 


Palomar 












2004 Oct 03 


Palomar 












2004 Oct 07 


Keck 


corona 


19.2 


20.9 


21.0 




2006 Dec 12 


Keck 


corona 








HD 22179 


2002 Nov 16 


Palomar 


corona 










2003 Nov 09 


Keck 


corona 


18.1 


19.4 


20.2 




2004 Feb 07 


Palomar 


corona 










2004 Oct 08 


Keck 


corona 








HD 23208 


2004 Oct 05 


Palomar 


corona 


14.9 


19.2 


20.1 


HII 120 


2003 Sep 20 


Palomar 












2004 Oct 04 


Palomar 












2004 Oct 08 


Keck 


corona 


18.6 


20.0 


20.2 




2006 Jan 18 


Palomar 










HII 2147 


2003 Jan 13 


Palomar 


corona 


15.2 


19.3 


20.6 


1RXS J035028.0+163121 


2002 Nov 17 


Palomar 


corona 


16.7 


19.9 


20.3 


RX J0354.4+0535 


2003 Jan 13 


Palomar 


corona 










2004 Feb 07 


Palomar 












2004 Oct 07 


Keck 


corona 


18.3 


19.9 


20.3 


HD 283167 


2002 Nov 18 


Palomar 


corona 


17.0 


19.6 


19.8 


RX J0357.3+1258 


2003 Jan 11 


Palomar 


corona 










2003 Nov 09 


Keck 


corona 


18.9 


20.3 


20.5 




2004 Feb 07 


Palomar 


corona 








HD 285751 


2002 Jan 31 


Palomar 


corona 


15.7 


18.1 


19.2 


RX J0442. 5+0906 


2003 Jan 11 


Palomar 


corona 


16.4 


19.5 


20.0 




2004 Feb 07 


Palomar 










HD 286179 


2002 Jan 31 


Palomar 


corona 










2004 Feb 07 


Palomar 












2004 Oct 07 


Keck 


corona 


17.8 


19.6 


20.3 




2006 Dec 12 


Keck 


corona 








HD 31950 


2002 Nov 16 


Palomar 


corona 










2003 Nov 09 


Keck 


corona 


18.1 


19.4 


19.8 




2004 Feb 07 


Palomar 


corona 










2004 Oct 08 


Keck 


corona 








HD 35850 


2002 Feb 01 


Palomar 


corona 


11.7 


14.3 


16.9 


1RXS J053650.0+133756 


2002 Feb 28 


Palomar 


corona 










2003 Jan 14 


Palomar 


corona 










2003 Nov 10 


Keck 


corona 


18.4 


18.0 


20.6 
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TABLE 7 

Shallow Sample Observations and Sensitivity Limits 



Star 


Date 


Observatory 


Optic a 


Limiting A"s-Band Magnitude at b 


(UT) 


U. O 


1 


2 


w' 
o 


1LLJ ZZIO 1 O 


9009 Alio- 31 


P a 1 nrn q v 
L cLlUlllcll 


NDl 












9003 91 
ZUUO >0"P 


PciloLTiax 


NDl 


11.8 


13.5 


10.0 


14.7 


HD 6963 


2004 Jun 04 


Palomar 


NDl 


10.3 


11.9 


13.8 


14.1 


11U OrfcU I 


9004 Tim 04 

ZUUrl: O Llll Urfc 


P a 1 r"i vy i o t- 

L tl.lUlllclil 


NDl 


10.8 


12.4 


13.6 


13.9 


HF) 8Q41 


9004 Tun 04 


P £l 1 t~l W"l Q "P 

IT Cl.lUlllcll 


NDl 


9.8 


11.3 


13.5 


14.0 


tin Q4.79 

11U C/'rt ( Z 


9009 Nnv 1 8 

ZUUZ 1>UV lO 

2003 Sen 20 

ZL/UiJ kj" kJ ZrU 


P ^ 1 t~l W"l Q "V 

IT ellUUlcll 
P pi 1 t~i m ci v 

1 CLlLJllldl 


NDl 
NDl 












9004 Ort 0^ 
ZUU'i WCLr UO 


Palomar 


NDl 


10.8 


12.6 


15.1 


15.8 


RE T0137-H8A 


9002 Tan 31 

ij U U ii O elll v> 1 

9003 ^pn 90 


P a 1 t~i m n v 

1 CLlLJllldl 

P p* 1 ntn a v 

IT fllUlllcll 


NDl 


11.7 


13.2 


12.0 


14.9 


un 1 1 Qccn 

11.L/ llOUU 


9003 ^pn 91 

ZUUO JtU Zl 


P Q 1 TYl CI T" 
L Cl.lUlllcll 


NDl 


10.7 


12.4 


14.8 


15.3 


un i 3389 


9004 Tim 04 

ZUUrl: O Llll U"i 


P Q 1 t~l VY"1 O "P 

IT ellUUlcll 


NDl 


10.4 


12.2 


14.7 


15.0 


HD 1 3^07 


9009 Alio- 98 


P p* 1 t~i vn T" 

IT dlUlllcll 


corona 


11.7 


13.8 


15.6 


18.1 




9009 Alio- 98 
9003 ^pn 91 

ZUUO OcL) Zl 


P a 1 c^ vn v 
IT cllUlIlcll 

P Q 1 C^ 1 11 O "V 

IT cllUlIlcll 


NDl 
NDl 












9004 Opf 0^ 


Palomar 


NDl 


10.4 


12.1 


14.7 


15.6 


HD 13974 


2004 Jun 04 


Palomar 


NDl 


7.6 


9.2 


11.8 


13.5 


IRY^l T09^993 ^4-379014 


9003 91 

zuuo oep zi 


Palomar 














9004 Ort O 1 ^ 

ZUU^i Utl uo 


P C* 1 t~l VY"1 O "P 
IT Cl.lUlllcll 




14.1 


12.8 


17.6 


18.2 


2RF 102^+474 

Zl Vi-J U UZ OO n^^t 1 ^ 


9009 Fob 28 

iUUiJ 1 LU ZO 


Pa 1 nrn civ 
i ci.iu iiicii 














9004 Ort O 1 ^ 


Palomar 




12.5 


14.9 


15.8 


16.7 


HF) 1 8Q40 

11J_^ lu.7 1U 


9009 Alio- 9Q 

iUUi Jt. Llg Z ,y 


T-^ ci 1 a v 

IT Cl.lUlllcll 1 


NDl 


9.9 


11.9 


14.8 


15.1 


un i am a 


9004 Tim 04 
ZUU4: J Llll Kj^i 


Palomar 


NDl 


9.9 


12.0 


14.6 


14.9 


HD 1 QR39 

11.L/ l;?UOZ 


9009 Alio- 30 

ZUUZ I\ Llg, OU 


P Q 1 t~l VY 1 d T" 
IT Cl.lUlllcll 1 




9.8 


12.0 


14.2 


14.5 


vB 1 


9009 Alio- 9Q 

ZUUZ -£l Llg Z i7 


P C* 1 C^ VY"1 O "V 
IT CllUlIlcll 


NDl 


10.3 


12.2 


14.6 


14.8 


HP 3^0 


9003 ^pn 90 

ZUUO OcJU ZU 


P Q 1 C^ TT1 CI T" 
IT CllUlIlcll 














2004 Oct 04 


Palomar 




13.7 


15.9 


18.1 


18.3 


HE 373 


9003 ^pn 90 

ZUUO Ot,JU ZU 

9003 Nnv 1 

ZUUO 1>UV 1U 


P C* 1 C^ 1Y"1 Q V 
L CllUlllCll 














9004 Opt O^i 
ZUU^i ULL UO 


P Q 1 t~l TYl CI T" 
IT dlUlllcll 




13.8 


16.1 


18.0 


18.2 


HP 38Q 

lll_j OOj 


9003 ^pn 90 

ZUUO UCL) ZrU 


P C* 1 C^ TT1 CI T" 

IT cllUlIlcll 














9004 Opt 04 

ZUU^i ULL U'i 


P Q 1 C^ VV1 CI V 

IT dlUlllcll 




14.3 


16.4 


18.3 


18.5 


HP fiQfi 


9003 ^pn 90 

ZUUO OvrJU ZrU 


ST dlUlllcll 




14.4 


15.6 


17.5 


17.7 


HP fiQQ 


9003 ^pt^ 90 
zuuo oep zu 


Palomar 




14.1 


16.2 


17.5 


17.7 


HE 750 
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Candidate Companions in the Deep Sample 
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zUUz-ll-lb 


1 




Z 


Q 

o 


iok ± n nn7 
loO ± U.UU* 


1 OQ QQ _L n 

lyo.oo ± u 


1 7 
1 / 


D 


a a in 
44 ± U 


O/l 
z 1 


i -1 

14 


nQ _i_ n 
Uo ± U 


OA 

z*± 


n c; /i _i_ n q a 
U.04 ± U.o4 


onno 11 ift 
zUUz-ll-lb 


1 




3 


6 


7yio j_ n no a 
(4o ± U.U24 


o i ft i o _L n 
olb.lo ± U 


22 


8 


28 ± 


24 


15 


87 ± 


24 


n n7 _i_ n o /i 
U.y / ± U.o4 


onno 11 ift 
zUUz-ll-lb 


1 




1 


1 n 
1U 


oo*7 _L n no/i 
yz / ± U.UZ4 


qi c; ftQ j_ n 
olO.bo ± U 


1 7 
1 ( 


D 


01 J_ n 
zl ± U 


OA 

z4 


1 Q 

lo 


on _i_ n 
oU ± U 


O -1 

Z4 


n c; k j_ n q /i 
U.OO ± U.o4 


onno 11 ift 
zUUz-ll-lb 


1 




5 


2 


Tod J_ n nn7 
/24 ± U.UU ( 


OI 07 _L n 
zl.o7 ± U 


13 


11 


55 ± 


24 


19 


14 ± 


O -1 

z4 




onno 1 1 nn 
zUUo-ll-Uy 


IV 


G A f~\ 1 7,nft7ft 
oAU lOUb/b 


i 
i 


o 
o 


Q7cr j_ n noo 

o to ± u.uzy 


q c 1 q 1 i n 
oOl.ol ± U 


1 A 
14 


y 


Qn J_ n 
oU ± U 


on 
zU 


1 ft 
ID 


77 _i_ n 
/ / ± U 


on 
ZU 




onno 11 17 
zUUZ-11-1 / 


1 


HIJ 4yiy/ 


1 


6 


ncn _i_ n ni ft 

yoz ± U.Ulb 


o A c; o o _L n 
o4o.OZ ± U 


18 


6 


7K i n 
zo ± U 


06 


12 


82 ± 


06 


n i k _L n i o 
U.lo ± U.lz 


onno no oq 
zUUZ-Uz-zo 


TD 
1 




Z 


n 
U 


CiA Q _L n flQO 

y4o ± U.Uoz 


/ / .OU ± 1 


n q 
Uo 


Q 

o 


oo _L n 
zz ± U 


1 A 
14 


1 A 
14 


oo _i_ n 
zy ± U 


1 1 
14 


1 ftQ _l_ 1 01 
l.bo ± L.ZL 


onno no oq 
zUUz-Uz-zo 


1 ) 
1 


"D T7 1 Tn70Q i on 
KHj JU ( Zo+zU 


1 


8 


i nft _i_ n ni o 
lyb ± U.Ulo 


on oft _i_ n 
oU.ob ± U 


03 


7 


80 ± 


20 


14 


68 ± 


20 


n i ft _i_ n oo 
U.lb ± U.zz 


onno no oq 
zUUZ-Uz-zo 


1 




Z 


O 


crqo _L n ni Q 

Ooz ± U.Ulo 


QOO Qft _L n 

ozy.ob ± u 


no 

uy 


Q 

o 


/in j_ n 
4U ± U 


on 
zU 


lo 


oq _i_ n 
zo ± U 


on 
ZU 


1 nft _i_ n oo 
l.Ub ± U.zz 


onno no oq 
ZvuZ-UZ-Zo 


1 


HD bU/oY 


1 


7 


ftc^T _L n non 
bo/ ± U.Uzy 


1 07 OK j_ n 

Iz / .zo ± U 


"1 o 

lo 


9 


40 ± 


20 


15 


65 ± 


20 




onno m qi 
zUUz-Ul-ol 


1 


nn CO/1 /I Q 
rilJ oz44o 


1 


o 


a c^o _L n ni n 

4oy ± u.uiu 


1 on qn _i_ n 

lyu.ou ± U 


O Q 
ZO 


i 1 
ii 


77 _i_ n 


1 A 
14 


i ft 
ID 


oo _l n 

oy ± u 


1 A 
14 




onn/i no n7 
zUU4-Uz-U / 


1 




2 


8 


i c j j_ n non 
154 ± U.UzU 


no 7ft _i_ n 

yo. to ± U 


15 


12 


59 ± 


21 


"1 7 

1 7 


71 J_ n 
71 ± U 


21 




onn/i no 07 
2UU4-U2-U/ 


1 




Q 
O 


I 


i a o j_ n no7 
14z ± U.Uz i 


OC^Q 71 _L n 

ZOo. / 1 ± U 


O Q 
ZO 


1 Q 

lo 


q a in 
o4 ± U 


Qn 
oU 


1 Q 

lo 


oft _i_ n 

yb ± u 


Qn 
oU 




onn/i no 07 
zUU4-Uz-U / 


T") 
1 


G A f~\ 1 70070 

oAU 1/oZ/Z 


1 


1 n 
1U 


Uoz ± U.Uv5z 


q c^ft ft 1 in 
v>0b.b4 ± U 


"1 Q 

lo 


y 


ft7 _L n 
o ( ± U 


1 c 
10 


"1 7 
1 ( 


nft _i_ n 
Ub ± U 


lo 




onnQ m 1 q 
zUUo-Ul-lo 


1 ) 
1 




2 


8 


i o /i _L n n /i ft 
184 ± U.U4b 


274.53 ± 


15 


10 


75 ± 


22 


18 


14 ± 


22 




onno m 1 o 
zUUo-Ul-lo 


1 > 
1 


xi ij yuyuo 


i 
l 


O 


q 1 ft _L n no7 
olb ± U.Uz i 


1 01 77 _L n 

iyi. / / ± u 


O Q 
ZO 


1 1 
ii 


Qn J_ n 
oU ± U 


i n 
1U 


i ft 
lb 


go _i_ n 
oz ± U 


i n 
1U 




onno no m 
zUUz-UZ-Ul 


1 




2 


12 


a Ad _L n nQ 1 
44b ± U.Uol 


1 7ft 7Q _L n 

1 (b. zo ± U 


13 


13 


49 ± 


19 


19 


01 ± 


19 




onn/i nft nc 
zUU4-Ub-Uo 


IV 


rrn m 7qo 

ri ij yi ioz 


l 


l 


nno J_ n nno 
UUz ± U.UUo 


QQ ft7 _L n 


/1ft 

lb 


A 
1 


Qn _l n 
oU ± U 


nft 
Uo 


i 1 
ll 


no _i_ n 
Uo ± U 


nft 
Ub 


n on _i_ n i q 
u.yu ± U. lo 


onno nQ no 
ZvuZ-Uo-uZ 


1 


TTn nooKt; 

lll) yzooo 


i 
l 


Z 


OQ/i J_ n nnt; 
yo4 ± U.UUO 


OOI QQ _l n 

zyi.oo ± u 


"1 Q 

lo 


A 
1 


K7 i n 
/ ± U 


no 

uy 


i n 
1U 


/ip in 
4b ± U 


no 

uy 


n 7C i n i 7 
U. / ± U. 10 


onno no m 
zUUz-Uz-Ul 


1 




Z 


1 O 
IZ 


oi ft _l n noo 
zlb ± U.Uzz 


1/17 70 J_ O 

14 / . / y ± u 


O K 

zo 


o 
o 


on J_ n 
yu ± U 


on 
zU 


i -1 
14 


70 J_ n 

/y ± u 


on 
ZU 




onno no m 
zUUZ-Uz-Ul 


1 


Lrlc^ IjCO 


1 


U 


oa q _L n nno 
z4o ± U.UUz 


07Q oo _L n 
z / o.zz ± U 


"1 1 
1 1 


n 
U 


-i o i n 
lo ± U 


nft 
Ub 


o 
o 


ceo _i_ n 
Oo ± U 


nft 
Ub 


1 no _i_ n no 
l.Uz ± U.Uo 


onnQ m 1 
zUUo-Ul-lZ 


1 




Z 


1 n 
1U 


noo j_ n nno 

uoo ± u.uuy 


QOC ftC _L n 

ozo.bO ± u 


no 

uy 


D 


/in i n 
4U ± 


nft 
Ub 


1 A 
14 


ip in 
lb ± U 


nft 
Ub 


n iq i n no 

u. io ± u.uy 


onnQ m 1 
zUUo-Ul-lZ 


1 


HD 104576 


1 


10 


455 ± 0.028 


19.66 ±0 


21 


11 


00 ±0 


50 


17 


68 ±0 


50 




2002-06-22 


1' 


HD 104860 


1 


3 


803 ± 0.027 


287.01 ± 


28 


10 


92 ±0 


25 


17 


42 ±0 


25 


0.00 ±0.47 


2002-06-23 


P 




2 


11 


961 ± 0.033 


260.09 ± 


19 


12 


09 ±0 


18 


18 


59 ±0 


18 




2004-06-05 


K 


SAO 15880 


1 


2 


176 ±0.018 


293.93 ± 


72 


8 


98 ±0 


17 


16 


27 ±0 


17 




2004-02-06 


P 


HD 111456 


1 


3 


783 ± 0.010 


117.45 ±0 


30 


12 


72 ±0 


16 


17 


27 ±0 


16 




2004-02-06 


P 


HD 139498 


1 





311 ±0.002 


3.39 ±0 


21 





00 ±0 


02 


8 


26 ±0 


03 


0.50 ±0.05 


2003-07-15 


P 




2 


11 


246 ± 0.033 


123.98 ± 


19 


8 


48 ±0 


30 


15 


98 ±0 


30 




2004-06-26 


P 




3 


8 


801 ± 0.026 


61.50 ±0 


21 


10 


98 ±0 


30 


18 


49 ±0 


30 




2004-06-26 


1' 


HD 142361 


1 





705 ± 0.001 


236.41 ± 


13 


2 


01 ±0 


10 


9 


19 ±0 


10 


0.85 ±0.14 


2002-06-21 


1' 




2 


11 


207 ± 0.046 


164.99 ± 


17 


5 


85 ±0 


17 


12 


88 ±0 


17 


0.77 ±0.28 


2002-06-21 


P 


HD 143006 


1 


8 


355 ± 0.026 


130.27 ± 


25 


9 


28 ±0 


16 


16 


33 ±0 


16 


1.18 ±0.23 


2002-06-23 


P 




2 


6 


626 ± 0.028 


0.32 ± 


23 


10 


40 ±0 


16 


17 


45 ±0 


16 


1.40 ±0.23 


2002-06-23 


1' 




3 


8 


502 ± 0.029 


268.41 ± 


23 


10 


66 ±0 


16 


17 


71 ±0 


16 


1.27 ±0.23 


2002-06-23 


P 




4 


7 


698 ± 0.023 


357.97 ±0 


12 


12 


11±0 


16 


19 


16 ±0 


16 


1.30 ±0.23 


2003-05-18 


K 



no (a) 
yes(c) 
no (a) 
yes(c) 
no (a) 
no (a) 
7 
? 

no (a) 
no (a) 

? 

? 

? 

? 

? 

yes(c) 
no (a) 
no (a) 
no (a) 
no (a) 
no (a) 
no (a) 
no (a) 
no (a) 
no (a) 
no (a) 
no(c) 

yes(c) 
no (a) 
7 

no (a) 
no(c) 
yes(a) 
no(c) 
? 

no(c) 
? 



no (a) 
no (a) 
yes(c) 
no (a) 
no (a) 
no (a) 
no (a) 
no (a) 
no (a) 
yes(a) 
no (a) 
no (a) 
no (a) 

? 

? 

? 

? 

? 

no (a) 
no(c) 
yes(a) 
yes(a) 
no (a) 
yes(a) 
no (a) 
no(e) 
no (a) 
no(c) 
no (a) 
7 

yes (a) 
? 
? 

yes(a) 
no (a) 
7 
7 

no (a) 
? 



MH04 



FM00 



WDS 



G93 



32 



Metchev & Hillenbrand 



TABLE 9 

Candidate Companions in the Shallow Sample 







P 


9 


AK S 


K s 


J - K s 










Primary Star 


N c 


(arcsec) 


(deg) 


(mag) 


(mag) 


(mag) 


(UT Date) 


Tel. 


Assoc. 


Ref. 


HD 224873 

illy _ i - ) 


I 


i ofiS -1- 002 

1 . Zj U O _l_ \J .\)\J jLj 


171 44 + 1 2 

A. 1 A-.TTTT 1 Ui 1^ 


25 + 02 

U.ZjO _1_ U.UZi 


7 K7 + 03 

1 . O 1 _l_ U. UO 


Kn + 05 

U.OU _l_ U.UO 


2002-08-31 

ZjUUZjUOOI 


p 


yes ^a j 


WDS 


HD Q472 


I 


2 703 + n n25 

Zi . 1 CO _l_ U.UZjO 


343 fiq + 30 

'Jt:'J . UC _1_ U.JU 


k 7q + 09 

O. 1 C _l_ U.UC 


11 83 + 09 

11.00 _l_ U.UC 


1.04 ± 0.14 


2002-1 1 -1 8 

ZjUUZj 1 1 10 


p 


yes(a) 




RE 101 37+1 8 A 

A LJ-J iA U A O 1 |^ J. Oai. 


I 


1 «qi + n nnfi 

A.UtJA _l_ U • UUU 


24 fiq + 1 2 

iirr.VAU 1 V7 - 1 Z. 


05 + 01 

U.UO _l_ U.UA. 


7 4q + 02 

1 .TiC _l_ U.UZi 


7R + 03 

U. 1 U _l_ U.UO 


2002-01 -31 

ZjUUZjUI OI 


p 


yes (a) 


WDS 


HE 13531 

A A I ' IfJUrJl 


I 


71 7 + 003 

u. 1 a. 1 _i_ u . uuo 


1 k 7q + 43 

A_U. 1 f _1_ U.rtiJ 


4 20 + 08 

rl.ZiU _1_ U.UO 


q oa + 08 

COO _1_ U.UO 


1 04 + 1 5 

l.Urt _l_ U.IO 


2002-08-28 

ZjUUZjUOZjO 


p 


yes(a) 




IRY^l Tfi9^99^ ^4-379Q1 4 

±1 l__/vO U U^iJiiO. UyO i ZjC Art: 


I 


U.UO 1 □__ U-UUO 


01 90 _|_ n 90 


1 Ar> _|_ n no 


in 77 -1- OS 


n f\7 -u n 19 

U.u 1 1T_ U. 1-i 


90(1^ HQ 91 
-jUuo uy j-i 


p 


yes (a) 






2 


O . Zj U U _l_ U.U AU 


7« or: 1 q 1 

1 U.OU _l_ U.A.O 


t: . O 1 _1_ U.UC 


lO.rtO _1_ U.UC 


n cn _|_ Q 1 

U.OU _1_ U.IO 


9on i >-nQ-9i 

ZjUUO UC Zj 1 


p 


no (a) 




2RE T0255+474 

ill til 7 O l/ilUU | ^ 1 




2131 + 004 

i:. AO A _1_ U.UU^ 


272 63 + 16 


08 + 05 

U . UO _l_ U . UO 


7 20 + 06 

1 .ZjC _1_ U.UU 


72 + 08 

U. 1 Zi _1_ U.UO 


2002-02-28 

ZjUUZjUZjZjO 


p 


yes la j 


WDS 




2 


1 1 4«q + n 033 

ALIUii _1_ U.UOO 


4« 40 + 1 1 


7 00 + 1 

1 . UU _1_ U.AU 


1 4 21 +010 

Art.Zjl 1 U.1U 


k.7 + 1 4 

U.Ol _l_ U.lrt 


2002-02-28 

ZjUUZj UZj ZjO 


p 


no (a) 




HE 18940 

A A I ' IOt/tU 




1 fi7 + 002 

U . 1U 1 _1_ U.UUi: 


a Kq + 1 is 

_1_ A.A.CJ 


7a + 03 

U . 1 O _l_ U . UO 


ft 71 + 04 

U. 1 1 _1_ U.Urt 


Ki + 08 

U.Ol _l_ U.UO 


2002-08-29 

ZjUUZj WI. i — Z -7 


p 


y es l c i 


Hip 




2 


4 321 +|)(iio 

1»(JZ1 1 U ♦ U A Z 


207 38 + 1 2 


4 cs + 03 

rl.OO _1_ U.UO 


1 08 + 04 

1U.UO _l_ U . Urt 


02 + 07 

U.CZi _l_ u.u 1 


2002-08-29 

ZjUUZj WO Zj.V 


p 






3 


4 1 20 + 01 


203 78 + 1 3 


k 21 + 03 

U.ZiA. _1_ U.UO 


1 71 + 04 

1U. 1 1 _1_ U.Urt 


21 + 07 

U.Ol _l_ u.u 1 


2002-08-29 


p 


. X 




vB 1 




2 47n 4- n n06 

i.T 1 U _1_ U . UUU 


200 63 + 1 4 


2 K3 + 03 

Z. . UO _l_ U.UO 


a Ro + 04 

O . UZj _1_ U.Urt 


qi + 06 

U.C1 _1_ U.UU 


2002-08-29 

ZjUUZjUOZjC 


p 


yes^c j 


WDS 


HE 350 

AAA_J OOU 




8 464 4-0 016 


1 0q 22 + 14 

1U7 . £i £i I U • A-TT 


K OK _l_ 21 

O.OO _l_ U . Zj A. 


1; 11 + 21 

lO.ll 1 U • Zj 1 


1 50 + 37 

l.OU _1_ U.Ol 


2003-09-20 

ZiUUO UC ZiU 


p 


no (a) 






2 


« oqfi 4. 01 1 

U.UJU _l_ U.U11 


3a 37 + 1 q 


7 «« + 30 

1 . UU _l_ U.OU 


1 a no + 30 

IU.CZj _1_ U.OU 




2004-1 0-04 

ZjUUtt 1U Utt 


p 




HE 373 


1 


2 081 4- 005 

iiUui i_ u - \ j u u 


193 77 + 18 

a_co. 1 f _1_ \J*-Lyj 


5 24 + 1 

',7 . — 1 1 U.AU 


1 4 59 + 1 

ATT. OC _1_ U.AU 


98 + 14 

U.CO _L_ U.l^t 


2003-09-20 

ZjUUO UC ZjU 


p 




MH04 




2 


1 1 cno _i_ n 031 

ll.Ui/O _l_ U.UU1 


2fic: 01 1 25 


7 ki + 30 

1 .OA. _l_ U.OU 


1 « a« + 30 

AU.OU _l_ U.OU 




2003-09-20 

Zj u uo u c Zj yj 


p 






3 


2 47c 4. 034 


55.82 ± 0.22 


a 37 + 30 

O . O 1 —i U . OU 


1 7 72 + 30 

-L 1 m i £j _1_ U.OU 




2003-09-20 

-— V 7 17 - ) UC ZjU 


p 






HE 389 

AAA_J OOi^ 


I 


q 023 4- 01 6 


1 33 30 + 1 2 

A.OO.OVA _1_ \J • 


K 47 + 13 

O . 1 _l_ U.AO 


1 4 q« + 1 3 

A1.C7U _1_ U.IO 


1 10 + 018 

1.1U _1_ U.IO 


2003-09-20 

Zj u uo u c Zj yj 


p 






HE 696 

AAA_J Uv/U 


I 


448 4- 001 

U . t:t:0 _1_ U.UU1 


3K7 22 + 1 8 

Jul _1_ U.A.O 


2 72 + 08 

Zj . 1 Zj _1_ U . UO 


1 2 50 + 08 

IZj.OU _l_ U.UO 


70 + 1 2 

U. 1 U _l_ U.lZi 


2003-09-20 


p 


yesi^aj 


P02 


HE 935 

AAA_J CO U 


x 


026 4- 025 

U.UiU _1_ U -Uiu 


247 44 + 21 

Ai'-t 1 . 'A'A _1_ U.iiA 


00 + 20 

U.UU _l_ U.ZjU 


q 21 + 20 

C . Zj 1 _l_ U.ZjU 


40 + 28 

U.rtU _1_ U.ZjO 


2003-09-20 

ZjUUO U C ZjU 


p 


yes(c) 


P02 




2 


3 116 4- 025 

o.iiu _i_ vj.yjAxj 


1 qq 45 + 21 

1UJ . t:0 _1_ U.iiA 


a 70 + 30 

O . 1 U _l_ U . OU 


1 7 1 « + 30 

1 ( . 1U _l_ U.OU 




2003-09-20 

ZjUUO UC ZjU 


p 


no (a) 




HII 102 


I 


3 Koq 4. 009 


21 3 29 + 1 4 

Zj A O . Zj C7 _1_ U.A_rt 


3 07 + 1 

O.Ul _1_ U.AU 


11 72 + 1 

11. 1 Z _l_ U. 1U 


1 05 + 14 

l.UO _l_ U.lrt 


2003-09-20 

ZiUUO UC ZjU 


p 


yes^aj 


B97 




2 


q qc;q 4. q 027 


240 21 +0 16 

£i'-t\J . Zj A. _j_ U.A.U 


K 75 + 10 

O.IO _l_ U.AU 


1 4 40 + 1 

Irl.rtU _1_ U.1U 


1 05 + 14 

l.UO _l_ U.lrt 


2003-09-20 

ZiUUO UC ZjU 


p 




1RXS T034423 3+281 224 


x 


425 4- 002 

U . 'tZiU _1_ U.UUi 


202 20 + 10 

Zj U Zj . Zj V 7 _1_ U.A.U 


1 13 + 010 

A..AO _l_ U.AU 


a «2 + 1 

O.UZj _1_ U.1U 


4a + 22 

U.rtO _1_ U . Z/Zj 


2002-1 1-1 7 

ZjUUZ 11 1 1 


p 


yesi^aj 


WDS 




2 


5 711 +0 006 

• 1 11 _i_ vj.uuu 


313 30 + 12 

O A O . OVA _1_ \J • A_ Zi 


7.11 ± 0.11 


14 27 + 1 1 

A Tt . Zj 1 _J_ U • 1 A 




2004-10-05 

— 1 / 17 1 1U UO 


p 




HII 571 


1 


3 903 -1- 005 

• 7 . < 7\ f 1 ) 1 U.U U U 


66 10 + 08 

' ' v7 - 1 V ' 1_ V 7 - V ' 17 


3 84 + 08 

O . L' 1 _1 U . UO 


1 3 07 + 08 

AO.UC _J_ U.UO 


95 + 1 2 

u.co _j_ u.izj 


2003-12-10 


p 


/ \ 
yes^aj 


B97 M9 


HII 1101 


1 


q 167 4- 016 

C.1UI _1_ U t u 1 u 


1 04 93 + 1 2 

Au^i co _i_ (j. az 


5 70 + 09 

O. 1 U _L_ U.UC 


1 4 46 + 09 

Irr.^TU _L_ U.UC 


40 + 16 

y f . 1 y f 1 u.iu 


2003-12-10 


p 






HII 1182 


1 


1 113 4-0 009 

ltllU _1_ U-UvM 


219.69 ± 0.26 


4 54 + 1 9 

Tt.U^t _1_ U.AC 


1 3 48 + 1 9 

A O.^O _L_ U.1C 




2003-12-10 


p 




B97 


HII 2106 


1 


240 4- 010 

v / . — t y t i_ v / . v / 1 v / 


31.09 ± 0.59 


1.71 ± 0.12 


11 29 + 12 

AI.ZjC _J_ U.IZj 




2003-12-10 


p 


yes (a) 


B97 


RX T0348 9+0110 


1 


047 4- 007 


41 50 + 3 64 


00 + 05 

U.UU _1 U.UO 


q 02 + 06 

C.UZ/ _1 U.UU 




2003-12-10 


p 


yes (a) 




HII 2278 


1 


331 + 005 


1 79 20 + 32 

A_ 1 C/.ZjVA —1— VA»OZj 


03 + 02 

U.UO 1 U.UZj 


q 57 + 03 

C.OI _L_ U.UO 




2003-12-10 


p 


yes (a) 


B97 


HII 2881 


1 


099 -1- 005 


335 73 + 1 20 

OOO. 1 O _1_ A » ZiO 


26 + 09 

U.ZjU 1 U.UC 


9 94 + 09 

C . C^t _1_ U.UC 




2003-12-10 


p 


yes (a) 


B97 


HD 285281 


1 


770 4- 001 


1 88 34 + 05 

1 L7 ',7 . v7 1 1_ V7 - '7 7 


1 20 + 10 

A_.ZjU -J-^ U.AU 


9 12 + 010 

C ■ A Zj _J_ U.1U 


66 + 14 

U.UU _l U.l^ 


2002-02-01 


p 


yes (a) 


KL98 


HD 284135 


1 


3«7 4- 002 

v / . • > y f 1 i_ v ' - y j y y ^ 


253.23 ± 0.21 


12 + 01 

U ■ A Zj _J_ U.UA. 


8 58 + 02 

O.OO _J_ U.UZj 




2002-01-31 

ZjUUZj UI Ol 


p 


yes (a) 


WDS 


HD 281691 


1 


6 768 + 014 


1 38 91 +0 13 

lUUi C A_ _1_ VA • A_ O 


1 90 + 05 

A_ . CU _J_ U.UO 


10 30 + 06 

AU.OU _1 U.UU 


74 + 08 

U. 1 r 1 U.UO 


2002-11-18 


p 


yes (a) 


KL98 


HD 26182 


1 


818 4- 002 


1 75 1 1 + 1 1 

A_ 1 O . A A_ _i_ VA • A_ A 


92 + 08 

U.CZ/ 1 U.UO 


9 09 + 08 

C.UC I U.UO 


0.42 ± 0.11 


2003-12-10 


p 


yes(c) 


WDS 


HD 284266 

aaa_/ ^o^ijUU 




t;cq 4. 006 

U.JUC _1_ u.uuu 


3K« q2 + 1 1 

OOVA.CZi _1_ VA.A.A 


1 qo + 1 

A_ . CU _l_ U.AU 


1 fifi + 1 

A U.UU _l_ U.1U 


aq + 41 

U.OC _l_ U.rtl 


2002-01 -31 

ZjUUZjUI OI 


p 


yes (a) 


KL98 


HD 26990 


1 


123 4- 004 

u. ii(U i_ v / - y ) v '1 


1 63 56 + 1 40 

A_ UO . O U _1_ A • rtVA 


38 + 20 

U.OO _1 U.ZjU 


6 81 + 20 

u.oi _i u.z.u 


67 + 36 

U.U 1 _l U.OU 


2003-12-10 


p 


yes (a) 




vB 49 


x 


2 1 3q 4- 01 7 

ii.ioi) _i_ y> .ui 1 


2^R 86 + 1 6 

Z/OU . OVJ _1_ U.A.VA 


4 «0 + 1 4 

rt.UU _1_ U.A'i 


11 40 + 1 4 




2003-01 -1 2 

ZjUUOUI IZi 


p 


yes(c) 




vB 52 


1 


1 115 + 002 

itiiu _i_ u-uUij 


236 40 + 18 

^ • >y f . 1 v 7 i_ V / - 1 c 7 


2 73 + 06 

Zj . f O I U.UU 


9 10 + 06 

C.1U _J_ U.UU 


85 + 08 

U.OO 1 U.UO 


2003-01-12 


p 


yes (a) 


P98 


vB 176 


I 


227 + 003 


307 06 + 38 

JUI . UVJ _1_ U.JO 


28 + 09 

U.ZjO _l_ U.UC 


7 K7 + 09 

1 .UI _l_ U.UC 




2003-1 2-1 

ZjUUO IZi AU 


p 


yes (a) 


Hip 


vB 66 


I 


q 701 + 023 

3. ( OA _l_ V / . U Zj • > 


248 88 + 1 1 

ZIO . (JO _l_ U.A.A 


1 75 + 1 

AU . 1 O _l_ U.AU 


1 « qi +010 

1U.C1 _l_ U.1U 




2002-1 1-1 7 

ZjVV'Z.111 i 


p 




vB 91 


I 


1 33 + 002 


1 72 qa + 2 79 

A 1 Z. CO _l_ Zi . 1 C 


37 + 1 4 

U.Ol _1_ U.A'i 


7.72 ± 0.14 




2003-1 2-1 

Zj U U O lii AU 


p 


yesi^aj 


WDS 


vB 96 


1 


171 + 003 

u » a j a _i_ y_> * \_f yj * j 


264 05 + 78 

V 7 1 . V 7 7 1_ V 7 - 1 O 


36 + 1 

U.OU 1 U.AU 


7.41 ± 0.10 




2003-12-10 


p 


yes(p) 


P98 


RX T0434 34-0226 

a V../V t.j u ao a . o u z. zaa 


I 


1 340 + 022 


271 7fi + 30 

Zj / A . / U _1_ U.OU 


2 3a + 05 

Zj . OO _l_ U.UO 


11 qq + 06 

ll.CC _1_ U.UU 


0.77 ±0.11 


2003-01 -1 2 

ZiUUOUl 1-j 


p 


yes(c) 




vB 106 


1 


7.230 ± 0.012 


76 50 + 44 

/ V7 - 7 V7 1_ U - A'A 


9 50 + 30 

C.UU _1 U.OU 


1 5 94 + 30 

1 O . C^t _J_ U.OU 




2003-12-10 


p 


no (a) 




HD 282346 

A A 1-7 Zj O Z > 1U 


I 


4fii + 001 

W.WA _1_ U.UuA 


272 14 + 18 

Zj 1 Zj. A.^ _l_ U. A.O 


1 1 3 + 04 

A. . A O _l_ U . Urt 


a qi + 04 

O.C1 _l_ U.Urt 


fji +017 

U.UA _l_ U.ll 


2002-1 1 -1 8 

Z. 1 7 W Z 1 1 IO 


p 


yesi^aj 


Hip 


vB 142 


1 


6 070 + 013 


123.82 ± 0.16 


11 30 + 20 

A A_ . OU _J_ U.ZjU 


1 8 04 + 20 

1 17 • V7 1 1 U.Z/U 




2002-11-17 


p 


no (a) 




1RXS T051111 14-281353 


I 


4qt^ + 001 


211 51+010 

ZiAA_.UA. f U.A.U 


3q + 04 

U.OC _l_ U . Urt 


a 77 + 05 

O. 1 1 _1_ U.UO 


fiq + 07 

U.UC _l_ u.u 1 


2002-02-28 

Z 1,7 1 'j_ \ I Z Z (7 


p 


yesi^aj 




HD 36869 


1 


8 230 + 014 


1 52 30 + 1 2 

A_ UZj . OU _1_ U • A_ Zj 


3 1 + 35 

O.AU _1 U.OO 


9 95 + 35 

O . CO _l U.OO 


93 + 35 

U.CO _l U.OO 


2003-01-14 


p 


yes (a) 






2 


a 043 + 0016 

. yj'-t'j _i_ u.uiu 


24q 72 + 20 

Z.rrC . 1 Z. 1 U.Z.U 


7 sq + 1 5 

1 . OC _l_ U.AO 


1 4 44 + 1 5 

Irt.rrrt I U.IO 




2003-01 -1 4 

ZiUUOUl Irt 


p 




HD 61994 


1 


5 210 + 008 


77 00 + 08 

t I .UU _1_ U.UO 


7 32 + 13 

1 .ozj _i_ u.ao 


1 2 67 + 1 3 

A Zj . U t _L_ U.IO 


28 + 30 

U.Z.U 1 U.OU 


2002-11-18 


p 


nofc) 




HD 69076 

IIJ-" UCU ( U 


I 


1 232 + 005 

A . Ztjil _1_ U.vUU 


101 06 + 1 1 

1U1.UU _l_ U.A.A 


3 qi + 05 

O.CA. _1_ U.UO 


1 38 + 05 

1U.OO _l_ U.UO 


1 04 + 21 

l.U^r 1 U.ZjA 


2002-1 1-18 

Z 17 \7 Z A 1 IO 


p 


> < 
yesi^aj 




HD 71974 


1 


383 + 014 

V / . - >' " - > 1_ u • VA A rt 


87 34 + 63 

O 1 - ■ 7 1 1_ U.UO 


42 + 05 

U.^tZj _1 U.UO 


6 45 + 06 

V ' . 1 7 1 U.UU 


35 + 21 

u.oo _i u.z/A 


2002-03-03 

— ' J 17 -— UO UO 


p 


yes(c) 


S99 


HD 72760 

111/ 1 — 1 UU 


x 


qfi4 + 007 


21 5 08 + 38 

Z. A O . UO _l_ U.OO 


4 04 + 01 

rt.Ort _1_ U.UA. 


1 28 + 02 

IU.ZjO _l_ U.UZj 


1 01 + 04 

A.U A _1_ U.Urt 


2002-1 1 -1 6 

Z1UUZ1 1 1 1U 


p 


yes(c) 




HD 77407 

lliv 1 1 t:U i 


1 


1 659 + 004 


353 36 + 04 

■ 7 7 ■ 7 . '7 y ' 1_ V7 - V 7 1 


2 00 + 1 

Zj . UU 1 U.AU 


7 60 + 1 

1 . UU _J_ U.1U 


90 + 14 

U.CU _L_ U.l^ 


2002-01-31 


p 


yes(a) 


M04 


HD 78899 


I 


a 1 74 + 01 3 


7K 7K + 1 2 

1 'J . 1 U _l_ U. A_Zj 


3 3« + 08 

O . OU _l_ U . UO 


q 1 7 + 08 

C . 1 1 _1_ U.UO 


7q + 1 3 

U. 1 C _1_ U.IO 


2003-1 2-09 

ZjUUO A Zj UC 


p 




HD 91 962 

A A 1 ' ClCUZi 


I 


«42 + 003 

yj . Ortij _i_ yJ •yjyj-J 


1 7fi 00 + 1 1 

1IU. UU _l_ U. A. A 


1 37 + 06 

A. . O 1 _l_ U.UU 


7 03 + 06 

1 . UO _1_ U.UU 


0.73 ± 0.13 


2002-03-02 

Z1UUZ1UOUZ1 


p 


yes^aj 


WDS 




2 


1 42 + 004 


5fi 1 7 + 1 7fi 

JU.l 1 _1_ A.. 1 U 


i 25 + 011 

A. . Zj O I U . A A. 


« q4 + 1 1 

U . Crt _l_ U.ll 


Q2 + 1 9 

U.CO _l_ U.1C 


2003-05-1 

Zj U U O U O AU 


p 


ycs(c) 




HD 99565 

A A A_/ C'/OUO 


I 


408 + 001 

V 7 . 1 V ' (7 > \> .y '\ > I 


6 13 + 18 

U.A.O _l_ U.A.O 


09 + 05 

U.UC _l_ U.UO 


« kk + 05 

U.OO _l_ U.UO 


42 + 06 

U.rtO _l_ U.UU 


2003-01 -1 1 

ZiUUOUl 1 1 


p 


yes (a) 


WDS 


HD 108799 

Ally AUG 1 CC 


\ 


2 070 + 006 

Z. . V 7 / W 1 V 7 - V 7 v » V ' 


33a 4fi + 09 

OOO. rtU _l_ U.UC 


1 47 + 02 

A_.rtl _1_ U.UZj 


« 30 + 03 

U.OU _l_ U.UO 


ao + 06 

U.OU _1_ U.UU 


2003-05-1 

ZjUUO U O 1U 


p 


yes (a) 


WDS 


HD 108944 

A A 1 ' AUO\/ 1 1 


\ 


1 q4i + 006 

A_.CrtA _1_ U.UUU 


345 48 + 1 8 

OrtO.rtO _l_ U.lu 


3 4q + 02 

O.rtC _1_ U.UZj 


q k« + 03 

COU _1_ U.UO 


as + 09 

U.OO _l_ U.UC 


2002-03-03 

ZjUUZi UO UO 


p 


yes (a) 




HD 1 1 21 96 

A A 1 ' A 1^1 i7U 


I 


1 ko1 + 001 

1.UU1 _1_ U.UU1 


w K2 + 09 

OU.OZj _1_ U.UC 


2 07 + 01 

Zj . U 1 _l_ U.UA. 


7 77 + 02 

1 . 1 1 _l_ U.UZj 


73 + 04 

U.IO _l_ U.Urt 


2002-02-01 

Z1UUZ1UZ1UI 


p 


yes (a) 




HD 115043 


1 


1 639 + 003 

A_ . VAO C ™t_ U * V7 W ' J 


358 61 +0 05 

■77 v ' . Vf 1 L_ U.UU 


4 87 + 08 

rt.UI 1 U.UO 


10 22 + 08 

AU.ZZ I U.UO 


86 + 12 

U.OU _L_ U.IZj 


2003-12-09 


p 


yes (a) 


L05 


HD 129333 


1 


71 7 + 009 

u» 1 a_ 1 _i_ yj *\jyj *y 


1 72 77 + 1 1 

A 1 Z. 1 1 _1_ U . -A A_ 


2 83 + 05 

Zj . OO 1 U.UO 


8 82 + 05 

O.OZ/ I U.UO 


81 + 07 

U.Ol _L_ U.U 1 


2003-01-11 


p 


yes (a) 


DM91, MI 


HD 134319 


1 


5 356 + 020 

7 . ' 7 7 V 7 1_ V ' - V7 — V7 


260 77 + 1 

z.uu. 1 t _i_ u.a_u 


4 00 + 10 

rt.UU 1 U.AU 


1 79 + 1 

AU. 1 C _J_ U.1U 


85 + 1 1 

U.OO _J_ U.ll 


2002-03-02 


p 


yes (a) 


L05 


HD 135363 


1 


0.251 ± 0.003 


121.35 ±0.46 


0.68 ±0.10 


7.34 ±0.10 


0.57 ±0.14 


2002-02-01 


P 


yes(a) 


L07 


RX J1541. 1-2656 


1 


6.261 ± 0.018 


82.05 ± 0.13 


3.13 ±0.02 


12.05 ±0.03 




2003-07-15 


P 






2 


6.250 ±0.015 


224.11 ± 0.15 


7.19±0.11 


16.11 ±0.11 




2003-07-15 


P 


no (a) 




PZ99 J161329.3-231106 


1 


1.430 ± 0.002 


91.41 ± 0.05 


2.70 ±0.05 


11.28 ±0.05 


0.87 ±0.08 


2003-05-10 


P 


yes (a) 




HD 150554 


1 


11.595 ±0.023 


183.44 ± 0.08 


3.06 ±0.10 


9.37 ±0.10 


1.03 ±0.14 


2003-05-10 


P 


yes(p) 


WDS 


HD 152555 


1 


3.819 ±0.008 


56.86 ±0.15 


3.78 ± 0.02 


10.14 ±0.03 


0.81 ±0.05 


2002-08-31 


P 


yes(a) 




HD 155902 


1 


0.062 ± 0.007 


0.28 ±6.05 


0.50 ±0.30 


6.26 ±0.30 


0.39 ± 0.42 


2003-09-21 


P 


ycs(c) 




HD 157664 


1 


0.036 ± 0.002 


118.76 ± 3.21 


0.00 ±0.10 


7.46 ±0.10 


0.31 ±0.14 


2003-05-10 


P 


yes(a) 




HD 166435 


1 


2.653 ±0.022 


273.69 ± 0.26 


10.67 ±0.20 


15.99 ±0.20 


-0.15 ±0.28 


2002-06-23 


P 


no (a) 






2 


10.376 ± 0.030 


281.28 ± 0.12 


11.90 ±0.20 


17.22 ±0.20 


0.54 ±0.22 


2002-06-23 


P 


no(a) 






3 


9.496 ± 0.020 


183.40 ± 0.19 


11.48 ±0.20 


16.80 ±0.20 


0.18 ±0.28 


2002-06-23 


P 


no (a) 






4 


3.293 ± 0.009 


239.04 ± 0.45 


13.50 ±0.30 


18.82 ±0.30 




2002-08-30 


P 


no (a) 
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TABLE 10 

Color Companions and Their Chance Alignment 
Probabilities 



Companion to 


Nc 


P 


AK S 






CAP 






(arcsec) 


(mag) 


(mag) 




(%) 


HD 1 5526 


X 


0.0770 


0.00 


8.01 ± 0. 


10 


0.8 


HD 1 55902 


I 


0620 


0.50 


5.73 ± 0. 


10 


14.3 




I 


0.1670 


0.78 


6.28 i 0. 


10 


1.5 




[ 


2.0790 


3 97 


9.70 i 0. 


08 


9 


vB 1 


I 


2.4700 


2.63 


8.62 ± 0. 


10 


0.9 


HD 245567 


1 


0.3480 


1.79 


9.38 ±0. 


10 


0.8 


HD 26182 


1 


0.8180 


0.92 


8.71 ±0. 


10 


2.3 


vB 49 


1 


2.1390 


4.56 


11.36 ±0. 


10 


1.3 


HD 71974 


1 


0.3830 


0.42 


5.89 ±0. 


10 


1.7 


HD 72760 


1 


0.9640 


4.84 


10.26 ±0. 


10 


0.9 


HD 91962 


2 


0.1420 


1.25 


6.64 ±0. 


10 


0.9 


HE 935 


1 


0.0260 


0.30 


8.76 ± 0. 


10 


6.7 


RX J0329.1+0118 


1 


3.7610 


3.62 


12.82 ±0. 


10 


0.7 


RX J0354.4+0535 


2 


0.2050 


2.10 


10.77 ±0. 


10 


1.1 


RX J0434.3+0226 


1 


1.3400 


2.38 


11.88 ±0. 


10 


1.2 


RX J2313.0+2345 


1 


1.4060 


1.79 


10.41 ±0. 


10 


1.4 


ScoPMS 27 


1 


0.0790 


0.60 


8.64 ±0. 


10 


0.7 



NOTE. — Most columns headings are as for Ta bles QQ and \p\ CAP is the chance 
alignment probability from Equation^in g |5.2.2| 
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TABLE 11 

New and Confirmed Secondaries 





M K 


Projected Separation 


M 2 








Companion 


(mag) 


(arcsecj 


(AU) 


(M e ) 


<3 


Sample a 


Ref. 


NT) 994.°s7^R 

L1LJ £ZfO I O-Lj 


4 1 o _|_ n 99 


1.27 


62.13 


0.84 


0.98 


A Don 


WDS 


HD 9472R 


q 94 + n 1 1 


2.79 


92.17 


0.11 


0.10 






RP tai QV-U 1 8R 


^ 4fi -i- n 97 


1.69 


108.22 


u.uo 


n Q7 


A Don 


WD^ 


un 1 R 

L1LJ _LOUOJ-J_> 


701 _i_ n 1 9 


0.72 


18.64 


0.19 


0.20 


A Don 




TTD 1 ^9fiR 

L1LJ 1 JOiUiJ 




0.08 


8.16 


0.90 


1.00 


AD 




109^993 c »4- < ^79Q1zLR 


A K9 4- 1 HQ 


0.64 


108.29 


71 


0.66 


ADon 




9RF T09 c > r i-UzL7zLR 


q on _|_ 1 HQ 


2.13 


106.55 


0.91 


0.99 


ADon 


WDS 


RY T09^& 4.4-9Q4.7R 


C 1 C _|_ 1 1Q 


0.09 


8.60 


0.59 


0.78 


AD 




un 1&Q40R 


a r\x. -u n ns 

t.UtJ ^ u.uo 


0.17 


5.68 


86 


0.83 




Hip 


vB IB 


o.*±o ^ u.uo 


2.47 


106.46 


U. UU 


n 47 




WD^ 
VV i-JJ 


RF ^7'^R 


on 4- n i fi 

O. ^U ^ U. 1U 


2.08 


395.39 


0.10 


0.08 




MH04 

1V111U4: 


UT? fiQfiR 


« 1 1 ±n ik 

U. J. -L ^ U. J-U 


0.45 


85.12 


39 


0.38 




P02 


RY T0 4; i9Q 1-U011&R 


7 qo _|_ 1 HQ 
( .oz ^ 1 .uy 


3.76 


376.10 


n 1 1 
u. r 1 


12 


A Don 


MR 04. 


RF Q'^R 


09 _i_ n 9zt 


0.03 


4.94 


1.20 


1.00 




P02 


HTT 1 09R 


f: 1 n 4- n 14. 

U. ±U ^ U . -L'-t 


3.60 


478.67 


43 


0.39 


ADon 


B97 


1RYQ Tfl'344.9'3 ^-U9£1 994.R 


c 1 7 4. n 9n 

0. ± t ^ u. ^u 


0.43 


20.83 


57 


0.71 


AD 


WDS 


HTT 571 R 


7 4. A 1 Q 
1 .rttJ _1_ U.-LiJ 


3.90 


519.10 


0.20 


0.20 




R97 MQ9 


TTTT 1 1 89R 


7 4. n 91 

i . ou ^ u. ^ ± 


1.11 


148.03 


0.15 


0.14 




R97 


RTT 91 OfiR 


O.U / ^ U. J-vJ 


0.24 


31.92 


51 


0.59 




R97 


rx T0348 q+oi i or 


4 09 4- 1 HQ 


0.05 


4.70 


0.87 


1.00 






RTT 997SR 

1_±±±. __ 1 OJ_> 


Q QC; 4. A 1 A 


0.33 


44.02 


0.89 


0.99 




R97 


HTT 2881 R 


4 39 4- 1 3 


0.10 


13.17 


0.79 


0.92 




R97 


RY T0^4 A-Ufl^^^R 

iwv J UOO i ±.'±T^UOOO-D 


koo-Li in 


0.21 


20.50 


n 48 


U. UU 


AD 






c; «7 4. a 9fi 


0.77 


37.73 


0.19 


0.42 


ADon 


KL98 


un 9SA1 ^R 


9 qc: _i_ 1 nq 
^ . 00 ^ ± . uy 


0.37 


51.38 


98 


0.93 


AD 


WDS 


un OKI fiQ1 R 


a c;7 4- 1 HQ 


6.77 


947.52 


68 


0.60 


ADon 


KL98 


un O^l «9R 

1_±_L/ U J. O J_> 


4 0Q4-1 HQ 
i.uy ^ ±.uy 


0.82 


81.80 


0.79 


0.72 


ADon 


WDS 


un 9&4.9fifiR 


A QQ 4- 1 flQ 

4i.yo ^ J. .uy 


0.57 


79.66 


n 4*3 


u.oo 


A Don 


KL98 


un 9^QQnR 


4 aq 4. a 94. 


0.12 


4.31 


0.84 


0.91 






vB 49B 


7 fin + fl 1 ^ 

1 . UU _1_ KJ . A-iJ 


2.14 


122.99 


0.22 


0.19 






V J_> tj^ J_> 


r: qa 4- fl 07 


1.12 


49.95 


0.49 


0.45 




P98 


vR 1 7fiR 


4 9fi 4- 1 

'rl.ZiU _i_ u. j_u 


0.23 


10.90 


0.80 


0.99 




Hip 


vB 91B 


J. A1 4- 1 A 


0.13 


6.10 


n 7fi 

U.iU 


n Q4 
u. 




WD^ 


vR QfiR 

V J_> C/UJJ 


4.12 i 0.11 


0.17 


7.76 


0.84 


0.97 




P98 


RX J0434.3+0226B 


5 Qfi J_ A QQ 
'-J * *J _1_ VJ . t-J'-J 


1.34 


215.74 


0.38 


0.39 


ADon 




HD 282346B 


4.65 ± 0.43 


0.46 


32.73 
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This work 



NOTE. — Surveys arc listed in approximate chronological order. For > 0.3 Gyr-old stars, masses arc es timated fro m the models of |Baraffe et al.| 11993 )- For solar 

neighborhood- aged ( 5 Gyr) stars , mass estimates follow the spectral type— mass correspondence from |Cox| |200C| ) . The median outer projected separation p ut 

is obtained as the product of the median sample distance d and the half-width of the FOV (i.e., the OWA) of the imager used in the survey. The sensitivity of 

each survey, in units of the limiting companion mass, corresponds to the median sensitivity to substcllar companions at the widest probed separations, generally 

well outside the contrast-limited regime. Where this sensitivity was not explicitly stated, it was estimated based on the published survey depth and on substcllar 

evolutionary models from Chabrier ct al. i 2000.) and Baraffe et al. i2003). iV BD and / BD arc the number of detected brown dwarf comp anio ns and the fraction 

of survey stars with brown dwarf companions, respectively. The label in the penultimate column refers to the survey identifier in Figure |TT| The comments on 

the adopted parameters for each survey from the last column are as follows: 1. I ncluded for comp arison to the radial velocity brown dwarf desert. The median 

spectral type and stellar mass have been estimated approximately. 2. The work of |Gizjs et aX\ <2001| ) analyzes brown dwarf companions to <M0 stars within 25 pc 

in the 2MASS Second Incremental Data Release (IDR2). The outer probed separation range is likely > 10 4 AU. The detection rate has been obtained from the 

ratio of the number of bound brown dwarf co mpanions to the estim ated number of brown dwarfs within 25 pc in 2M ASS IDR2, assumin g a field mass function that 

is flat across the stellar /substcllar boundary (Metcliev^^lJ^200a). 3. The medi an spectral t y pe in the survey of IPotter et all 120021) has been estimated based 

on the sample of nearby young solar analogs cT TBSdoi^iHaf 120001 ), from which [Potter et al"| <2UU2fl borrow to form their sample. Also, the detected substellar 

binary c ompanion, HP 130948B/C, is cou nted as a single compani o n obj ect. 4. We estimate a median spectral type of G5 for the young Sun-like stars in the 

survey of INeuhauser & Guentheij f 2004ft . 5. IMcCarthy & Zuckermanj 12004|) search for substellar companions onl y to the 83 apparently single stars in their 102-star 

Keck survey. 6. The parameters of the solar analog survey of young southern associations by Chauvin et al. I 2005b) have been guessed. 7. The spectral type 

distribution of the sample targets in |Kouwenhoven et alj (2^0*51 , 120071) is approximate. A median primary mass of 2.0 MQhas been assumed. The listed brown dwarf 

companion is a candidate pending astrometric and spectroscopic confirmation. 8. The work of | Luhman et all IZOOTb] ) surveys two distinct populations of Sun-like 

stars, which have been listed separately hero based on the samples of their two Spitzer programs (PID = 34 and PID=48). The sensitivities of the two sub-surveys 

are estimated approximately. 9. The Palomar AO survey sample of Tanner et al. i2()07) contains 15 stars in Taurus and 14 stars in the Pleiades. Definitive proper 

motion associations arc available only within \" of the primaries. Here and in Figiire |l3l we h ave shown on l y the Taurus subset because only that attains sensitivity 

to substcllar objects within l" from the primaries. 10. Half of the sample observations of |Kraus et al.| 120USI ) are sensitive to companions below the 13 Aij up 

deuterium-burning mass limit, and half arc not. Therefore, wc have adopted 13 M Jup as the median sensitivity mass limit of the survey. The listed brown dwarf 

companion is a candidate pending astrometric and spectroscopic confirmation. 

a b 

The target source sample for each work. Median value for the primary stars in the survey. 
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Fig. 1. — Distribution of the sample stars as a function of effective temperature (a) and mass (b). The non-shaded histograms refer to 
the entire sample of 266 stars, whereas the shaded histograms refer to the deep and young sub-sample of 100 stars. All stars fall in the 
F5-K5 range of spectral types and the majority are between 0.7 Mq and 1.3 Mq. 
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Fig. 2. — Age distributions of the complete survey sample (non-shaded histogram) and of the deep sub-sample (shaded histogram). 
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Fig. 3. — Heliocentric distance (a) and proper motion (b) distributions of surveyed stars in the complete sample (non-shaded histograms) 
and in the deep sub-sample (shaded histograms) . 
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Fig. 4. — Empirical _ft~5-band contrast limits as determined from artificial star experiments in images of the program star HD f72649 
(V = 7.5 mag), taken under good AO performance (ft:50% Strehl ratio). The solid and long-dashed curves delineate coronagraphic 
observations at Paiomar (24 min) and Keck (6 min), respectively. The short-dashed line shows the non-coronagraphic component of the 
Paiomar survey. The dotted line represents the 4<r r.m.s. deviation of counts in the PSF halo as a function of separation, normalized to an 
aperture with radius O'.'l: equal to the FWHM of the i\g-band PALAO PSF. The vertical dash-dotted line shows the edge of the occulting 
spot at Paiomar and Keck. The slight decrease in contrast in the Paiomar coronagraphic limits at >5" separations is due to an additive 
parameter used to mo del the decreasing exposure depth toward the edge of the PHARO field, because of image mis-registration among the 
different CR angles (£j 13.1. It . The contrast degradation is set to vary between mag and 0.75 mag in the 4"0— 12'.' 5 separation range. The 
bumps and spikes in the r.m.s. limits correspond to bright features in the image of HD 172649, such as the corners of the waffle pattern at 
1.0" and projected companions to the star at 2.1", 4.8", and 8.6". 
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Fig. 5. — Contrast (a) and depth (b) of the deep survey at Kg. The solid lines represent the 10%, 50% (thick), and 90% completeness 
of the combined Palomar + Keck AO survey. The median (50%) sensitivities of the Palomar (dotted line) and Keck (dashed line) surveys 
are also shown. The gradual decrease in i maging contrast and depth at Palomar between 4"— 12''5 is partially due to mis-registration of 
images taken at different CR angles (§ 13.21 1. and partially to the sometimes smaller depth of observations at 11"-12'.'5 separations because 
of a 0'.'5— 1'.'5 offset of the coronagraphic spot from the center of the PHARO array. 
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Fig. 6. — Magnitude difference AK$ vs. angular separation p for all candidate companions discovered in the deep and shallow surveys. The 
various symbols denote: "•" — astrometrically associated companions; "x" — astrometrically unassociated background stars; and for objects 
with insufficient astrometric data: "o" — companions associated based on their JKg photometry; "+" — objects with JKg photometry 
inconsistent with association; "A" — undecided objects. The encircled points show the two brown dwarf companions from the survey: 
HD 49197B (at p = C//95) and HD 203030B (at p = 11"92). Detection limits for the shallow (dashed line) and deep (solid and dotted 
lines) components of the survey are also shown. The solid line shows the median contrast AK$ of the deep survey, while the dotted lines 
delimit the 10—90 percentile region (cf. Fig. [5^.). Binaries with separations smaller than the PALAO i^g-band diffraction limit (CK'10) were 
resolved only at J-band. Correspondingly, the plotted magnitude difference for these companions is the one at J. 
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Fig. 7. — Mk s vs. J — Ks color-magnitude diagram of candidate companions for which J-band photometry was obtained. The symbols 
are the same as in Figure [6] The additional small dots denote M0-T8 dwarfs (Mk s ^ 4.5) with parallaxes from lDahn et al.l 1120021) and 
IVrba et al.l ( 120041^ . The points with errorbars represent the two confirmed brown dwarf companions from our survey. The errorbars on 
HD 203030B are representative of the photometric precision for the faintest (J > 18 mag) objects in the survey. Brighter objects typically 
had J — Ks errors < 0.3 mag, except for the large J — Kg uncertainty of HD 49197B, which is unique because of its relative faintness 
(AJ = 9.6 mag) and proximity (p = 0'.'95) to the primary. The vector in the upper right corresponds to Ay = 2 mag of visual extinction, 
equivalent to a distance of ~3 kpc, or a distance modulus of 12 mag along the galactic plane. 
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Fig. 8. — _R"s-band image of ScoPMS 214 and its candidate companions taken with NIRC2 and the Keck AO system on 5 June 2004. 
The image is the median of nine 60 s exposures. ScoPMS 214 is occulted by a partially transmissive l"0-diameter circular coronagraphic 
mask. The seven p < 12'.'5 candidate companions listed in Table[8]are pointed out with arrows. The candidate proper motion companion 
CC1 is the brightest of the seven and closest to the star. 
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Fig. 9. — Proper motion diagram for candidate companions to ScoPMS 214, spanning the At = 4.81 yr time period between the first 
and last epochs of observations, between 30 August 2002 (to = 2002.66 yr) and 23 June 2007. The dashed line denotes the expected 
relative motion of a stationary background object with respect to ScoPMS 214 between the initial epoch (marked with x) and the final 
epoch (marked with thick errorbars without a solid point). The solid points with thin errorbars denote the observed changes in the relative 
positions of candidate companions 1, 2, 3, 4, and 6. Candidate 5 was outside of the field of view of the medium (20" X 20") NIRC2 
camera during the last epoch of observations, and candidate 7 was below the detection limit during the initial epoch. Candidate 1 is 
ScoPMS 214 "B", which shares the proper motion of ScoPMS 214 during the 4.81-year time span within 3<r limits and is inconsistent with 
being a stationary background object (at the her level in declination). Candidates 2, 3, and 4 (and 7, based on observations at intermediate 
epochs) are all inconsistent with being proper motion companions to ScoPMS 214 and are consistent with being background objects. 
Candidates 5 and 6 are consistent with being either bona fide companions or unrelated background objects, i.e., their status is undecided. 
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Fig. 10. — if-band spectrum (red) of ScoPMS 214"B" ( CC1), compared to spectra of M3-M6 field dwarf (in black) and giants (in 
blue) from the IRTF Spectral Library (Cushirig ct al. 2005; Rayncr et al. 2008), smoothed to the same R 1200 resolution. Dominant 
absorption features by Na I at 2.21 /im (doublet) and 2.34 fim (doublet), Ca I at 2.26 ^m (triplet), and CO bandheads at A > 2.29 fim 
are identified. The comparison dwarf spectra are of Gl 388 (M3V), Gl 213 (M4V), Gl 51 (M5V), and Gl 406 (M6V), and the giants are 
HD 28487 (M3.5III), HD 214665 (M4+III), HD 175865 (M5III), and HD 196610 (M6III). 
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Fig, 11. — HR diagram of the candidate binary ScoPMS 214A"B" with evolutionary models for 0.05-1.4 Mq objects from lBaraffe et al.l 
(1998). The continuous lines are isochrones and the dashed lines are evolutionary tracks at constant mass. The thick (1 Gyr) isochrone 
approximates the main sequence, and the thick evolutionary track corresponds to the minimum hydrogen-burning mass. The positions 
of ScoPMS 214A and "B" under the assumption of equidistance and membership in Upper Scorpius (the "young" ScoPMS 214"B" 
scenario; § 16. 3.21 1 are shown with solid points with errorbars. The shaded region represents the range of effective temperature allowed for 
ScoPMS 214 "B" if it were an unassociated field-aged (1-10 Gyr) M dwarf. Since in the "young" scenario the candidate binary components 
do not lie on the same theoretical isochrone, ScoPMS 214"B" is probably not a member of Upper Scorpius. Instead, it is most likely a 
foreground field M dwarf. 
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Fig. 12. — Probability density distribution P(fJ.\x = 2) for the expected substellar companion detection rate in our survey per 100 stars, 
given x = 2 detections. The curve is a Gamma distribution (Eqn. 0, with a peak at fi = /^ml = 2, but a mean value of (fi) = x + 1 = 3. 
The minimal 2a (95.4%) confidence interval on (fi), 0.3 < (fi) < 6.5 is indicated by the shaded region under the curve. 
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Fig. 13. — Substellar companion detection rates of the published direct imaging surveys listed in Table IT2l Circular symbols denote 
surveys with at least one detection; crosses denote surveys with no detections. The filled circle denotes the present work. The size of the 
circular symbols is proportional to the survey detection rate prior to corrections for survey incompleteness. Black symbols denote the least 
sensitive surveys, with > 30A/j up median companion mass sensitivity in the background-limited regime. Blue symbols denote surveys with 
median companion sensitivities between 13-30 Mj up . Red symbols mark surveys with the highest sensitivity, < 13Mj up . The survey labels 
are as listed in the penultimate column of Table 1121 The locus delimited by a dotted line contains only surveys with non-zero detections, 
with detection rates ranging from 0.5-5%. All surveys outside of this region have detection rates < 0.6%. 
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Fig. 14. — Mass ratio distribution for the 30 <0.5 Gyr-old binaries in our AD30 sample (see § 19. Il l in terms of log(dN/d log q) (a) and 
dN/dq (b). The dotted histogram traces the observed data, while the solid histogram delineates the incompleteness-corrected data. Further 
incompleteness due to bias against near-equal binary systems exists in the highest mass ratio bin, but has not been taken into account 
in the present incompleteness correction. The long-dashed line is a power-law (PL) fit to the data, dN /d log q oc with an index of 

(3 = —0.39 ± 0.36 (lcr limit). The short-dashed line represents the log-normal MF of field objects from Chabrier (2003, C03) in units of 
Mq, normalized to the incompleteness-corrected data. We note that because the primary masses for stars in our sample are Ril Mq, then 
q = M2/M1 PS M2/MQ. The log-normal field MF peaks at fi = 0.08Mq and has a width of a = 0.69 (in logarith mic mass units). The 
Salpet er index in these units is a = —2.35. The solid points in pa nel (b) are the incompleteness-corrected data from Duqucnnoy & Mayor 
(1991), normalized to our data. The Duqucnnoy & Mayor (1991) data have been offset slightly to the right from ours for clarity. 
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Fig. 15. — Projected physical separations probed in the deep sample survey. The vertical dashed lines delimit the region, 22-1262 AU, in 
which each 1 AU-widc projected separation interval was probed around at least one third of the stars in our deep sample. The geometrical 
incompleteness factor for this region is 1.40. That is, 1/1.40 = 71.4% of all companions in the 22—1262 AU projected separation range 
should have in principle been detected, had their visibility not been limited by contrast. 
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Fig. 16. — (a) Observational (circles) and total (triangles) completeness of the deep survey as a function of companion mass. The 
observational completeness at a given mass is the fraction of companions of that mass that would be detectable within a projected separation 
of 22—1262 AU from all sample stars (§ IA,3,2| |. The total co mpleteness is defined similarly, but for a 28-1590 AU range of semi-major axes, 
and after consideration of orbital incompleteness IA.3.31 . Both sets of completeness fractions are calculated assuming a logarithmically 
fiat distribution of companion semi- major axes a (§ I A ,2ft . The horizontal lines delimit the maximum possible observational (continuous 
line) and orbital (long-dashed line) completeness at any given mass ove r these AU range s. Our defini tion of the orbital completeness 
coincides with the "single visit obscurational completeness" (SVOC; see § IA. 3.31 1 defined by Brown (2004). The vertical dotted lines mark 
the deuterium- (D) and hydrogen- (H) burning mass limits, (b) Same as Figure [T5l but for the expected semi-major axes (rather than 
projected separations) of substellar companions and for a range of companion masses. The dotted lines are labeled with substellar masses 
in units of Mq/100. The solid curve delineates the geometrical completeness limit and the long-dashed curve, the SVOC limit (cf. panel 
a). The vertical short-dashed lines have been adjusted from their positions in Figure [T5l to correspond to the expected range of semi-major 
axes, 28-1590 AU, corresponding to the 22-1262 AU projected separations probed by the survey. 



